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ABSTRACT 

Analysis and design of linear feedback control systems by application 
of the basic Mitrovic method is well established. However, little has been 
done in applying the basic Mitrovic method to nonlinear systems. The objec- 
tive of this work was twofold. First, Mitrovic’s methced as applied to 
systems with a single gain-variable nonlinearity predicts limit ecvcles with 
results comparable to existing analysis methods, Secondly, an analog comput - 
er study was made on 42 system with two gain-varilable nonlinearities and the 
results interpreted by Mitrovic's method, 

The analog computer results were obtained by a Domner Scientific Corpor: 
tion analog computer, Model 3100, at the U. S. Naval Postgraduate School, 
Monterey, California. The writer wishes to express his appreciation for the 
assistance and encouragement given him in this endeavor by Professor Ceorge 


Thaler of the U. S. Naval Postgraduate School. 
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PART I 
PREDICTION OF LIMIT CYCLES BY MITROVIC'S METHOD FOR 
CONTROL SYSTEMS WITH ONE GAIN VARIABLE NONLINEARITY 
Introduction. 

Mitrovic's method is easily applied to nonlinear control systems if 
the nonlinearity can be described by an input-output characteristic in 
which the output is an odd single-valued function of the input. Such 4 
nonlinearity ene be expressed as an instantaneous variable gain 
with zero phase shift. This property allows the nonlinear effect to be 
expressed as a variable coefficient of the system characteristic equation. 
For the nonlinearity in the error channel, such a system can only have limit 
cycles if the linear system is unstable since the nonlinearily contributes 
no phase shift. 

Mitrovic's method determines absolute stability for a linear system, 
which has fixed coefficients in the characteristic equation, by assuming 
two of the coefficients to be variable. The method derives equations 
for two coefficients with 4% and Wr of the S-plane as parameters. 
With the JW axis (¥ = 0) as the S-plane mapping contour, the curve of 
one coefficient versus the other coefficient is plotted using different 
values of Wy. If this curve encircles a point defined as the working 
point or M-point in a clockwise direction for increasing Wry, the 
system is stable. The working point is determined by the fixed values 
of the two coefficients which were assumed variable. 

For a nonlinearity, as described above, one or more coefficients of 
the characteristic equation are variable. This manifests itself on the 


variable coefficient plane as a locus of possible working point locations. 





te, = Ce 


Thus a limit cycle for the nonlinear system exists for the conditions 

at the intersection of the working point locus and the stability (¥% = 0) 
curve. A typical Mitrovic plot is shown in Figure Il. The frequency of 
the limit cycle is read from the stability curve at the intersection. 

The amplitude is determined approximately from the gain of the non- 


linearity obtained from the working point locus. 
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Figure I 


Coefficient Plane of Mitrovic's Method 
The amplitude determined in this manner will necessarily be less 
than the actual amplitude. The gain of the nonlinearity 'is defined as 
the ratio of the instantaneous output to the instantaneous input and the 


actual amplitude could be determined accurately if during a cycle the 





minimum/maximum instantaneous gain could be determined. However, the 
non-linear gain determined from the Mitrovic plot is a "psuedo - average" 
gain over a complete cycle, since this "psuedo - average’ gain is deter- 
mined from an "average" M-point position over a cycle. 

The shape of the stability curve on a coefficient plane depends on 
the two coefficients chosen and the order of the characteristic equation. 
If the coefficients of the S° and st are chosen variable, then the stability 
curve for a third order system is a straight line on the Bo versus & 
plane, for a fourth order system the stability curve is a parabola, and 
so on. If different coefficients are chosen variable, then the stability 
curve has a different shape on the corresponding coefficient plane. 

The working point locus depends also on the linear system and specti- 
fically the number of zeros in the loop transfer function. For no zeros 
in the loop transfer function and the nonlinearity in the forward path, 
the working point locus is as straight line segment parallel to the Bo 
axis since the variable nonlinear gain appears only in the constant co- 
efficient of the characteristic equation, For a linear system with one 
zero in the loop transfer function, the working point locus is again a 
straight line segment with a specific slope on the RB. versus iE, plane. 
Other systems can be conceived with different placements of the nonlinear- 
ity but the basic working point locus is always a straight line segment 
on the appropriate coefficients plane. 

The use of Mitrovic's method for a nonlinear system for which the 
linear loop transfer function has two or more zeros becomes more compli- 
cated, [n this case the nonlinear gain appears in three or more ¢o- 


efficients of the characteristic equation, With only two assumed variable 





coefficients, the stability curve becomes a family of curves with the 
nonlinear gain a parameter. To determine the limit cycles for such a 
system the intersection of the straight line working point locus and 
the family of stability curves must be found. 
The various systems chosen to analyze in the chapters of Part One 

are completely arbitrary. The author had no specific system or specific 
device in mind when the block diagram was drawn. Consequently, all signals 
in the various systems are referred to in volts. The following guidelines 
were followed in order to simplify the analysis without disrupting the 
arbitrary nature of the system: 

1. The poles and zeros of the linear loop transfer function 
were in general chosen one octave apart. This was done in order to 
sketch the actual Bode magnitude curve in the describing function anal- 
ysis from the asymptotic plot by applying the usual corrections at the 
corner frequency and two octaves either side of the corner frequency. 

2. The magnitudes of the poles, zeros, and gain constants 
were chosen to eliminate magnitude and time scaling procedures in the 


analog computer simulations of the system. 





CHAPTER I 
SATURATION IN ERROR CHANNEL 
I-A Loop Transfer Function with No Zeros. 


' The block diagram of the system chosen to analyze is given in Figure 


bd 


I-A. 
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Figure I-A 


Saturation in Error Channel 
Loop Transfer Function with No Zeros 
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I-A-l. Mitrovic's Method, 


Assigning the variable gain N to the nonlinearity, the character- 


istic equation becomes, 


674354 25+ ION =O ; Ceaeee 


Assuming the last two coefficients variable the characteristic 


equation is rewritten, 


Sor S) s* + B,S+B, =O . (I-A-1,.2) 





From Appendix A, Mitrovic's equations for Bo and By are, 


Bp = -[30,7 (4) + Wn? hol) ] (I-A-1.3) 
and 


B,= 32wn0,(4) + wn Da (S) : (1-A-1.4) 


Substitution of the values of the @ functions Appendix A for ¥ = 


gives the parametric equations of the stability curve as 

= Say (T-A-1.5) 
and 

6) CORA (1-A-1.6) 


From comparison of equations I-A-1.1 and I-A-1.2, it is seen that 


the equations describing the M-point locus are 


cw ION (I-A-1.7) 


and 
b, _ a. (I-A-1.8) 


Thus the set of four parametric equations Tf-A-1.5 through I-A-1.8 
describe the two curves whose intersection is desired for limit cycle 
determination, 

At the intersection, Equation I-A-1.5 equals Equation [-~A-1.7 and 


Equation I-A-1.6 equals Equation I-A-1.8. Thus, 





and 
a 
Solution of these two simultaneous equations yields, 
, Lon = I, Al4 rad/sec ; (I-A-1.11) 


and 


Nave = 0.6, ‘ (TA 1 


Since the nonlinear gain N is defined as the ratio of the output 


‘to the input, 


Sa 


E may = ie (I-A-1,13) 


Thus the frequency of the limit cycle is given by equation I-A-1.11 


and the approximate amplitude by I-A-1.13 or 


Ewe One) Nolts we  (I-A-1,14) 
es The above solution.is presented graphically in Figure I-A-1.1._ = | 
wal ag STABILITY CURVE (450) 
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© N=0.6 ae LIN IT CYGIEE 
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Figure [-A-1.1 


Graphical Solution for Saturation in Error Channel, Loop Transfer 
Function with no Zeros 


7 





The limit cycle is a stable limit cycle as determined by the follow- 
ing reasoning. If the input signal, 0; , is large enough such that N 
is small enough to place the instantaneous M-point at say My, in Figure 
I-A-l.1. The system is instantaneously stable and the oscillations de- 
crease, As the oscillations decrease, N increases moving the instantan- 
eous M-point above the stability curve. Here the system is unstable and 
the oscillations increase thereby decreasing N and moving the M-point back 
into the stable region. Eventually a dynamic equilibrium will occur in 
which the "pseudo average" M-point lies on the intersection of the two 
curves. If the input signal is small enough, N will be large enough 
such that the instantaneous M-point lies at say Mo at which point the 
system is unstable. By similar reasoning as above, the “pseudo - average" 


M-point will again eventually end up at the intersection. 





I-A-2, Root Locus Method, 

The root locus method of analyzing a system with such @ non-linearity 
consists of construction of the linear root locus and determining the 
point where the root locus crosses the JW axis. 


The system characteristic equation is 


5° * 35° +2st+tlON=0O0 , (l-A=22)) 


Substitution of $=JW and simplifying 


=o Ww? Spe +)\2W +10 =O ee (f-A-2.2) 


By requiring that both the real and imaginary part of equation I-A- 


22 go to zero independently, 

ao Ww? + ) 2) =. (1-A-2. 3) 
and 

-3w* lO Nesom. (T-A-2.4) 


Dividing equation I-A-2.3 by J& the frequency is obtained as 
UW) <= 1.414 rad/sec. 
By substituting W) into equation I-A-2.4 the "pseudo - average” 
gain over a complete cycle is obtained as 
Nave = 026 
from which the approximate amplitude of the limit cycle is determined 


as 


en = 5 = fee 1) Volts , 





A sketch of the root locus solution is given in Figure I-A-2.1. 
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Root Locus Solution for Saturation in Error Channel ~~ 


Figure T[-A-2.1 


Loop Transfer Function with No Zeros 
The root locus method also gives a stable limit.cycle by the follow- 
ing reasoning. If the input signal is large enough, N will be small 
enough such that the instantaneous root location is at rf, in Figure 
T-A-2.1. The system is stable and the oscillations decrease. As the 
oscillations decrease N increases moving the instantaneous root into the 
right half plane. The system is now unstable, the oscillations increase 


causing N to decrease, and the instantaneous root moves back into left 


° ' 


10 





half plane. Eventually a dynamic equilibrium is established in which the 
"average" root location is at the imaginary axis cross-over. If the 
input signal is small, N is large placing the instantaneous root in the 
right half plane at fg. By a similar reasoning the “average” root 


location over a cycle will again lie at the imaginary axis cross-over. 


11 





I-A-3. Describing Function Method. 

Application of jescxibane function theory to the given system entails 
plotting the linear loop transfer function and the negative reciprocal of 
the describing function on the gain-phase plane (Nichols plot). Limit 
cycles of the system are then given by the intersections of the two curves 
on the gain-phase plane. 

The describing function for a saturating element with unity gain in 


the linear region is 
Gp = =[Sin R + RV |= Rz | Wen (T-A-3.1) 


where 


Esat 


R= — (1-A-3.25 
E max e 


The linear loop transfer function in frequency response form is 


G Ww) = Sapa Re CAPS (1-A-3.3) 


Since G, has no phase shift thenl/G. has an associated phase angle 


D 
of -180°, and the 1/G, curve lies entirely on the -180° axis of the gain- 
phase plane. Consequently, the only region of the linear loop transfer 
function which is of interest is the region in which the phase shift is 
-180°. The values of gain in db and phase of the linear loop transfer 
function in the region of -180° which were used to plot on the gain- 


phase plane are given in Table I-A-3.1. These values were obtained from 


a Bode diagram plot of the linear loop transfer function, 


12 





TABLE I-A-3.1 


| 1Gowdlae | ou sos _| 





The values of 1/G, in db were calculated for various values of R 
according to equation [-A-3.1 and appear in Table I-A-3.2. ‘The gain- 


phase plane of the system is shown in Figure I-A-3.1. 


TABLE I-A-3.2 
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By interpolation of Figure I-A-3.1, the intersecticn occurs at 


Ud 


1.41 rad/sec 


| 


and 


R = 0.495 o 


Thus from Equation I-A-3.2 the amplitude of the limit cycle is 


erie 10,1 “volts. 


Assuming the 1/G, curve to be the locus of the critical point of 
the system, the describing function method predicts a stable limit cycle 
by the following reasoning. If the input signal, O21 is large RK will 
be small and the instantaneous critical point lies above the G (yw) curve 
say at point R, in I-A-3.1. The system is momentarily stable due to a 
positive phase margin and the oscillations will decrease, As the oscil- 
lations decrease, R increases and the critical point moves below the 
G (Jw) curve, The system is now unstable due to a negative phase 
margin and the oscillations increase moving the critical point above the 
G()w) curve. Eventually a dynamic equilibrium is reached in which the 
"average'' critical point over a cycle lies at the intersection of the two 
curves. If the input signal is small, R is large and the critical point 
is momentarily below the G Qu) curve at Re » «the system is unstable 
and the oscillations increase. Similar reasoning shows again that the 
"average'’ critical point over a cycle will eventually lie at the inter- 


section of the two curves. 
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T-A-4, Analog Computer Simulation. 
The simulation of the given nonlinear system is given in Figure 
I-A-4.1. The coefficient pot settings with associated resistances and 


capacitances for real time and real magnitude scaling are given in Table 


I-A-4.1., 


TABLE I-A-4.1 


Pot Setting ASsociated Elements 












For input Step of E. volts and 


Ry = | Meg, Rey + | Meg 





1 Még 





1.0 1 Meg 
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yr 
tk 


h 1 Meg C.. = sie 





EtG, 


x 
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1 Meg C = 


= 





a0 


vs) 
iy 


Fe O.1 Meg Ci, =e 
bald, 


mx 
i{ 





0.5 Meg C 
1.0 





Rg = I Meg C = 1 uf 










a = 0.05 
Initial 
Setting 


For a 5 Volt Saturation 
Voltage 





To accurately set the back-biasing pots, "sa", a gine wave gererator 


was connected tc the input of the saturation simulator. The twa pots 


were then adjusted until the saturated cutput voltage was exactly 5 volts 


on a recording device. The frequency of the generator was set 


and the static characteristic of the saturation simulation wae obtaines 


by an X-Y recorder. This characteristic appears as Figure {f-A-4.2., 


The system was first excited with a 2 volt step input 
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system was found to go in to the Limit cycle shown in recorder trace 
I-A-4.1. The system was then excited with a 20 volt step input and wae 
found to go into the same limit cycle (recorder trace [-A-4,.2). From 
the recorder traces the limit cycle is a stable limit cycle with an 
angular frequency of 1.395 rad/sec and 10.0 volts amplitude. The phase 


portrait for the system appears in Figure I-A-4.3. 


18 


ee 








. 


2009) nz NIS Of, 


I 


on hy 


FAYND DIASIYTLWEVHD 
GP eevs NOL en ts 
Bata aye | 


| 
| 








| 


; 
ie 
a4 
| 
| 


me nn oe tf eee 
' 











~ 


§ 


METS TAS Brie cis UIT J MEI /Ssaton ce a sit TOA a an : 
Q°UTSTART/ PS S 7 Wy T PUET/ISSTON $°O xX MUTT JIICA T - 3S SUTT/ITON - SF -SaTVOS 
@ 


? = : » 3 ij 72° a) > Saal - 7: = 7 + 
sNGUZ da3g sqzo, G $Od0Z ON YIIN UOT4OUNY sazsuerzy dooy ‘UOTZEANZRS LauueyD szcazsyq 
5 s =a) 4aw Ast 3 


L°v-Vel B87, Aapzovay 


ea ee a Se SE EEE EERE 
sey e QeEEY Gerster et Eesneebey cess ser eeucnteuacuaeuns 
Ta iatar ea gesuiarseasga Tay at geese ater eral pian a 
TEER A Corer A rN . 
| ; | . 
FEEL LLL Risers wen ned PEELE ALEPH 
Jee ees See ee CECE OEE rSEEn 
SLES ST EGE REESE! CEGERERE LeGeenenEncneee noo aececce 
ECR ACEC APE CCAR ER | 
FENSS25SS ENG eee ee eee Gee Ne eNee anne 
EEE REE ere 
EERE EERE 
EE EELEL ELI EEL EL Le reese easel cE 
eS ees EE) EEE ee EELL ULL LEE 
aGS=", SESS 36 =e). @SSo 2 See eee. SSS ese Sor eSeSneneeneee 








wee at be oe 


JTC CNEL ULLAL 
POA OEE EECA EEE ERNE 3 
COZ NGG REECE VENEERS EPL 
SEE ssp EEE EEE 
EEA EEL EPP PEPE 











a 
ry 
ITssen 1 - Z SUIT /ILOS Si SATYOS 
s 


STA YT . oC Se : 
{ an s ¢: ji ads wos -) Ps 4 ‘ & 
935 JIPA OZ “SO4a7 ON YQIM UCTioung asjsuezy, dooT ‘uotzeanges 


EPEC EEE EE EE EEEEEPEECECECEr Lee erEpeEe 
Ge -£5 2] 50086] = eS6S====5= 48 se Se eRe Seas eeeae 

BEC E RPE Cee pee 

ae6 SHHAIEEEHE 


- 


TLLLNLLUCLLLILLLIELLLULAL 
PEEEPE REE re 
SHEE 
BESsolsieeereeinien 


SFR ES SRR ee SSeS eee 
“HAHEI CE EEE EE = 
Ty Saeco 5255nmm 


BEBEs 52 =5s58abseoeSee es Senene' 

BD SDS SSeS Sag Seo SSS 22 oe Se ee ee ee eo ee 
Be ee a eee eT | Ty 
EAL Ie EIN EER EERE LEAL ACL ALL TIAL 

5] > Sa ER aataSee Me IGS SERS, Sar RGSS NER RGREee 
Selee eee aes St \= Sha) Seem |e eres rr Frei 

SRA EEE SEERA es ae a AE 5 

Jee eVa81 4225 7255\4e087e5 a eea eee eee eee 

FIP EEL LEE EEL grave Pepper | 
fee eae ee eee ee El a 

i = SSeS) SBSE=]= SSS Se2eSeaneeh—.. BBS ease Saeaeerecaek 


RE CCAEE AEE EACE Enc bbe HEA na a5 SRE CEBSEe 
Se ae0n0 am i 








FONG AACA EA ECLA 
Ae ENA SSN ERE 
MTA AAT A eet FP Pd PE ed PE 








ef ee ea ae =) fe 
2S Se eS eae 





Soe S4eH se E Easel 
=== =25G BESS =e ae eaeneae 
EEEEE ET PEEEEE DEES EEE 


Recotder Trace T[-A-4, 3 


Limit Cycle Measurement, Error Channel Saturation, Leep 
Transfer Func eo with No Zeros, - 


Scales: & - 1 Volt/line X =- 0.5 Volts/ Line 


y 


AN 2 i h ¥ ‘ rm a oY éy _ 
Llime <= oe Seca nc/Di IVLSiOon 





i ee 











a oe ee @; 
|  's0NIZ 
ON HLEM NOLLINAA Lt 
YAISNVYEL dooT "TANNYHD | 
vows NI NOL YOLYS | 

= LIVSLYOd | ei i Hd 





} 
i _ j ae 
Le (| 
Abi , 
a 





Ad re) nl (br 
HA O'OL =F 


Sfp 19 LET EM 
aL ELL PRT OAD LI 


— anevis 





| —_ mal : t | { 7 ' | 

| Sy Ue 

: ores deat 
rar ts cp ud ag a a 


re a I ree me 
} : : 








I-A-5. Comparison of the Four Analytic Methods in Limit Cycle Pre- 


diction for the System of Figure I-A. 


Method Frequency Amplitude 
(rad/sec) (volts) 


Mitrovic's 
Root Locus 
Describing Function 


Analog Simulation 
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I-B. Loop Transfer Function with One Zero. 


The block diagram of the system chosen to analyze is shown in 
6 


Figure I-B. 


% 
*¢" 


Go. + E 





_Figure I-B 
Saturation in Error Channel Loop 


Transfer Function with One Zero 
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I-B-l1. Mitrovic's Method. 
If N is the instantaneous gain of the nonlinearity the character- 


istic equation for the system is 


57+ 1357+ 445°+ (324 300N)5 + GOON=O, (1-B-1. 


Assuming the last two coefficients to be variable, the characteris-~- 


tic equation becomes, 
S44 19.56°444 62 +,s+ BL =O 28 >1'. 


where 


B, 


32 + 300 N (i-5-e 


il 


and 


Bo 


From Appendix A, Mitrovic's equations for BY and B 


th 


600 Ng (I-B-1. 


1 2tes 


Bo=-[44 uw, F (4) + [3w,°D, (4) i Lin4 Py (8) ] (1-B-1. 


and 


6, = 44 wy Q,($) + Swe D,(8) + Ww? D, (4) , (I-B-1. 


Substitution of the values of the @ functions from Appendix A for 
%g = 0 gives, 
Be = 44 Wr - Wy (I-B-1. 
and 
6, = 13 Wy, z - (I-B-1 
Thus the sets of parametric equations I-B-1.3, I-B-1.4, I-B-1.7, 


and I-B-1.8 describe the stability curve and the working point locus on 
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2) 


3) 


2, 


7) 


8) 





the B versus B, plane. To find the intersection of two curves solve 
O 


1 
Equation I-B-1.8 for Coe and substitute in Equation I-B-1.7 to obtain 


a 
Be = ah B, "eon (1-B-1.5) 


Solving Equation I-B-1.4 for N and substituting in Equation I[-B-1.3 


gives, 


EY = "S2anO;S ib, (1-B-1. 10) 
At the intersection, equations I-B-1.9 and I~B-1.10 must be satis- 
fied simultaneously. Solving Equation I-B-1. 10 for By and substitut~ 
ing in Equation I-B-1.9 yields, 


B,* — 2348, - 10816 = 


gives 


© 


(I-B-1.11) 
Solving for By 


By = 27 ons eand By == 39.56 
The negative value is meaningless, thus substitution of the posi- 
tive value into equations I-B-1.8 and I-B-1.3 gives the frequency and the 


"average" gain respectively, as 


Wy - 4,587 rad/sec 


and 
Nave = 0.8051 
The approximate amplitude of the limit cycle is given by 
e = Psat 
max N 
or 
a0 
anne = 9.8051 ~ 6.2) Volts 

The graphical representation of the solution is shown in Figure 

I-B-l.l. 
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Loop Transfer Function with One Zero 
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IT-B-2. Root Locus Method. 


The system characteristic Equation is 


4 


Substitution of s = )W and simplifying, yields 


Ud" -S13 wy? -44y) 2 +3(324300N)W+ BOONE O (1-B-2.2) 


By requiring that both the real and imaginary parts of equation 
I-B-2.1 go to zero simultaneously, then 


Ww 44 W*+ 600 N 


and ~ dwt 32W+ J300WN 


Dividing Equation I-B-2.4 by JW and solving for WN gives 


no Ge 
300 300 e 


Substituting this into Equation I-B-2.3 gives the following equa- 


N = 


tion for WU 


z 
uy -18W- 64 =0 - (1-B-2.6) 


Solution of Equation I-B-2.6 for w* gives 


To = NEC, 
from which 
QD = 4.587 rad/sec. 
Substitution of wo into Equation I-B-2.5 gives the "average" non- 
linear gain over a cycle as 
Ns 0.8051 


From this the approximate amplitude is obtained as 


5 | 
& “ase. ) ORSOB = 6,21 Volts 


A pictorial representation of the root locus solution appears in 


Figure [-B-2.1. 
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s + near + Wee + (32 + 300 N}s + 600 N = OG. (I-B-2.1) 


0 (1-3-2. 3) 


0. (I-B-2.4) 


(1-B-2.5) 





The root locus method also predicts a stable limit cycle by follow- 


ing similar reasoning to that given in Section I-A-1 of this Chapter. 
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Root Locus Solution for Saturation in Error Channel 


Loop Transfer Function with One Zero 
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I-B-3. Describing Function Methed. 
The describing function is the same 4s that given in Equation I-A-3.1 
of this Chapter. 


The linear loop transfer function in frequency response form is, 


G (Jw) = 18.75 (0.5 Jw +1) pean 


IWOW +O IW+N(OI2SIw+!) —« 


Table I-B-3.1 gives the magnitude and phase of Equation [-B-3.1 
which were used to plot the linear loop transfer function on the gain- 
phase plane of Figure I-B-3.1. These values were obtained from a Bode 


diagram plot of the linear loop transfer function. 


TABLE Sioa 


[600 aes 








Table I-B-3.2 gives the values of 1/6, in the area of interest. 


These values were obtained from Equation I-A-3.1. 


TABLE I-B-3.2 
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By interpolation of Figure I-B-3.1, the intersection and thus the 


limit cycle occurs at 
dy = 4.55 rad/sec. 
and 
R = 0,682 


From Equation I-A-3.2, the amplitude of the limit cycle is then 


5) 
a Tau SONGSD = 7.33 Volts 


The describing function predicts a stable limit cycle by following 


similar reasoning to that given in Section I-A-3 of this Chapter. 
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I-B-4. Analog Computer Simulation. 


The analog simulation of the given nonlinear system is shcewn in 


Figure I-B-4.1. The coefficient pot settings with associated resistances 


and capacitances for real time and real magnitude scaling are given in 


Table T-B-4, 1, 


Pot seti e 
ae 
1 100 
ay = eo 
a. = 1518, 
ae Ale 
ag = 0.4 
ag = 8, 
ag = 0 
ag = 0.3 
aga 1 
210 7 L518, 
2447 0.8 
419 = 1.0 
ai3 7 L.0 
a Ges. 


initial Setting 





TABLE I-B-4.1 


Volts 

17 «= 1 Meg 
: 1 Meg, 
1 Meg, 


= 0.5 Meg, 


O.1 Meg, 
1 Meg, 
O.1 Meg, C 
O.1 Meg, 
= 1 Meg, 


1 Meg, 


Rig 


For a 5 volt 
Saturation Voltage 


connected to the input of the saturation simulator, The two pots were 


then adjusted until the saturated output voltage was 5 volts on 4 re- 


cording device. 
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The static characteristic for the saturation simulator appears as 
Figure I-A-4,2 in seceon T-A-4 of this Chapter. 

The system was first excited with a 4 volt input signal. The 
system was found to go into the limit cycle shown in recorder trace 
I-B-4,1. The system was then excited with a 10 volt step input signal. 
The system was found to go into the same limit cycle (recorder trace 
I-B-4.2). From these two recorder traces, the limit cycle is seen to be 
stable with an angular frequency of 4.58 rad/sec and 7.2 volts amplitude. 


The phase portrait for the system appears in Figure I-B-4.2. 
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I-B-5. Comparison of the Four Analytic Methods In Linit Cycle Prediction 


for the System Given in Figure [-5. 


Method Frequency Amplitude 
rad/sec 


Mitrovic 
Root Locus 


Describing Function 





Analog Simulation 





40 





CHAPTER II 
DEAD ZONE IN ERROR CHANNEL 
II-A Loop Transfer Function with No Zeros. 
The block diagram of the arbitrarily chosen eFarpi wien 1s to be 


analyzed is given in’Figure II-A. 





Figure [T-A 


Dead Zone in Error Channel Loop 
Transfer Function with No Zeros 
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II-A-1. Mitrovic's Method. 


The characteristic Equation of the given system is, 


s*+lls°+26s*4 1654+ GON=O, (1t-A-1.1) 


Assuming the last two coefficients variable the characteristic equa- 


tion becomes, 


2. 
statis? + 2657 +8.54Bo=20 (1T-A-1.2) 


where the parametric equations of the working point locus are, 





Bo = 60 N (IT-A-1.3) 
and 
B, = 16 P (Ti-A-1.4) 
From Appendix A, Mitrovic's equations for BS and By are, 
2 
Bo= -[26w, P (8) + up, Do(¥B) + wn! (4) J] (II-A-1.5) 
and 
as 2 S ake nen 
8, = 26W, (4) + I] Wn h.(4) tT Why BD, (4) 6 (1T-A-1.6) 
Substitution of the values of the @ functions from Appendix A for 
‘A = Q gives the parametric equations of the stability curve as 
Bas Ze =} 4 (II-A-1.7) 
Oa Un Nn oti 
and 
B= Hw,* : (1I-A-1.8) 


At the intersection of the M-point locus and the stability curve, 
equations I[I-A-1.3 and II-A-1.4 and Equations I[I-A-1.7 and IY-A-1.8 
must be satisfied simultaneously. Thus at the intersection, 


2 it ec 
}| Wy =E6 (IT-A-1.9) 


42 





from which the angular frequency of the limit cycle is determined to be, 
Wy, = 1.206 rad/sec. | 
Also at the intersection 
otoe Wyn =eeON | (II-A-1.10) 
: Substitution of , BOA from Equation II-A-1.9 into Equation II-A-1.10 
gives the "average" gain of the nonlinearity over the cycle.as, 
Nowe 0.595 & 
To obtain the approximate amplitude of the limit cycle, reference is 


made to Figure II-A-1.1 which is a sketch of the nonlinearity static 


characteristic. 





Figure II-A-1.1 


Limit Cycle Amplitude Calculation from the 


Dead Zone Characteristic Knowing Nixays 
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= T-A-1,1 
ee 5 +c (1T-A-1.11) 
¢ = 4 cos 45° = 07.07a (II-A-1,.12) 
5 Sin A zs 
a =a (II-A-1.13) 
Ace Tan (1I-A-1.14) 
ave 
and 
B = 180 - A - (180 - 45) , (Ii-A-1.15) 
Entering equation II-A-1.14 with the "sverage'’ gain over a cycle, 
A = Tan) 0.595 = 30.8° , 
Substituting this in Equation II-A-1.15 gives 
B a 457 = 30.8° — Ve 2 e 
Substituting for A and B in Equation If-A-1.13 gives, 
sin 30.8 _ | 
—? Sie ee 
Substituting for "a" in Equation If-A-1.12 gives 
ec = 04707 (107469) = 724024 
And finally Equation II-A-1.11 gives 
e = 5+ 7,402 = 12.402 Volts 
max 
The graphical solution of the intersection is sketched in Figure 
II-A-1.2. 


reasoning. If the input signal is smail, N will be small and the in- 


From the figure the following equations are obtained, 


Mitrovic's method predicts an unstable limit cycle by the following 





Stantaneous M-point will be in the stable region such as M, in Figure 


II-A-1.2. The system is stable and the oscillations will decrease in ampli- 


tude, 


M-point even further into the stable region. Consequently the M-point 


As the oscillations decrease, N becomes even smaller driving the 
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le | 
ae “STABILITY CURVE ($0) 





Ole 142 286— «iB 
Figure II-A-1.2 
Graphical Solution for Dead Zone in Error cnaenet 

Loop Transfer Function with No Zeros 
can ase into the unstable region and all oscillations will eventual- 
ly die out. If the input signal is i N is relatively ares and the 
instantaneous M-point lies in the unstable region such as Mo: The system 
is unstable and the oscillations increase, N increases toward N = 1.0, 
and the M-point moves farther into the unstable region. Consequently, 


the oscillations increase indefinitely in magnitude. 
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II-A-2. Root Locus Method. 
Substitution of s =2W into the characteristic equation given in 


Equation LI-A-1.1 and evaluating the higher powers of J gives, 
4 > 2 , 
WwW -JIlu~-26W 4JI6W+ CON=O, ORR) 


By requiring that both the real and imaginary parts of Equation Ii-A- 
2.1 go to zero independently, 
4 2 
Ww’ - 26W + 60 N = 0 (TI-A-2.2) 
and 
= {fi 
+Jiéw-o . (11-A-2Z. 3) 
Dividing Equation ITI-A-2.3 by JW , gives 
2 
tee, - 16 = 0 (II-A-2,4) 
from which the angular frequency of the limit cycle is obtained as 
() = 1.206 rad/sec. 
Substitution of WO from Equation If-A-2.4 into Equation T{-A-2.2 
and solving for N yields 
ave 
N = 0.595 e 
ave 
To determine the approximate amplitude of the limit cycle, Equations 
IT-A-1.11 through II-A-1.15 are again used. Thus, 
E a as 12.402 Volts 
A graphical presentation of the root locus solution is given in 
Figure II-A-2.1., 
The root locus method predicts an unstable limit cycle by the follow- 


ing reasoning. If the input signal is small, N is small and the loop gain 





is small placing the instantaneous root at Y, in Figure If-A-2.1. The 


system is stable and the oscillations decrease in magnitude. As the 
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x mV, 


\ LINEAR SYSTEM CLOSED 7 
\ LOOP ROOT ; ! 
_ MS Via ra 
? Va 
< LIMIT CYCLE 
J. 4 W= 1,206 her 
~ 7 Nave 3 0.595 


aan +o 


~)wW) 


Figure IT-A-2.1 / 7 
Root Locus Solution for Dead Zone in Error Channel, 
Loop Transfer Function with No Zeros 

oscillations decrease, N also decreases, driving the root location farther 
into the left half plane. Consequently, the system is always stable and 
the oscillations eventually die out. If the input signal is large, N is 
also relatively large making the loop gain high. The instantaneous root 
location is at fq. The system is unstable and the oscillations increase, 
As the oscillations increase, N increases toward N = 1. The root location 


moves farther into the right half plane.. Thus the system is always un- 


stable and the oscillations will increase indefinitely. 
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II-A-3. Describing Function Method. 


The describing function for a dead zone nonlinearity is 


~ = | , = ral Be 0 re 
Gp = =[($ -Sin R-R 1-R2 | jon (1T-A-3.1) 
with 
Pe, oe 
a 2 Cae (ITI-A-3.2) 
where b is the total minus to plus dead zone. 
The linear loop transfer function for the given system in frequency 
response form is 


zo 3,75 | aoe 
G WW) ~ JWGW+1XKO5)W +!YOl25 Jw +!) Sa Cy 


e 





Table LI-A-3.1 gives the magnitude and phase of Equation I[-A-3.3 
for the different frequencies which were used to plot the linear loop 
transfer function on the gain-phase plane of Figure If-A-3.1. These 


values are obtained from a Bode diagram plot of Equation IT-A-3.3. 


TABLE IT-A-3.1 


1G0w) db 
2 ae 


-168 













Gye -174 







4.7 Snes, 





o> -184 





“189.5 


TABLE IT[=A-3.2 





These values were obtained from equation ITI-A-3.1. 
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By interpolation of Figure II-A-3.1, the intersection and thus the 
limit cycle occurs at 


(Ura acletad/ sae, 


and 
R = 0.34 e 
From Equation II-A-3.2, the amplitude of the limit cycle is then 
_ & 
E max 2R 
or 
& = 14.71 Volts 


max 


The describing function method predicts an unstable limit cycle by 
the following reasoning. For small input signals, R is relatively large. 
Assuming the “1/6, curve to be the locus of the system critical point, then 
the instantaneous critical point is at R, in Figure II-A-3.1 and the sys- 
tem is stable due to a positive phase margin. The oscillations will de- 
crease causing R to decrease, thus moving the critical point farther up 
the -180° axis. This motion of the critical point increases the positive 
phase margin. Thus the system is always stable for small signals and all 
oscillations will eventually die out. For large input signals R is small 
placing the instantaneous critical point at Ree The system is unstable due 
to a negative phase margin and the oscillations increase. As the oscilla-- 
tions increase R becomes even smaller moving the critical point down the 
-180° axis and increasing the negative phase margin. Consequently, for 
large signals the system is always unstable and the oscillations increase 


indefinitely. 
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II-A-4 Analog Computer Simulation. 

The analog simulation of the given nonlinear system is shown in 
Figure II-A-1,. The coefficient pot settings with associated resistances 
and capacitances for real time and real magnitude scaling are given in 
Table II-A-4.1. 

TABLE If-A-4,1 


Pot Setting Associated Elements 


For Input Step of E, Voits 





R, = 1 Meg, Rey 1 Meg 
Ry = 1 Meg, Rey = 1 Meg 
R, = 1 Megs Reo = 1] Meg 
Ry = 1 Meg, Coy uit 
Re = 1 Meg, Cry = 1 uf 
Re = 0.1 Meg, Ceo =" [ue 
Ro = 0.5 Meg, Ceo = lowe 
Re = 0,1 Meg, C3 = Laue 
Ry = 0.1 Meg, Ce = eon 
Rig = 1 Meg, Ces = 1 uf 
Ray = 1 Meg, Re3 = ] Meg 
For a 5 Volt Dead Zone Each 
Initial Setting Side of Zero 


To accurately set the back-bias pots "a", a sine wave generator was 
connected to the input of the dead zone simulator. An X-¥ recorder was 
then calibrated for 5 volts per inch on éach axis. The input and output 
of the simulator was then connected to the X-Y recorder and pots “a" were 
adjusted until both break voltages of the static characteristic curve were 


> volts. The static characteristic curve appears as Figure II-A-4.2. 
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The system was first excited with a 15 volt step input signal and 
the system was found to be stable as shown in recorder trace II-A-4.1. 
The system was then excited with a 17 volt step input and was found to be 
unstable as shown in recorder trace [I-A-4.2. The unstable limit cycle 
was found to occur for an input step of 17.07 and is shown in recorder 
trace I[I-A-4.3. From recorder trace II-A-4.3 the limit cycle is seen to 
have an angular frequency of 1.207 rad/sec and an amplitude of 15 volts. 


The phase protrait of the system is shown in Figure II-A-4.3. 
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II-A-5 Comparison of the four analytic methods in limit cycle prediction 


for the system of Figure II-A. 


Method Frequency Amplitude 
| (rad/sec) (volts) 


Mitrovic 
Root Locus 
Describing Function 


Analog Simulation 





a9 


II-B Loop Transfer Function with One Zero. 
The block diagram of the system chosen to analyze is shown in Fig. 


II-B. 





Figure [II-B 


Dead Zone in Error Channel Loop 
Transfer Function with One Zero 
II-B-1l. Mitrovic's Method. ' 
If N is the instantaneous gain of the nonlinearity the character- 


istic equation for the system is 
374 11.554t B15 $° $298 +(B45N)s + 20N =O. (I1-B-1.1) 


Assuming the last two coefficients to be variable, the characteris- 


tic equation becomes, 


5° + Wear ss peg s* +B, s +Bo=O (II-B-1.2) 
where | 
B,= 8+5N : (II-B-1. 3) 
and 
Gee Ze. | mt (LI-B-1.4) 
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From Appendix A, Mitrovic’s equations fer B. and b, are, 


By = — [24 09 7G,(#) +3) Swe Gal) +11, 5100 (@) + ung) ] — (EE-B-2.5) 


and 


B= 24 ur $1(4)+3).5in2Gyl9) + HS DylY) tun Ge (B) . aE-B-2.6) 

Substitution of the values of the @ functions from Appendix A for 
=O gives, 

Sas Ze Ww,” eae (LI-B-1.7;, 
and 
Ce ign Sen (LI-B-1.8) 

The sets of parametric equations [I-B-1.3, 1f-b-1.4, Yi-B-1.7, and 
II-B-1.8 describe the M-point locus and the stability curve on the Ba 
versus By plane. Direct solution of these equations for the intersection 
becomes fairly complicated; consequently, a straight graphical solution is 
performed, 

Table LI-B-1.1 gives the values of B and By for values of Lon 
plotted in Figure II-B-1.1. Elimination of N from Equations [[f-5-1.3 and 
II-B-1.4 yields 

B,= St Ge (Li-B-1.3) 
which is a straight line with a By intercept of + 8 and a slope cf + 4, 
The intersection of the two curves is seen to be approximately at G),. = 0.36 
rad/sec. 

To accurately determine the intersection more values of stability curve 


are calculated and the value of B, thus determined is substituted in Eque- 


j 
tion II-B-1.9. 
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TABLE 1i-B-t2 





The intersection occurs when the two values of By are equal. These addi- 


tional calculations appear in Table II-B-1.2 


TABLE II-~B-1.2 


8.72017 5.74897 


025607 110502355 8.13459 8.09419 


Ge566 (OR O5657 oalta)S) 116 8.23428 


10.12679 8.20816 8.50716 





By interpolation the intersection occurs at 
Wy = 0.5673 rad/sec 
and 
B_ = 8.140 
O 
From Equation II-B-1.4, the average nonlinear gain over a cycle is 


N = 0.407 
ave 
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Since the nonlinearity is the same as the one of Section A-1 of 
this chapter, Equations [I-A-1,11 through II-A-1.15 are used to obtain 
the approximate amplitude of the limit cycle. Substituting Noe into 


Equation II-A-1.14 yields 


A= Tan) mONcOTe encom 
Substituting A into Equation II-A-1.15 yields 
B= 92 ene oe 


Substituting A and B in Equation II-A-1.13 gives 
Sin 22.15 
a =) 2 Sain Se cc 
Substituting "a'' into Equation II-A-1.12 gives 


¢ = 0.707 (4.65) =53542 


Finally Equation II-A-1.11 gives 


e 


Applying the same reasoning as given in Section II-A-1 to Figure 


= 5 + 3.42 = 8.42 volts, 


ma 


II-B-1.2, Mitrovic's method predicts an unstable limit cycle. 
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I¥-B-2. Root Locus Method. 
Substitution of Sajw into Squation TI-E-1,1 and simplifying the 


higher order J terms gives 


yw HHS e133). & wo — 29 + G+ SN) WwW F2ON=O, (1y-B-2.1) 


By requiring that both the reali and imaginary pa: 


Gs 
Fy 
— 
La 
() 
+ 
§ 3 

2 
e 
{5 
ag 
bate 
'@) 
toh 


[f¥-B-2.1 go to zero independently gives 
5 3 Pee 
JOP) 31,5@"+)(84 SNW = o (17-3-2,2) 


and 
“ 2 ar. a a 
lS ww -29W + 20N - 0 (iI 4-28.35 


Dividing Equation II-B-2.2 by JW and solving for 5n giwes 


Substitution of Equation I[f-B-2.4 into Equation [EY-B-2.3 and sinplify- 


ing gives 


4 2 A 
vee, 3D) + 97 UV a = 0 PY 
Z 
Solving this equation for KJ) gives 
2 neu 
UD = 0.3219 
from which 
UY) = 0.567 rad/sec 
Substituting the value of W back into Equation I[i-B-2.4 gives the 


"average" nonlinear gain over a cycle 4s, 
N = 0.4072 
ave 


To find the approximate amplitude of the limit cycle, Equations 


TI-A-1.11 through II-A-1.15 are again used. These equations yield 


(¥) 
if 


&.41 volts 
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A ‘graphical presentation of the solution is sketched in Figure II-B- 
Zel. 


Applying the same reasoning as given in Section II-A-2 to Figure II- 


B-2.1, the root locus method predicts an unstable limit cycle. 
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II-B-3. Describing Function Method. 
The describing function for dead zone is given in Equation [{f-A-3.1. 
The linear loop transfer function for the given system in frequency 
response form is 


’ 25 JW +] 


JU) (20+ IYIWF!)(O.5IW+! ONZE IW +1) aaa 


Table II-B-3.1 gives the magnitude and phase of Equation Tf-B8-3.1 
for the different frequenctes which were used to plot the linear loop 
transfer function on the gain-phase plane of Figure {{-B-3.1. These 


values are obtained from the Bode diagram plot of Equation IY-P-3.1. 


TAELE {I-B-3.1 
\GQv)! db 








Table II-B-3.2 gives the valves of 1/6, plotted in Figure I[T-B-3.1. 


These values are obtained from Equation [{[-A-3.1 for warious valwes of R. 


0.4027] 
0.3910 





By interpolation of Figure II-B-3.1 the intersection and thus the 
limit cycle occurs at 
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0.56 rad/sec 


e 


and 


R O74S5 “s 


From Equation II-A-3.2, the amplitude of the limit cycle is then 


10 
a = FO -a5sj 7 > (10-1 volte. 


Applying the same reasoning given in Section LI-A-3 to Figure II-B- 


3.1, the describing function predicts an unstable limit cycle. 
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LI-B-4 Analog Computer Simulaticn. 

The analog simulation of the given nonlinear system 1s shown in 
Figure If-A-4.1. The coefficient pot settings with associated resist- 
ances and capacitances for real time and real magnitude scaling are given 


in Table IY-A-4.1. 


TABLE [f-A-4.1 


Pot Setting Assoctated Elements 


For Input Step of z. volts 


= 1 Meg, = 1 Meg 


Rey 


1.0 l Meg, ; 1 Meg 


Rey 
moO 1 Meg, 
0.4 
0.8 
Ln0 
lO 
0.5 , = O.1 Meg, C 
1.0 
1.0 
Exo 
0 
= U8 
= 1.0 


au—7 0.05 
Initial setting 





The back-bias pots "a" were accurately set in the same manner as 
described in Section If-A-4. The static characteristic of the non- 


linearity is shown in Figure II-A-4.1. 
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The system was first excited with a 10.2 volt step input and was 
found to be stable as shown in recorder tramce [T-b-4.1. The system 
was then excited with a 10.6 volt step input and found to be unstable, 
as shown in recorder tréce II-B-4.2. The unstable limizx cycle was 
difficult to obtain but was found to lie between step inputs of 10.3 volts 
and 10.4 volts. These two responses are shown in recorder traces I[I-~F-4.3 
and II-B-4.4. The average of the two traces show the angular frequency to 
be 0.552 rad/sec and the amplitude, using the first two peaks of each 


trace, to be 9.9 volts. 


The phase portrait of the system is shown in Figure TI-B-4,2. 
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Recorder Trace II-B-4.1 
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Loop Transfer Function with 
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Input. ' 
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Time - 1 Second/division 


r Function with One Zero 


Ti-B-4,3 


X - 0.5 Volt/line 
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Recorder Trace 


LIMIT CYCLE, Dead Zone in Error Channel, Loop Tra 
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II-B-5 Comparison of the Four Analytic Methods in Predicting Limit 
Cycles for the System Shown in Figure II-B. 
Method Angular Frequency Amplitude 
rad/sec Volts 
Mitrovic's 0.5673 
Root Locus 0.567 


Describing Function 0.560 


Analog Simulation 0.552 








CHAPTER III 
IDEAL RELAY IN ERROR CHANNEL 
III-A Loop Transfer Function With No Zeros 
The system arbitrarily chosen to analyze is ax in the block 


diagram of Figure III+A. 





Figure III-A 


Ideal Relay in Error Channel 


Loop Transfer Function With No Zeros 


III-A-1 Mitrovic's Method 
Assigning N as the instantaneous gain of. the relay, the character- 


istic equation can be written as ‘' 


or 4 
$+9.5S Page Revi 1S | che =O 7 (III-A-1.1) 


Assuming the last two coefficients variable, the characteristic: 


equation becomes, 


- 3 
San eG si t 30.5 5 oise = GB, Ss +Bo =O; (III-A-1.2) 
where 


B = 12N - (III-A-1.3) 
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and 


B = eZ > ( 


is 


LIY-A-1.4) 


¢ 


From Appendix A, Mitrovic's equation for B and B, ars, 
Bo = -[37un (4) +30.5t1°,(4)4 7.5un Pal¥) + Lon Bal) | (111-A-1.8) 

and 
B= 3%, 9,4) +20.5tin bs (8)+9, 5 bn Pl 4) +wn' 66) (EY TH-A-1.83 


Substitution of the values of the @ functions from Appendix A for 


¥ = 0 gives the parametric equations cf the stability curve as, 


3 2 \ = & ala fh = 3 
B . a 37 ZOM = ye Ln (LIT-A-1. vy 
and 
: 2 d —_ - 
B, = 30-5 Wn - Wy, (L1f-A=1.8) 


To find the intersection between the M-point locus (Equations ILI-A- 
1.3 and ITI-A-1.4) and the stability curve (Equations I[fI-A-1.7 and LII- 


A-1.8), the graphical approach is chosen. Table {YYf-A-l.1 vives the 


for various that are plotted in FPigewre TYf-4-1.1. 
n 


values of Band B 
fe) 1 


The M-point locus is seen to be a straight line parallel to the B_ axis 


at By = 12. From Figure IIf-A-1.1, the intersection is geen to occur at 


about Q), = 9.62. 


~=J 
00 


TABLE ITI-A-1.1 





To accurately determine the intersection more values of E and 
B, are calculated for Won close to 0.62. These values eppear in Table 
IIf-A-1.2. By interpolation from Table II{-A-1.2 the intersection and 
thus the limit cycle occurs at 
Gy = 0.6315 rad/sec 


and 


B = 13,24 
O 


TABLE IIT-A-1.2 


12.8191 


PJ aloo) 


132259 11.9854 


13.2631 12.0229 





The "average" nonlinear gain over a cycle is determined from 
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Equation III-A-1.3 as, 
N = 1,103 
ave 
From this the approximate amplitude of the limit cycle is 
eo = 4,53 volts 
max 

Note that the M-point of the linear system in Figure IYI-A-1.1 shows 
that the linear system is stable. It is well known that an ideal relay 
gives a limit cycle for a linear system of order three or abowe and the 
occurance of a limit cycle is independent of linear system stability. 


This is'easily explainable on the Be versus B, plane. The ideal relay 


1 
has an infinite maximum instantaneous gain whereas the nonlinearity used 
in Chapters I and II have a maximum instantaneous gain of one. Thus on 


the BY versus B, plane the M-point locus is unlimited for the ideal re- 


1 
lay but does have a maximum valve for the other noniinearities. Conse 
quently, an infinite straight line segment must necessarily intersect the 
stability curve at some point producing a limit cycle. 

Mitrovic's method predicts a stable limit by the same reasoning used 
in Chapter I, Section I-A-l. This same argument can be appiied to the 


ideal relay since this nonlinearity can be envisioned as saturation with 


an infinite gain in the linear region. 
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TtT-A-2 Root Locus Method. 
Substituting § =JW into the characteristic Equation ili-a- 


evaluating the higher order powers of J yields, 


Pe eon esr 4 Si2gwW+ti2N=o,; 


Requiring that the real and imaginary parts of Equation ITf- 


go to zero independently gives, 


su9* -$30.5 uf +)12wW=0 


and 


ip Sy) eee: hie = O 


¢ v bas) ¢ = Z 
Dividing Equation III-A-2.2 by )w and solving torW give 


by” = 0.3986 


¥ 


from which 


by) = 0.6314 rad/sec 


e =e 4,533 voles 
max 
The graphical presentation of the ssiution ty root Tscus 


is shown in Figure III-A-2,1. 


the rooce locus method predicts a stable linit cycle 
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Same arguments to Figure TITI-A-2.1 that ave given in Chapter 
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LII-A-3 Describing Function Method. 


The describing function for an idsal relay is 


Gp 2) ae / ‘ey (INI-A-3.1) 
We Spray f 


The linear loop transfer funccion for the given system in frequency 


response form is, 


_ 
Me) = Naa), FIL eae 
G (0) DUA JW+1 JO. 5 Jw-r 10.333 W+FINKOAWSIW+!) 5 prey ag. 25 


Table IIT-A-3.1 gives che magnituce and phase of Equation ILT-A-3.2 
for the different frequencies which were used to plot the linezr lesp 
transfer function on the gain-phase plane of Figure III-A-3.1. These 


values are obtained from a Bode diagram plot of Equation TYI-A-3.2. 


TABLE Lii-A-3el 





Table III-A-3.2 gives the values of 1/G in the area of interest as 


rh 
ab 
F 
ry 
© 
Fr 


Obtained from Equation [II-A-3.1 for various values o 


TABLE TIL-A-3.2 





By interpolation of Figure III-A-3.1, the intersection and thus 


the limit cycle occurs at 


iD = 0.635 rad/sec 


Ss = "5, c8.Olts ~% 
max 
The describing function method predicts a stabie limit cycle by 


applying the same arguments in Section [I-A-3 of Chapter { to Figure [it~ 
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ITI-A-4 Analog Computer Simulation. 

The analog simulation of the given nonlinear systex is shewr in 
Figure IZI-A-4.1. The coefficient pot settings with associated resist- 
ances and capacitances for real time and real magnitude scaling ara given 


in Table TYI-A-4.1., 


TABLE IYI-A-4,1 


Pot Setting 


GC 





Associéted Elements 


















For Input Step of E, 


a q = = = i] ne 


(See discussion) 
0.4 
Os 
0.3 
ae 
1.0 
Gino 
1.0 
0.4 


fe 


0.05 (Initial 
setting) 
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The simulation of the ideal relay is seen to be @ simple saturation 
simulation. However, the gain of the linear portion is made very high 
by the insertion of coefficient pot B 0 To set the back-bias pots ' 


and pot a sine wave generator was connected to the input of the re- 


gs 


lay simulator, The output was monitored and the back-bias pots "a" 


WELE 
then adjusted until the output was limited at + 5 volts. The output and 
input was then connected to an X - Y recorder to obtain a static character- 
istic curve. Pot a, was then adjusted until the slope of the Linear por- 
tion was as nearly vertical as possible without the amplifter overloading. 
The static characteristic curve appears as Figure T[IY-A-4.2. 

The system was first excited with a 1 volt input signal and was found 
to build up to the limit cycle as shown in recorder trace TIf-A-4.1. The 
system was then excited with a 10 volt step input and was found to go into 
the same limit cycle as shown in recorder trace IIY-A-4.2. From these two 
recorder traces the limit cycle is seen to be stable with an angular fre- 


quency of 0.628 rad/sec and an amplitude of 5.7 volts. 


The phase portrait of the system is shown in Figure Y{I-A-4.3. 
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IQTI-A-5 Comparison of the Four Meth 


the System of Figure ILI-A. 


Method 


Mitrovic's 


Koot Locus 


Describing Functions 


Analog Simulation 


© 


d 


Ss 


in Predicting Limit Cycles for 


Angular Frequency 
rad/sec 








III-B .Loop Transfer Function with One Zero. 
The system to analyze is shown in the block diagram of Figure 


III-B. 





Figure III-B 


Ideal Relay in Error Channel Loop 


Transfer Function with One Zero 
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TIf-B-1. Mitrovic’s Method. 
Assigning N as the instantaneous gain of the relay, the character-~ 


istic Equation can be written as 


Ee 3) EGNOS +1240 5° ¢ 1934 5” +(1024+8000N)s +4000N=O. 
(ITI-B-1.15 


Assuming the last two coefficients are variable, the characteristic 


equation becomes, 


5°45 319°4310s" $1240 S°+ 19845748, S$ +@o= 0 


(LTY-B-1. 2) 
where 
B =-4000 N (1¥Y-B-1., 3) 
and 
B, = 1024 + 8000 N (TIT-B-1.4) 


1 


From Appendix A, Mitrovic'’s equations for 3, and B, are, 


~~ 


Bo= -[1984 Wag, (4) + 1240wn Do (4) + 316 4d, 8) 


. 31 wd, (4) TORO o d-(4)| (111-B-1.5) 
and 
2 
B, = 1984 Wn Go(4) + 1240 wns) + 310 Wy? Dy CG) 
+31 WAG 8) 4 Wn PeB) EEE te 
Substituting the valves of the functions from Appendix a for 4 — 6 
gives the parametric equations of the stability curve as, 
i. Be , (eS PAE eat PN 
6,= 1984 Un ~31Ow,* tw, e (1TT-B-1.7) 
and 
ox 
B, = 1240 un - 3a," (111-B-1.8) 
To tind the intersection between the M-point locus (Equations IYfy-E-i 
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and Iii-B-1.4) and the stability curve, the graphical appreach is chosen, 
Table IIY-B-1.1 gives the values of BY and by for the stability curve as 


plotted in Figure III-B-1.1, 


TABLE ILI-B-1.1 





To obtain the M-point locus solve Equation ITI-B-1.4 for W and 


Substitute into Equation III-B-1. 3, 


B = 0.5 B, - 512 (IIL-B-1.9) 
Equation [TI-B-1.9 is seen to be a straight line with a Slope of 0.5 and 
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a B. intercept of 1024, 


It is seen from Figure [LI-B-1.1 that the intersecticn occurs 4t 


reughly Wn, = 2,24, To more accurately decermine che inrersection, 
values of By and By are calculated for (Uy close to 2.24, Tre value 
of By is then substituted into Equation {II-B-1.9 and the intersecticn 


occurs when the two B's are equal. These additional values are given 


in Table LII-B-1.2. 


TABLE III-B-1.2 


BA (Stability 8 (M-point 
curve ) Locus ) 





Thus from Table III-B-1.2, the frequency of the limit cycle is 
Oy, = 2.25 rad/sec p¢ 


Substituting BA at Wy, = 2.25 into Equation III-B-1.2 


y 


gain over a cycle is determined as 


_ 2229 


ee ~ 4000 ~ ei I 


from which the approximate amplitude is determined as 


5 
© a. = dee sa ee Cae oe. 
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the “average” 





Mitrovic's method predicts a stable limit cycle ty the same argu- 
ments applied to Figure III-B-1.1 as were given in Chapter I, Section 
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Iff-B-2 Root Locus Method. 
Substituting s = Ju) into the characteristic Pquation Lif-B-1.)3 
and evaluating the higher powers of J) yields, 


-w +) Z1w°+ 2idius ~JI240w> —|9840 > 
+)(1024 + Z000N)W)+4000N =O, (IT2-8-2.1) 


Requiring that the real and imaginary parts of Equation YIT-5-2.1 
go to zero independently, gives 
G “ q 2 i it 
—0 + 310w —!1984W9~+4000N=6 CLLT-B-202, 
and 
Ss s aye ate rrtiny 4 
J3) Ww? -)1246W ~+ 3(1024+8000N)wr=0 . (LLT-5+2, 5) 
Dividing Equation ITI-B-2.3 by Jw and solving this for 4000 K 


gives, 


~ 


oe 2 o 4 e717 wee 
4000N = 626 w°-S12=15, FW (ITI-B-265) 
Combining Equations YfIY-8-2,.3 and TLT-B-2.4 gives 
b 4 on arcuate 
W) - 294.5 Wd +1360 +512 =O. (LYE-B+2, 5) 
2 | 
To solwe for W” from Equation {Yi-E-2.5, a digital computer pro- 


gram was used and gave the roots, 


Z rs 

BY = 5,062 
Z 

bd) ae 0.455 
Zz FA) A iY i ~~ 

WwW =299,55 


; . 2: ae | 
Substituting pO into equation I[[l-#-2.4 gives the “average” vain 
of the relay over a cycle as 


Ve 00,5574 
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ee citrine + Aer te D 


from which the approximate amplitude of the limit cycle is determined to 


be | 


eS = 8.972 volts 
y max 


A graphical presentation of the solution is sketched in Figure 


ITI-B-2.1. The root’ locus method predicts a stable, limit cycle by 


applying to Figure III-B-2.1 the same arguments given in Chapter I, 


Section II-A-2, 


+) 


\ 
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en 
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Figure III-B-2.1 
Root Locus Solution for Ideal Relay in Error Channel, 


Loop Transfer Function with One Zero, 
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IQlI-B-3 Describing Function Method. 


The describing function for the ideal relay is given 


TYY-A-3.1. The linear loop transfer function for the piven system in 


frequency response form is 


G (ww) 3.91 (2Jwtl) _- 
~ JUXIW+1 (OS IW +1)G. 259 W+1)(0.12F Just! YOO6Z25IW+1) 5 


Table IIf{-B-3.1 gives the magnitude and phase of Equation ITY- 


for the different frequencies which were used to olot the linear loop 


transfer function on the gain-phase curve of Figure [T{L-k-s.1 


values are obtained from a Bode diagram plot of Equation TYi-#-3.1. 


TABLE YIiTt-B-3.1 





Table I{I-B-3.2 gives the values of 1/C_ in the ¢rea of inter 


D 


obtained from Equation II1-A-3.1 for various valves of err 
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By interpolation of Figure III-B-3.1, the intersection and thus 
the limit cycle occurs at 
) = 2.24 rad/sec. 
and 
E = 11.31 volts , 
max 
The describing function method predicts a stable limit cycle by 


applying to Figure III-B-3.1 the same arguments given in Chapter I, 


Section I-A-3. 
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III-B-4 Analog Computer Simulation. 

ra The analog simulation of the given nonlinear system is shown in 
Figure [II-B-4.1. The coefficient pot settings with associated re- 
sistances and capacitances for real time and real magnitude scaling are 
given in Table III-B-4.1. 


TABLE ILI-B-4.1. 
Pot Setting Associated Elements 
E For a Step Input of E, Volts 


i 
i 
1 160 R, = l Meg, R-, = Ll Meg 
1.0 















fl 


1 Meg, Rey = 1 Meg 







v2) 
ii 








(See discussion) Meg 


Meg, 





-09 Cinitial 
SEEING ) 


> Volt relay voltage 
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The back-biasing pots "a'' and pot a, were set in the same way 


3 
described in Section III-A-4. The static characteristic curve is also 
the same as shown in Figure I[II-A-4.1. 

The system was first excited with a 5 volt step input signal and 
was found to build up to the limit cycle shown in recorder trace IIi- 
B-4.1. The system was then excited with a 20 volt step signal and 
found to go into the same limit cycle as shown in recorder trace I[II-B- 
4.2. From these two recorder traces the limit cycle is seen to be stable 


with angular frequency of 2.22 rad/sec and an amplitude of 11.7 Volts. 


The phase portrait of the system is shown in Figure I[II-B-4.2. 
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IYI-B-5. Comparison of the Four Metheds in Precicting Limit Cycle 


the System Shown in Figure {{(1-8. 


Method Anguler Frequency 
rad/sec 


Mitrovic's 
Root Locus 


Describing Function 


Analog Simulation 


W 





ft, 


oe 


ry 





CHAPTER IV 
RELAY WITH DEAD ZONE IN ERROR CHANNEL 


IV-A. Loop Transfer Function with No Zeros. 


. The block diagram of system chosen i& shown in Figure IV-A. 


‘¢° | 
? 
| 











5O 


x = 
sie 5540.5 YS +1 St2NS+4 (SHB) 





Figure IV-A, 


Relay with Dead Zone in Error Channel Loop 


Transfer Function with No Zeros 


IV-A-l Mitrovic's Method 


The characteristic equation of the given system is, 


ee + I. 5° +7758" = pC 124s5* + 3225+50N =O. 
(IV-Arl1.1) 


Assuming the last two coefficients variable, the characteristic 
equation becomes 


Berens) 7). ent 15 5 6 chi Aces B,s +Bo=0 (IV-A-1,2) 
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where 


B= 50 N (1V-A-1. 3) 
and 
B, = 32 (1V-A-1.4) 
From Appendix A, Mitrovic's equations for 5. and By are, 


Bo=-L24 un O,(S) + SF Wn? ol 4) + 77.5 wn Gs) 


(IV-A-1.5 
+ 15,5 Wn? Gal) + Wn? Pel ¥) | 
and 
a 
B, = 124 Waa(¥) +155 wn Gal ¥) +77. Stn By (4) 
4 5 (IV-A-1.6) 
+15 5 Won Pcl4) + Wn? Ge (4) - 
Substitution of the values of the % functions from Appendix A for 
(4 = 0 gives the parametric equations of the stability curve as, 
B= 124 ro. 77,5 Wy, ee (1V-A-1.7) 
and 
= 155 Wn : G : 
B, = 7 =) Ios (One (1V-A-1.8) 


Table IV-A-1.1 gives the points of stability curve plotted in Figure 
IV-A-1.1. These points are calculated from Equation [V-A-1.7 and {V~A-1.8 
for various Wy. 

From Equations IV-A-1.3 and IV-A-1.4, the M-point locus is seen to 
be a straight line parallel to the By axis with a maximum of 8, = 50. 

From Figure [V-A-1.1, the intersection of the M-point locus with the 
stability curve is seen to occur at roughly = 0.46 rad/sec. To more 
accurately determine the intersection more points in the area of the 


intersection are calculated in Table IV-A-1.2. 
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TABLE IV-A-1.1 





TABLE I[V-A-1.2 





By interpolation of Table IV-A-1.2 using Equation IV-A-1.4, the 


intersection and thus the limit cycle occurs at 


tt 


Un 0.4607 rad/sec 


and 


B = 22.84 , 


LS 





From Equation IV-A-1.3, the "average" nonlinear gain over a cycle 
is determined as, 
N= 0.457 


ve 


from which the approximate amplitude of the limit cycle is determined to 


ec E = 10.95 Volts, 


max 
Mitrovic's method predicts a stable limit cycle by the following 
reasoning. The input signal must be larger than the dead zone for the 
system to respond at all. Once the input signal becomes larger than the 
dead zone, the characteristic of the relay is for practical purposes 
similar to saturation. Thus the same arguments given in Chapter I, 


Section I-A-1 can be applied to Figure IV-A-1.1. 
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IV~A-1.2 Root Locus Method. 
Substituting S = JW in the characteristic Equation IV-A-1.1 and 
evaluating the higher powers of } gives, 


Ey an ies = Gio = sees [D2 pe eee +5ON =O, 
(IV-A-2.1) 


Requiring that the real and imaginary parts of Equation IV-A-2.1 yo 
to zero independently, yields the following two equations, 


= re ae 1g 
ae 77.5 uy 4 ey tOON =O TV-A-2.,2) 


and 
= iS SiGe = (1¥-A-2.3) 
lore bs ~\I5S5w Pyo2e = OO ‘ ye 
Dividing Equation IV-A-2.3 by |i) and using the quadratic formela 
to solve for W : yields, 
7a 
Q) = 0.210899 


as the only meaningful root. Thus the frequency of the limit cycle is 


JW = 0.4592 rad/sec. 
° * Z . o 5 “ ray =) ’ 
Substituting &) back into Equation [V-A-2.2 and solving for N 
as the "average" gain of the nonlinearity, gives 
N = 0.4543 
from which the approximate amplitude is found to te 
& = 11.01 voits., 
m 
A graphical presentation of the root lIccus solution is shown in Figure 
IV-A-2.1. 


1 > 


The root locus method predicts a stable limit cycle by the following 
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arguments. For an input signal larger than the dead zone, the relay 
characteristic is similar to saturation and the arguments given in 


Chapter I, Section I-A-2 can be applied to Figure IV-A-2.1. 
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Figure [V-A-2,.1. 
Root Locus Solution for Relay with Dead Zone 
in Error Channel, Loop Transfer Function with 


No Zeros 


Lig 





IV-A-3. Describing Function Method. 


The describing function for a relay with cead zone is 


ma AV \| fs - O 
Sp era, | (e— \ ee 


\V is the signal controlled by the relay 


(IV-A-3. 1) 
where 


Emaps the relay input signal 
cq) is the total dead zone minus to plus. 


The linear loop transfer function for the given system in frequency 


response is 


ee ee ee 
G (Jw) = JU 2IW+) YIW+ 1)(0,5.9W+1)(0,25Jud+41)(0,125 Jw +1) ’ 


Table IV-A-3.1 gives the magnitude and phase of Equation [fV¥-a-3.2 for 
the various frequencies that are used to plot the linear loop transfer 
function on the gain-phase plane of Figure [V-A-3.1. These values are 


obtained from a Bode diagram plot of Equation IV-A-3.2. 


TABLE IV-A-3.1 





Table IV-A-3.2 gives the values of 1/G, which are plotted in 


Figure I[V-A-3.1. These values are obtained from Equation IV-A-3.1. 
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TABLE Iv-A-~3.¢ 


Oi 5043 


0.4823 


0.4247 


0.4026 


i 





oN 
3 
ew 


The describing function method predicts a stabie iimic cyc.a tL, 


fi 


the toilowing arguments. For an input signa: grester chau che cezl 
zone, the relay characteristic is similar te saturaticn and che argu- 


ments given in Chapter I, Section T-A-3 can de applied to Figure [y-2-3.1, 


r= 
RQ 
be 


@aaA 8 @ basa 
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IV-A-4 Analog Computer Simulation 


The analog computer simulation of the given nonifmear syttem te 
shown in Figure [V-A-4.1. The coefficient ret settunge wit wesrvceisired 
resistances andé capacitances for real time and real magnitude soslivg 
are given in Table IV-A-4.1. 


TABLE [V~-A-4.1 
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To accurately set the back-biasing pote ‘a' and “a " a 


gute 


Ye Wave 


t 
(€p) 


generator was connected to the input cf the relay simusation. 
An %-Y recorder was then calibrated to ¢ wolts per inch om each axis. 


The 


tv 


nput and output of the relay simulator were them ccmnected to the 


; : : a, ae 
respective X and Y axis terminals. The pocs "a" and ‘a 


at 


justed until the X - Y recorder plot closely approximated the static 


characteristic of the chosen relay with dead Zore. The Static character- 


c 


% 


istic obrained by the X - Y recorder is shown in Figure Pv-A-4.2., 


tes 


the system was first excited by a 10 volt step input signal and the 
system was found to build up into a limit cycle 4s Shown in recorder trac: 
IV-A-4.1. The system was then excited by a 20 volt step inpet signal anc 
found to go into the same limit cycle 4s shown by recorcer crace [V-aA-G.-. 
From these recorder traces the limit cycle was found to have an anguler 
frequency of 0.455 rad/sec and an amplitude of 14.3 voits 


a 


The phase portrait of the system is shown in Figure ivV-A-4.3. 
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IV-A-5 Comparison of the Four Methods in Cirit cy-bes Pleci.eion fae 


a Nonlinear System of Figure IV~A. 
Method 


Mitrovic’s 
Root Locus 


Describing Function 





Analog Simulation 











IV-B Loop Transfer Function with One Zero. 


The block diagram of the system chosen is shown in Figure IV-B. 
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Figure IV-B 


Relay with Dead Zone in Error Channel Loop 


Transfer Function with One Zero 
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TV-B-1 Mitrovic's Method. 


The characteristic equation of the given system is 


61+ 12.55 +388 + (16 +5O0ON)s + SOON =O. 


(IV-B-1.1) 


Assuming the last two coefficients variable, the characteristic 


equation becomes, 


614 12,.6824 3857 +5,S+ Bo=o 


where 
B = 500 N 
O 
and 
By = 16 + 500 N a 
From Appendix A, Mitrovic's equations for Bo and By are 


Ba =—L 38 en $14) + 12.5 Lan Go (4) + or $e(4) 


and 


By = 38 urnbs (6) +12. lon G3(b)t lon®h4 C4) - 


(IV-B-1. 2) 


(IV-B-1.3) 


(IV-B-1.4) 


(IV-B-1.5) 


(IV-B-1.6) 


Substitution of the values of the @ functions from Appendix A for 


% = 0 gives the parametric equations of the stability curve as, 


a = 38 We + Te 
and 


2 
By = 2 25 WD, 


(1V-B-1.7) 


(IV-B-1.8) 


At the intersection of the M-point locus and the stability curve, 


equation {[V-B-1.3 equals Equation IV-B-1.7 and Equation IV-B-1.4 equals 


Equation [V-B-1.8, thus 
500 N = 38 Ww, + oe 


and 


16 + SOON = 12.5 Wa- 
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(1V-321.9) 


(IV-B-1.10) 





Substituting Equation IV-B-1.9 into Equation IV-B-1.10 and combining 
terms yields 
4 — 
Wy, ~’is Wy —-IG=O . (VS = a) 
Using the quadratic formula to solve Equation IV-B-1.11 for W 
yields 


Ww” = 26.11 


from which the frequency of the limit cycle is determined to be 


OQ), = 5.11 rad/sec. 


h 
| 2 
Substituting (y, back into Equation IV-B-1.9 the "average" gain of 
the nonlinearity over a cycle is determined as 
Ne 0.6208 
from which the approximate amplitude of the limit cycle is determined 
to be 
£ = 8.054 volts. 
max 
A graphical presentation of the intersection is sketched in Figure 
1V-B-1.1. 


Mitrovic’s Method predicts a stable limit cycle by applying the 


Same arguments given in Section IV-A-l to Figure IV-B-1.1. 
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Figure IV-B-1.1 


Graphical Solution for Relay with Dead Zone in Error 
Channel, Loop Transfer Function with One Zero 
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IV-B-2 Root Locus Method 
Substituting S = JW in the characteristic Equation IV-B-1.1 


and evaluating the higher powers of J gives, 


1-312. WS -3F wr FICO +5OONW + SOON=O. (1V-B-2.1) 


Requiring that the real and imaginary parts of Equation IV-B-2.1 


go to zero independently yields the following equations, 


Ly 4 = BipwES SOG IN EG (1V-B-2.2) 
and 
Aen bss NONE scores NED =O (1V-B-2.3) 
Eliminating N between equations IV-B-2.2 and IV-B-2.3 gives, 
Ww i DERIK Sb eae (1V-B-2.4) 
Using the quadratic formula to solve Equation IV-B-2.4 onto 
gives, 
wy = 2G 


from which the frequency of the limit cycle is determined to be 
lw = 5.11 rad/sec. 
Substituting mye back into Equation IV-B-2.2 gives, the "average" 


gain of the nonlinearity over a cycle as 


N = 0.6208 
ave 
from which the approximate amplitude of the limit cycle is determined 


to be 


é = 8.054 velts. 
max 
A graphical presentation of the root locus solution is sketched 
in Figure I[V-B-2.1. 
The root locus method predicts a stable limit cycle by applying the 


arguments given in Section IV-A-1 to Figure IV-B-2.1. 
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Ly 


LINEAR SYSTEM CLOSEO LOOP RooT- Cr. 


<-LIMIT CYCLE 
W= 5,10 ed/see 
Naye = 0.6208 


Figure IV-B-2.1 
Root Locus Solution for Relay with Dead Zone in 


Error Channel Loop Transfer Function with One Zero 
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{V-B-3 Describing Function Method. 
The describing function for a relay with dead zone is given in 
Equation IV-A-3.1. The linear loop transfer function for the given 


system in frequency response form is 


iuj=— —_ 3.25 OW) PS. 
G Ww) JW(2)Wt! KO, 255W F1)(0.{25) w+!) ° eee 


Table [V-B-3.1 gives the magnitude and phase of Equation IV-B-3.1 for 
the various frequencies used to plot the linear loop transfer function on 
the gain-phase plane of Figure IV-B-3.1. These values are obtained from 


A Bode diagram plot of Equation IV-B-3.1. 


TABLE IV-B-3.1 


OO) |G ow) | AS JG w) deg 


-176 





Table IV-B-3.2 gives the values of 1/G, which are plotted in Figure 


IV-B-3.1. These valves are determined from Equation IV-A-3.1. 


TABLE IV-B-3.2 
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Inspection of Figure IV-B-3.1 shows the intersection and thus the 


limit cycle to occur at 


O) = 5.1 rad/sec. 
and 
Se = 8.95 Volts. 
max 


The describing function method predicts a stable limit cycle by 


applying the arguments given in Section IV-A-3 to Figure IV-B-3.1. 


138 














COO rr jou jucea J 6G 
© 2860 (SUE 8 CORRE IDSC AHSPOCATTY 1 BODE 3OS0 SHERTORERE CORPS ORC! INGER REE BO OEE BUSES BUS bE Oe DEER seen OOOO 00000 CGB00 GEROD SERRE S00R VOERUAUEOE BORER BOS88 GPO BB O8 SERBS ESOSE DCSRERBoE: {i IOeeeOws Feeweeeea 
498 80 800 G0 G0 000 DOGL0 O08 0OORE 0B OBR 08085 GRRL TPSSHD 0OER DOSE BUS 00s288S CRORE SOEs OeOeR Dene DOSS CSG0 00 SSSE8 FG 000 SERSBSEB00D CECE OGRE) PAROD COURT POOOR RP ESs BE ORE BO OSE CRCOR REESE POSSE EREE! i 11S DEeeOs eeeeeeeeeO 
SEOEE SHS 88 FOE Ss Bae OE SHS Se FHS 28 FSONS COREE PRVSEE ~ ANE POET ASS ESE CESES BES SS BER BESE ESE PEERS BESS AEE OE SE BOO ES PRESSES ES SESE ESBEEE sees SEeEeeeeee GORGE DAGH8 CUDEE FEES 0 OER SSSR PARSTSSRT | 1108 OBeE8 Gees eeees 
SS FES OSHS ROETR AG DOE SASS OOORE FSUOT CHSPE FAUHE CS ~ 15 SHOT A OC ASP t CONSE BOS BE OEE BE BS BES POSE SOS OE BOP REESE Oe BEERS ES Oks CREE Bes Bea eEBeeR B00 G0 D0 25 GOS GSS O8S0 RODS OTOR0 BAO VHIOBEBE PASSA SOO.. | 1 1OVOBSCS eae Beea0 
DEEDS S2ORG CRORES PLDD PRORTRAOME FOORPOOCSSD FODSE EEE EEE EEE eee eC LLL CLCLELLOULLLLT te 
ees ss geeene! Seng Feens senes Sees eeessees DOO 80 UO800 GU800 SEOET SERED SON 0B OSE S0 Fees ees e SRSCC SESE SECS ERR EE SORE RROE eRe: ARERR ERBGRESCE | (CORP EOSS BOSE Bee ee 
DESEO FUSE JE SOL CEBSR ARGS ROOTE IP RCODEPES BODAe Gr eepeteeses 


OOS SSBEE) De. ~ SERCO SECO SELES PRU SS SORS SOS OR0 DHGSH f OORE FPORD CABS PERSE COUR SSRSOSERS) ROGET BECRE CESS BOSRE Bee Reese Sap— HESES SBSSS OSREL 17 100888 
BRS RS SS SHE CV SES SSTEs FSTEN SSS8 8 246 OS PSESR PARSE OBO RAGOE! BRP E PSBOE CES 88 BEBE BSE SS PORE EEE SS POSES BSE Ss Fees eee BRS CORSE DO CSRAURRDR! CORRE GREE FOCUS CROES SERED E SARE OS (O. TORDSSSORASESE 7 1 IDOSPOPs CEEeD PeBA 
Fd 300 OWES CLSORESST CUED BEOUE LOSS SPIED LeDeR CTL S0D800/ 236 . JBBOG COB00 COLE EF SOU RE BUSS SOUS BS PEOOU GOURD SPUOE FRURE POSES PODER GORE REASE) FAT S5 POSE BUSS0 SOPs BOSS eBEs GE BBEGS SBORGSREE. | i TUDES RTE BEBE BeeeE 
table Ieee Pee Tir rT Tier Pier iT ITiltilil) emi l lll miirriiitiiiriittiiiitt tw 60S Oe CO GSS SE OM GG 0 CESST BOS HE FOS RS C0 SE8 CORRS ROO! PRORS BA BES PSO ES FOOSE BEBE E OSE E AS oc DEERE CECE ESERE. |) tS RR RE Bee ee 
SUSSS CSL SE CRSRS FB OSS SESRK FES81 O08 CS SSP OS BSARSE Roe BSGGE! VO * POSED PORES CERES Pee Oe ee eee ee BBOGS FOCOU DOGS S CORRS SROR) LECCADOSAT BOARS LARA SARS SRR ORS OS) \.t PROSE SERCC SBOE. | 1 180 CBOSS BOSS eeee 
27080 GES BF BOER TTAPO QORBTDPBAE TADOT ABB pA fF 22 -- [UDR GSC F DURRT RS DHL DRROA CASES FRESE FOR BE OES SR EERE eee SROQP “Ag Gi SHSRRECCRE OHSS POD O8 FORTE OOSR8 OE:1Gd SERS COSRSSSEE! (1 ISRRSOHS BOBeeeBAeE 


HESES SSCS BPE RS OSD HO OOORO CEDOE DBBe8 FOBeB ae ee anene Bb OBESD0 GS G00 00800 EUS 08008 8 O00) DBO R02 S 0000 PSUR 0 CS000 FBGER FUSOT OEP 1c~ OD BODRE FEOee PESPPAAABs BERD BBEGS SERRORSr VS SEBS ESEE F12GCCRREE BRRee eee 
SOCGSG ED ESC PERG 0 08000 CORRE SCORE DBRRSOSRRS BBR888 COGEP) 10 (BSSESSE RE PARTE BOON 8 BEBE OSES CSE PES Es BOOED SAREE BES SSGGA <2" ...3 EERES Been OO OE S19 08 CESE8 BESS8eeee 
LESS CES JESSE IOG0 0 GHOSR FOOT POSPOBHS HS TOSS PS 4: PRERT!) Fi BO PORES SESS PERSE ESE E POSE eee BE BSH SPEDE SHPS80 PSS SE FOSRR CESTE SERRRFaSE) GHRRT DAOOE CODCOD OHO FEBOR RSENS 88 Oss OO. -*SBESREEOO | 1 CR RSeRE BEeeeeePee 
SEO EES SS Sl COSRE SESE ks SERGE HHS OE FESSRREEHE BOSPHS USS PORES! BE. » DSOSTRRERE BEES BOOBS CES OU eee s BSE Bees eee eee ODGO0 00000 OBBS0 SE0CC FOSR CCE HE CUD GE BORE FRELESE SSE OREO TERRES CPSESESOA . ) 1IQUUSSOT BOeeasBeBes 
UGE BO OS SSS GHSS0 FODEH CPOS20 0558 ICBSTHOSBOE FBSASE 1) (DO ESOSS1 CHT PDEOU CE ESE DBEEG EBREO PPP DE CODE SE DES SEE es eee bwEEE GEG 00 DEERE FR0S0 SEER) SUBS RS SE BE POSSU ESE RE FSS TROR 28 er—| — PPUBE BPSSASBOM i ADORE GOSS ene0 
PERSE SSB Hs GEASS SSERS 68S E0 BONE 6 OREO SORE DEERE — ¢ BEEP BOY) 6BORE SESE BSS ES Ese eS SESE BEE es BSE Es ee BESS 0 SSS BE SOO R0 SRBSE SESHS GES) P00 C8 COREE FER O08 FUGUE PERSE SE ES8 OS. VOR (AS PSSSS RRS 13108 BSee8 Seesseeeen 
£P GUE FOGOT SSS SE S08 B88 CBOSS E5700 D0RRS CSRS BDAES OBese LORRE) POce CODES REDS CORSE BEE Se Oeeee Pees eee QE USE GOR G0 CORRS GUC SEOD SN CECRS BRED) POORUDORSE URE RUR OSHS PEROC SEES OS ClO 4 y PUBS EAOE | 1 ISR OEes Bees eeses 
TOOL OBC ED CBSBOSCO8H02 QOBRS00 000 OON08 BERS BDABB O20 OGRE! GOe:, CBGDS FEU U8 DED SE GU GEE BOSSE FSBG8 PEE OS COUUD SHROUD CR ESS SAUOD BESS SORR8 BORE) BERS OSES FRB Peee8 88 FOeee eEees ee aU BSB. BERBER EREED FERES CREE Bese eeeses 
1ODER COGS 0 CAG E8 F88U8 FOGG 208281 DODO POD BBeBeD © OOO80) 6Ces DEBS 0 20000 DOORS FEU SE DER SE BOREE DOR RE COG H0 CED ET CAPRE SHED ASRS BERRA CREED BORCS DOCS REPRE E Os OTeee BEEBE | Tae ooo 
36000 COOCE Of BS32 E880 OBBBSBEDOCS IAGO 0 TERS BBB eD COGB00! 0O*5 OBER 0 OSE BE GE00 0 GR ESE E08 PERSE SOUT DOU UU DR UEN CUURE CEOS PROSEE BARRO ESS8t FRCS SC 88D SUOES ES O88 OEOSEGUERERR cu. Ga. Ui 8 BOSE BREE PASSE ESSE eee asee8 


BE SSE PSS! SPER SSS RE APSSHSGh8 ISHCHS ARH ASSREEt 5) BOSE Ot OR: = BERET PORES BSCR EEE eee 2016 BESS 2 SO SSl TARO 8 FSS CERES SS SESS H AREAS FROGS CS OOR PARES OO B88 ORB BEES OR Be SORES SREP THEE | |} Ieee ees Bese eee 
SEGRBSOGERE SE888 GDRS0 DOBEDERERDAIBAH QOSS0 GBRRRR- 240004 2U0 BBLT OR ESO O80 G8 POUT CCUDE DERR GRR E0E PRE DE BOROE BORO CORES POCES SEO RS BBORSESO88 GEO8s Bebe SEGRE BESS CRORE SOBER ER .“XEE..~ 4 BOSE S BOGE 68 2S8 BERS eeees Bee 
BEGGR 8 BS S8 CERES SESER BESERR BE G8 SOROS CRR88 BERR dben es OEE bE Cc ay OOD00 00800 DEDEE BEU008 BBSRDUEOL0 DOUCET DEBUC CEE SEP BE OS COC SS DEBAD PASCO D ESE POUR SORBET BESS) FOR Oe Sees eeeee GEES EB BE SOBER HARE) | 1 118 RB ESE OOSSSGReee 
bleh eteafeetet tbat tad Nolet at aor  er eT eTrTriiIfTiiiiliiiiitilitiiitiltitt ttt Pe Le G8 G8 8 BOSSE DEES BEDE) FEGS8 DEORE FOE ES DERE FRUSS RE REE BBen JOG U,”..8 BES ERE 6) deR BESS BESS Peses 
{ER 50S S80 200 BARBER THRREOTCllD DOSS CSRSRE Gen EE est Se PO=9 bi Bb. BU EE OBE RS POUR e eee (SCOTT Ee GOS OS SS SS” CES He SESE BESS ROS HRL CORRS AER E OP“ 8s BOS BE BOSS CES OD BERR ERESE BS TEEE.  -SBSSBOESF! (1 ISRO BRE BEE 
20800 COSCOUD FO008 SEHR OO BBBOD005 1000 C0888 BBBEBB-— OOBlOV Bece GOEEECUEES GRR00 DERS0 FER O0 COS OS FOOSE GUSUU BES COS U08 BFPO SDE SR GORO8 DGEE8 O (Bw UE ORES SOE 08 POORE EROS eEeen SSS 108! (BS SESER BE ERE £4388 SESR8 Sees sees 


Seead! OME OGG 80 SUG58 SUERBR DESED GHOBE DEORE BELEU CHEER PSRSRHS SOU CHORE FOROS CO8R8 FeReD ESS (ROSE CSE S8 BEEBE Bees eseEe San. SEB! ‘S82 G87. E. BEGGS F2SSR RG s BeBe eeeee 
v2 6 GU *— DBE OU CES8 0 COURTS ONE CORON DUDES PORBR OES RE DOSES SERED BOERS PEO RE FeeeRBORa Bee SRS GR BEGHE CEOS BASES SSRES OH 28 FE) DEO 188 BOREL 013580 80888 Beeekenese 
S808 009550 258 BE SRR8C OBES COE Al OOPeDE ees SaEnem aL enG@ jc meh? BEBE88 SEES SEES eee 


EESG8 BS S58 00008 BE SS8 CES 008 0008 88 PRR SESES BOR RO8.-_— ae BR2R8E F ace VORES BEES OES O8 BEES 8 CU UED PROGR SESS POEs BORE Beek Pee R ee BES SESRRES FOR 8 FES O8 SORE CORRE BER SS REESE SER BRC GE PRSRe BEE Lf PCR RSERREPee8 
BOSBOE GPOOE POH 20 BFEDT GASSER ROTRl DEBORSROBD pe tt ey. J OOOl ND Br @c DEBE FOGED FUBUR CBS SN G8 S00 DRSES FERRE PORES SOS C8 BRED.) DUAR ALBERS GUO DEESERE FBO 40000 U. BBRROSOERT SPE eBeEe BOEF' Ht St Ub SEGRE G6 G8 SEG88 Seeeeeeeee 
BESGE VSel s SSCS S GSES (HSSRHSSEL PODOCRBSSSROCTT TH Bh COBH) PSPS BURRS EPS S Oe OE PORES BORED Pease eee ee SOS SR SSGOP NOGE8 SOSES SASS 0 BEORE PUREE SS. 8 CHORE SERRE ORERESRSEE eRe. JEES— 28) B SOROr Fede8 Sees seee8 
TECSCOM SS TORS SERS A (fF PSSSCLSS COS SSO SOARS OBSSSR «FER? 38 BS) TSE SEER P= OBOE OBESE BROS EB Ee PRESETS OES SES ER ESSE (2 ORES Pees OEE ee eEeEs BS SOTEDErE DEBAS SEDGE SERBS SERES OE BOX” AS ODEBARASA C7308 BOee8 eeee8 based 


JSGD8 G0005 GB00 000000 00025 00 050 COD HO OLR QOS RD HORE BOGXA8 08 § 4 PRBOREO0d | COBO DBCS CUES DOSED PEP DEP RBA eeDeeTBAen aH BOCO0 GOB R0 CELSO FO 2 UO SO044 FOE00 P2000 GCBOSR SU T00 GUEES BA OBEOSQe2 BBO JT ABS OORS BOeeeeeeeo 
ab 18 BEGEE COOPER BOGS OF 2. VOOR 1S BEBE S SORES BEERS CERES SEES eeese SBer BEGEE FEOS8 BGe8s Geek eeee8 
SG84 SSSR SE888. SERRE SS6! U SOS 88 BOSSE BESE RES E8E BEE SRSREs FG AG BOSSE CARES BRSOE FO Ue BES) ~ HOSS BORE BORER BEERS SERIE ESES. — tr -CREO Es F058 8 Bees Bees eee 
PN BS ee TT CSCO CO oo ooo 2B. 4 BSSS) CESPEOSEBE BSOSRRRBGE BK i108 B80 Reese eeeee Ci0BS BRR SS POSER BBSS8 
STO BEE~6 40“ S BES 6 BSUS SBEESE SOSEE SSK 1:81) DERSESE SSSSSSOSSS 66 . SR EES 0... SEES SEEES SSE eSees Bee Ie GESES BOS —~VESER Lr TD EOOORRE GBOUBAee8e 
G..=— GERBER BE Bba BEEBE Gee t , JTRS SESE BEES BOUOE BERSE Feed & ~~ 408 CRESS See eeseE GHd OBO. FRO BSS EE EE OBOE SBEEE SOSUGEET WT BS Be Veen Lt IB OO ERE Ge aeenee8 
DEGSREOCEBBSR COBBBBED OM = VOSS 80 C00 G8 BURUE FOGSD DEGR0 SLUR U BOSSE SSOOD DORR 0A CED OOs 4 BOOB ER BEEN BE OBSCURES 88 WORE) 2D OUe VO GEOR [ PUES Reed Bees sesee 
eee Oe EO OO OOOO 8... JE88 6) SESS RE EBS SERBS SRRE8 eB Lig ‘Bae QSSboRR ESL £43270 GC BS00 GUBER See8D 
SEOSS Se oe! IS SSSRGR El SASS RSOPE 1SDODOCSRE BERSE Rs ERSRRI ET: Ine aREeD! § FOREE ees CeeeR BOSES SEGUE BORER BOSG8 FES8 BOS SPR SERS CORRS RRRES CERES EES.) SERSSRESRE SEES BSESs OBK-* Bat = SHOS8 BEES 2 see Beee8 Sees sees 
¥ 9 0 BE BS POSSESS SERS RERES IO SRE Re 1 ROSES PER Sh AOE SEO SE OCS BORER REES 1 0 RO RO SS Bee eee 
SEE BS CERES PESS0 SURES C8 BEN SHES FDES SEE OR OHSRRSF—_2 PRO RS ESSE |) PSSeeBSOF< GOERS CSESA 2OP—-UEE S BV PEBE TEESE SSSSSSSRS 6 PORES EOSES CRASS EBESE FOR 8 O's VW EB OOE BERET SEEGERS es OEP. Be GEGEO GO08C (7 IDGBES00 SEBS 0200380 
IBESS DASSR OBPDTODRSD SESRDESDDS 10 DOORS BTBARRe . J] BRBBR €SREB sO SOROeeay SEBBERGEF cs B°3 0ERe dé. BE B88 Beeee BOSGEESES 8 BROS BOSE SEBCE BOSSE CHOLEES Ch) SOR 08 BRSSE BERS8 OBESE BOSe.. 88 3 L 
BOSS 82 CESS BE088 BOS08 BEERS BONES FES BOORS Os SOREN Bb 1 BRSRREY 608 SUBSE OES. SEGESGEE € 100_ E0008 SRESS SS08s Bees sees eee BGDER SPO0R2 DESE8 0URS0 BU... DOB BR POOEB SES eBBeORe SOBC’ GEShe £0—"" SEE88 BO280 FOSS8 Been eeee8 
ID SESOGLNC POSER 8 SBOd POSSD CGNGL 10550 BPReeRBeED Bear Benne0a i 198 SSG eeeeh SESE Eeeee.eeseeese 4 00000 SESR0 UBBD00CS8 0 PO808 00000 FB8GE DHURE FEBOD0 CORG0 COR OODADEBeBeeeOrT po ee Bo UOEEEE Ff DBBESER8 PSSSSSBeaB 
Se a eee TT tr mr iii iiiiifiibapiirreliitiiiiiigiilititiitittt tert GBC 08 SEOS0 SSPER2 SORES BECCRTGRD2 486 ~y BRGSL #2 1 7BSUR08 BEDDLSBRL0 
JS 5S 00000" D668 BRE OE00 22005 PLOEL (ONL 5 BEAD POSER LED BS DOOl = "GE DOSES ER O8 BEER BEDE E OES 2D O80" SOB SBES SHS SERFS 80208 0 CPR S0 02000 GRODSTSE08 PCC DDSRSCS DOBRO GSBORSRSGEESEES 2Het) Hae ay bid +--+ See enna tt 
JO508 00S) 1008300000 0SB00 COOGE (LODSR DSBS DOGGOREI BE ROR8R: 100 SEBSROndd “1800 Fas PeP=2BESR. QO SBERR VODCOERRN & PBURR0R20000 FOCRSRGRCU HUBER ER RED FBR FO8SR Feo at N set ee FICSESSeeB 
O08 PS SSS REEFS Oa BS BOSS (BbE ISON) SRRRRESEEE CHAAR SASRS £ CARS SERERE BASSE EORS PORE 8 OE388 COREE BOSSE SABES Ss hb 8. "Tee 86 BOP ee eee Cee eee 
SBE BS BEE GS BEB SE SEEBRE BREESE EDOELS £6 BEES EES OEE EBs ce 38 _——- S285 Ba *ESSRS. ON VEE A) Cees BEE BRO SS CSRS 6 SES 28 SSRs ee EBORSEEES6 SERGE RRRER ESSERE! ws ASeeE RE ORR Re 12d Rees 
SES S8 CES E8 SES 88 COS FE CERES BONNE FB BOS BASES OA SH POSSE PHOS UAL. 18 SHUNT ASSES Bar ..w.. ~SRe 48h 408)? 2 2OS.. * AOR OE SAGES ESSE £ BOOS BEES Beeee Beeee BD _. BS SOD BORO8 PODER OEBBEORLEUE PEBDOSANs.. SERDASEBE Tt 1DeaBeae 
BEEBE S BEBE SEGRE PORES SRSRROSHRE CHSRS ROSE OS, Us SROES SEBS. 48 2ESESROEEE Se Ee, 14 CRORE CORSE BO > \S UO EBe SORES BOS0 F OSSES BAGSR SRSSSRSSEE TF EX) CRESS RARER ESRRTRRREEREF WP WeeReEs te) Se Ce 
SOSS8 £ BS BE S808 SOS SS FEROS FOOSE CSOST BERO8 E> “~ SRRR8 GOR O86 Mo WS OSES BESS8 Balt CREBe. WE REOEe Be.~  SS08s FOS0e Bees OSes Bees eee sees mer} tt SEERSSCRSREEREERSR 10 8 OB UN ORi ee Ber cheese Bees 
JCOCCCSSEL SEG 80 GOG80 0O0C 080000 02200 BODES SPSSHOT DEN BROBD>? 7 DOSDSSESOS88 fi -71h8ne88 Ht BBRCO SRE OS BOSRR SR fi'!) TOCRPRODEE GBAOSDASes TT SGGG0 BERGE SEER SREBREE “ 6 Be. 'BOC_~ SES> BESTE CSSER Reee8 Bees eeee8 
SESSRESSRE PVT ERE Feo S SERRE A DSCE (RODE OS EBS PARE PARES SEER OPS Rs ee 
See Be ee eee ee eee eee ir ee | « IGE GEC Sage tO Eos OCU SE PbGSSE DUS TREES EY CP GU BED PRBS SDD PEPPER EE rE ‘SSBRE GEER EBDOL Lt SDSS O08 OBE BREDOBE 
3000005000 000 [BG 0000 COB88 BOODF 1000 8 O00 88 BEER BOG B0 FROR2 DERE) CERO8 BER01 FO00bEn0 eer Yar =" @COB8BB feetnanean GSRRRR2GR00 DEBBOBEL. 4 FGem<BOSOE SHORE ERSERSBEESEE ‘ENGR. year BSSSR SRHSL 1 ISHS BORE SSB eeeR8 
HOU OSH HS GERBOBOORL GORLO OB BSC e8eO8ER ae TEESE BSSES ESSER SESE Esser s e% BSBOG OGRE FT 138 BORGS BSSREBBEES 
DBO F05B0 DOES G8 BHU GHOST HDD D0 LODDE BES DR SOLS DPOORS BOOeRleseD SOS0Re000( , Weeeees/ noon) weeee BOGaSED Rt et ff eo fo on ‘Sanne BOS O85 BOSSE BORE Se ie) BE. UBER. CREE SEOSE SESS CORSE SEeEe seeks 


= OSB! DOBBES SEE SS RSS SSDSORD i, G55 4. 55852 SSC" S088 §se88 8880868 Beea8 BGR88 


BCOBS B30 D5 PRS SRCT EDO OOSE20 03 BRU PRDDBBRRBAR CBRORDDA th +H GBGE0 GRE BER0R8e sIBED Ht K I Seat wasere oe BOQRE BABY - U DEBE ES EES BERS BERS SE FOS) DOBO OS EROS BEBO BOSSE VAaBR. OSC“ BEF oo USP — BEBEE CLS08 GEEEE BEeE8 wees 
SES SSRHH SS BOSSE FFO88 SSOTTRERSE FTSSRRPRER BAPAROES SS ROS SSOES GER SESS RRSRE EEF. € BOs) BENE.) t LAER e~TEe BABS EeR8 . eee Bee OOOO OEEST ISGUC COR 00 GESTS E0000 O88 DO OBORL Bers BED a8 URBu. DESR0 1 IOS ESO SS GSSRROeeAD 
BOSCO COORG 00G 520 24080 08000 BROS 100 88 CORRS BADER DUEBd. ZO BDRT ET BBU02REORE BORER | ND S685 Sens je po + BN GOR BESO8 Be SRSRE FORGE PROCES SASE FREES BOCES BUSES POSRE SRSSRSEERE CE Nn Abd” BES... CREEs Cees Bees seen eeees 
PSC! JEU. FS CSSSHSRSE £ SOQCRRHS OF PAROS BORES BEERS 8 OER POPPA TREE EE BERER ROSE Bee SGSR8 FEG_ * JOSS EROERn (BERR ESSAE SERRE SRREE SHER S BOSSE CEROE BOSSE BASES FERS SSSRSOR VEROs ib Boe’-ti wee Sh SEGRE P2088 BORE Bee eee 
Jf BOF (OSS COSSESS SUS GEESE OSE O01 C7 PROPOSE CEE RE ESR PSP PE SEEE Pe MOs CRUE E PERE S. Cee ee SE SES SBOPY 5. (OS COREE POSES SEES ET CESERSOEED CRORE CORES POSSE PEOES BOSSES 1 BE Cae hh’. SBS. BESeL 10008 88GR8 Beek aeeet 
be Eee eee) en iomrlrilllllli tlt ltt ibey Lia) TT TT) toa GO. 7 68 0B0 G50 BE BREESE GBEOR ODED VBS DAODORA COADARAERE PERBEE8 Sf pBUTAeD_— HBR 4 GEce0 SaEE8 PBReREeRna0 


BDCC ER EHN SPCRSEEB BRC DEGR8 BB8D0 0588 ? ] aBWJBe J S8880Ge080 0029 80D SSSR CSE SEE DER SS 80800 LESS 0 DOSRD PERRO 00008 OBSRO S88 - SREFr REO A OBO OR GSERD 01008 Oe888 SOeeeeee80 
= HH tt Gsoeeseeeeesecmecen® ar @ JBOSS SESERSeE | eESE 8 Ha SOeeeaee *s SeeGene-SeGneonneeee O@G508G000 foeelGenae BGU0DGRLR8 SESGS008/7,\pannt BES) sBar *"ESSR8 S8SGS88888 Gebetaeene 
0050 GERE0 PURER PRGLD GADRO COBDL 1IDDAD ODORS PPROTERPAD POLH8 Be8eD @’ MOSSE 88 CERES SEOP= RSS) BASF See RS BOER BEOSE Eee: SReee eeeen La | COO eet oo 


HESS REOEE 0500000 G800S 00080 10008 03080 C8000 FD000 8 PDR CESE0 0 SaaS E0080 OBSER BSS CSS DDT0t DELDU ESE. e000 BEa—’ 40080 Seeee SQS00G000 F080 0000020 25088 FRD00 PSBDEGD SOF @28 BEVSE SOF SF ORES £4888 BROER BeBe Beene 
BEBE CEE EEE EO Oe BUS” JB FU BS SF SG08 £9 BBSESRR8 Feee8eee8 
JE SSSSUSSE 18 Ul GOP ES PARES OOHRE 170 DS OO ORS COR OS DORNS FURER DEE ES SORES SRS BEES EBA Ge) ONES. E85. SEB0 SEShas CORRE BES) F NSRS8 PORES CBOSS BESSS 68" TESS SE FOEes Bees Beees eee Bi. S55 FEB. +E ROSS BOER Beles Beeee sees 
Cee S Ec eens Ce Gee Seen Gre Se BOOEs FOOSE Sense eee eee e ee SORES SE Ens SESS 0 S000 SERGE Sse. SEE re Reeh, SESSESSS) | RORERESEL.) REGS SSSSs SESE SSESS US ie) SEEEs BOSSE Rees FOES Benn, OSEEROSEER Fer seneees rsseeseeseeeaen sense 
LESS TSSSS50 80 SHELAGH SSGSR HHI 1 (ROSAS OPD GOSS 8 SO OER BERR SR ROER PORES FP ehe CBRE ER Bee oe. sees BOSS EDS S0 2000 DOG SUPT DABS &Oy “1 CHASE SSCS FRRD VERSO BOC~“SER88 BBP". 40008 GE88888e88 Feee8 eeee8 
ESSeCOSSSEt § SPSSBEOCSSEROSI6L 1) DSR SR DES SPREE GSES SS EE SESS EE CES EE 68 ESE ee 8 BES CO Oe OEE er BEEBE e RE 5 Oe REE Be ee C§HUDC CROUSTESFS SBBSR88\ BEED 408 BBE OEr > SEER 19808 BOSR0 BOCRSGGenr 
PREECE OEE EEE EEE EE Ee eee eS 
SOSS5 SOSGS GOSSSESESs CESSES GSEs COD eR GE EOS Corer cee ee eee eee Sean Sean FRESE SSESE SeeSeeEeen sEeeeense: Seees sees) €000 DOCED DB0E0 BB80E GO: {2 UBUOE GBSSSSES00 BEBBBOGB el lee eT 
Coan ee ee ea SS eee Cer Le LL LE LL) {988 BBO SEe BBS Se FOS SE OB... ~ IDE SE CORSE Bees eEeEe ee SSBEBSEBASSBE Ira RBOES CERES OOSOR SERB 28888 
ptt tee rt tk ff! Heed 5.01 SSR E B8ee8 Beees eeee beees 
Re eee eens E885 SS SS8 F588 SSE S8 FSS SS SSE S8 PSS 08 BEE SE CESS RSS EEE CESE8 CES E8 CRESS FS SEs FEES FESS CEES CESS. OSS8 SESS8 BOSSE SOSEE SES! ASSESS SESE CASES BSSSSSSs | 8 O8.U es: sie monn PRREE SECEDE SBEASeBeee 
pl tlt ee EEE EE TOO Le LL Sev Cire Li r LE Lee “Lei 
Ree Saree B58 SESS E ES BE R08 SSSN8 £6 805 8008 FOSS SERENE BEE CERES E08 SERN s CRONE PESOS BSE Es SESE BESS PORES CEES BOSSs COONS GOS ES CEES BEES COR V REESE S PESOS BOSE s CESES SESEs COSEs SOEs Bees sees Oe eee ees eee eens 
feet ere ete eee lead eee rep etna et eee en abe) alee illo oe EL iti Lec LLL LL ELLLLLECLLULACLCO LLCO LTT LLLUELCCOLLLD BELLE LLELICCLLIILE LIE 
pad ae td a kde ene ee relegate ee el ee eee ee LLL rt) LL el tt tr LC LELCLLLPCLLCCLLELCLLLELLLLIELCLLCDLOLCLLLDRELECLELLCCLELLLEICLI La 
ce ed etl tet oe) od ot LL ot) tC) bbe LE LELECLLOCLLCLLUCEMLLICLLLLLLDLGICL OL GERELTIEL EDEL ELLIE LL tat 
8G 00005 0508080000 C0008 80000 135900 20500 2Ub0 080885 8020200080 F250 0002020 BBERBE 02500 G0 DEF GE O85 POSS SEERSE PPPS E6 GHeS0eROR0 Besse eeeBEe BEES SEO 00 SERRBRR SESE BERSD OBO UE GHOSE GOSS GASSR POEs .GUESSSOR0 FER DERe8RR0 
SCO’ SSS COM SSC SS 85 COSORSLOL LVR AR ESRS BOSE BBOEE BEC TE BEBE SOD O0 8 BOR BO FSS E RPE e SPs OS ORE ee PR es BER RS BOSS Bee eee BSEE EGU SESGUSRROHS FORGE CUSHE GEREEDSSSS PASRSRSLUl COCA SSSE FBAROURRGR0 
S000 FOSEG GOSSn BOGGS GEOR RGGGES LOGE CTOEOS Cote Sree ee eee eee CORR SSS SEs FESS R SSS Es FOES SOSEs FESES ROESE RESR EERO BOSS SSS SES FESES BE EEE CES ES SESEE SSCS E SESS SESE SSeS SESS SSees Seees esas esas eee Seeee sees 
Meee Ene BOS Se OSS GSGSe G SSEs COS ES CO ree eee eee eee FOR Ss SESS S SRS Es BESS E SESS BORO SESE BOSSE BORE S RE RER SRS 8s BOE S ESE SEs SES S8 FOSS E CERES SSS es SSS E SESS SESE BSSEs SESEs CRS s Sees eeees Se ee8 Fenes eee ee eee H 
SESE DSSS 58 FOS BS BS S85 SSSS0 FSSE 0 FON SS P8 SRE SURES ORE DE COS SE OSS OG GOSS GS DERE GEOR CORES CREED BODES POSES CREE BEERS OBE BS OE OER OE SS 8 BBO SR ESRB ES PERSE Be BE OR Ee OO ees BO hee Beek eee oo 
ee Ee OLLIE LLL LELIELILLLE LILLE LLL LL 
eee eee eterno eee ee eect CLL LLL CLICLLLLLLLLCACLLLCL LD Se CLL LCL 
pt eee TTT TTT TOL TLL LLL LT. LILLIE TILL LILI T_ 
Nae ee ee ee See CSS eG ees eee eee eee SSE DS SEE Ss SE08s SE008 Sabes FOS EE SSRs FES RG RESEs PERSE SES S8 SESE SCHR SE BOSS BOS Ee SESE SSS Ss CORSE BSEEs SOS SS SS eSs CSRs BESe PRS FSRes SESS eeees es eee eens eeeee sees 
eke aNd de es) a ede te tel ee en I eet) ie Lice Lie LeeLee LLC LLL OPELLLURELLLLLELCLLLL LILLE LLELLLLLL LCL ELIE ELLE LILLLLLILDLLD 
SSR SORES COs SECO ROE BOO SSE SOEs CO CSe COE re ee ee eee cee Oe SB EE SER On CR ns FS SES CERES SESS S CESSES ESSk FON Oe SESE ESSE Se POSES BSS s SESS ESSE Ss BOE Ss BSS Es SOS EE Sees CRESS Sse Ns PEEes Besse Fees eee es Feeee eeees sees eeees 
GSE 1 GG RS CSCSS SS CSR BE IBSBSRRGL f } . CURCPHES ORORE CERES BES SD FE DEE DOSER CE BEE SESE E OBESE PCC e eee ee OBC GE CSES800 DOCOTBSHETN PORRE DEORE DURES DEORE BUSSE GOOG PPORD COBES COSTS COOES PORSBERSCE 000088002 88 Benes beee0 
S006 SS056 COGS SSSE Rs GOSS e GOES: FOGG CECE s Se eee ee ee eee ees 208 SSSR S88 BE SES SEE Ss CERES SESS FOES SESES CERES FSRRs PERE SRSES FOSS BSS88 FSS ES BESS E SESE SE Ses BOSE Sees BEES SeSes See ee sees See ee sees ceeee eee 
Dee CES ne CESS E DOSES Sane Sn Ee Fee Se eee ee eee OOOO BER Ss SESE SSR SSE SERS R SRS Es FOES ESESRs FROGS ROR SS POSER SESS SSESE SERS S SESES CES E CSS S Sees CSRS SSeS SESS CSEEs SESES CSEes Seees eeees esses eeees sees ces H 
OOS Se SSG CORES R SOEs Benne eee eee Se SOORE SSS85 FOS08 SSS E8 SOS Rs FE SSS FSO 8 8 FES Re FOOSE SESE RES Se SESS SEs CES es FE es POSSE SORES ROSS S SES FSSNG FO OEe COR e SO See OSCE eee es BOSSE Benes Fees eeees Beeeeseees seeee BBB 
ADDERS 70 DOL CRBESESSELl DBSPPOSOCT LUORSOSOHE OSSSDRORED FOP A EOD 000s 00050 COGDE C0DGE SOBRE BB0DE GRD0RDS0G0 DEGES OS B00 GBE00 CGR SS FOS SU C8880 (LOGS DS0 080 FEECR SU800 CURBS FESR GAOSE CORES SPRDSERPRC COCSR FPSE8 BeBe seeee 
Ree eee eee ee eer a Etre eee e eee eee eee ee SORES SEES SESSS SES ns FOS S8R0E85 SERSU EERE FE SES SSE88 SRRES ROSES BEERS BSS Re BOSSE BESSs SES SEESEEs BESS ReneS SOSRRSSEEs Benes Senses Feesesenes seees senes seeessaees 
SSCS S BOGS COSESSSGSS CESS SO On: Fee SESE Scere eee ee eee ae eee GOSS SE S88 FERS FES85 FERSe FS Res FEES FESR SES RS SURGE ROSES RESEE CERES CERES CORES SEES ESSE ESE es Pees Se ees POSSE SRE es Sees sees sees eeees Seeee se as 
SSESCE2E0 CSS ES ROSAS IBSPSERAN! 1 1 DROS DE OBA O8 SQPST PERSD SHEDS FESS CIDGE LAH E SOOO NT DOBEE (O88 FRESE SUES SEs BBE CAE Owe Ce eee Oe ee He He ES 
BESF SOS SS SES SS OF SES PSSSSE SESE FESS REESE POSES CSSE EF PSSES CESS S COS SS OES ES CER EE BEBE PEERS CROSS CERES COR OS PORE SR SSG BESS EES EE PERE BSE ES CORES CRESS OSE Ss BS Ose Bees Bees EES e BOSS Beeesseess Seeee eeees eee a 
SC DBD 2BSS0 SRE S8 SOG88 FRSSER OF 1s. DEO SS FORE O BS OEE 88 CODE CERES BEER BREE ee Oe OBE es Be Pe Bees OOO BSR BE B68 8 CE SNs CES 08 BS0E8 O88 88 SENS BOSSE FON Ns SEES BESns Pesan shares Senessenss eeeaceeess 
JCORE SEG 88 GSE OS CERO PERRO RUST Th 1 SS RERET ODEO S BRORS BERET TST G OSE OS BRUGES CREEP BOE SS FORO E OES OS ee ee eee eee OSBG0 DOUSE GHS00 FEGRS FUCUC RBA GES EU DASUE FUSS ROUSE PEO OU CORU0 BOGS RDOOL OOGRRGEBeb eee 4 Hh 
ttt tated EE Oe ee LLL LLELLEOLLCLLELLLLOELELLILILELILD LIC | ph 
EE TTT TOLLE LLL“ LLL CTT TTT eT 
die ate te teat De ee ae elena) (Slane tole) oll netic el) tl ee Leen eet LCL LLL Leni LELLELL LICL LOULELLULELLLILLCLLELILCLLLLLE LIL op I EE 
SOGGSC008 SSSC5 C5008 SEGR0 C0056 COGS Re seen CEC SE Set Se Sere eee eee eee SR ES SOR SS SSRs SESS SSEEOE BESES BRes FESE8 FREES FORRES ESRE LORS EERSRE PORES ES Es RES SESE CES EE RO RRS RE SRe seees Seee eee eee eeeee eees 
a SS TOOT OO Et tH 
BSC CUeHT Fees8 Be8aa ee tele are eee ae aed ele | nee on Joe LCI Thiel itl Thi) ltd to ie Th) en TCI ILLITE LEE LEC Diet mol PT Le a H 
DOSES SREBD FEGHe SSHGL 18058 RH0LF CODRRRAOSB ee eee eee ee ee BES S8 FOES SSS SEs FE S88 CASES CORES CREE SOROS SESS C8088 SSESs POSES SS O8 ES Es BOSS OOORS RERes Fee es Sees Beene Feees Sees eeees 
JDBES SBOGS DEERE SHES SBDB0 BOLT FHOAV ODED 8850 COU00G 82008380000 00000 C0088 £4605 00080 CODDE SO0S8D0 DOSEN FSEH0 DESDE G6 800 OBG00 CO0S0 C8R00 FOF SSSC0 GSES ODES CO000 0 C886 F8D00 FUSES SBOE COSE0 Coc Bees oeesee 
EEE EEE EEE EEE EEE ——————— ee ——————————————e——————————————————————— OOO 
SBE GS BSE CHR00 FRUSE SRCRDOOHLL JOBD0 02000 FOU0N DQ808 BBOOS FOR 00 SERB EDBZD BEE Eee 
SOO CSSRF SSSS SAR ASR ASE SHOOT 14 41 DB RSE2S OOOe CCR BRS BOSSE DE O8S CARO R OES CBSE OD FRESE CREE OSES OB ORE OE OSS PRES PSOE OR REE BOG EPS OSE POSES See? CREO OSE COREE OR ERE OSES ee Oe SoeeReees Beeeeeeees 
PES RR ASREE BERS OS SSS TBSSSRHE 1+) SURO RESDE (VOD ODEDE PE LO8 BOTS DRRSROBEOD 0 LORE PORE ORE ee EO PRS BOOS OBSESSED OE UPR OOGEEE Ee oeeBoL St 
SSGe ) ROSE FOR SE BDSSR PERE TORRE REARS PORE GOSS, COREG RSENS FPSO SPORE POR EF FURS TAPE a BER ERR RRO PSPS Ie RRS CO Se CS STO SE HPT IRS OOP we eee 





IV-B-4 Analog Computer Simulation. 

The analog computor simulation of the given system is shown in 
Figure IV-B-4.1. The coefficient pot settings with associated resistances 
and capacitances for real time and real magnitude scaling are given in 
Table IV-B-4.1. 


TABLE IV-B-4.1. 


Pot Setting Associated Elements 


For an input step of a volts 
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To accurately set the back-biasing pots "a" and "a", the same pro- 


cedure as given in Section IV-A-4 was followed. The static character- 


istic of the relay is the same asthat shown in Figure IV-A-4.2, 
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The system was first excited with a 6 volt step input signal and 
was found to build up into the limit cycle shown in recorder trace 
[IV-5-4.1. The ssytem = then excited with a 12 volt step input signal 
and the system was found to go into the same limit cycle ag shown in 
recorder trace [V-B-4.2,. From these two recorder traces the limit cycle 
was determined to have an angular frequency of 5.07 rad/sec and an 
amplitude of 9.0 volts. 


The phase portrait of the system is shown in Figure [V-A-4.2. 
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TV-B-5 Comparison of the Resuits of the Four Methods for the System 


of Figure 1v-5. 


Frequency Amplitude 
rad/sec. (Volts) 


Mitrovic'’s 
Rsot Locus 
Describing Function 


Analog Simulation 
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CHAPTER V 
NONLINEARITY IN A FEEDBACK PATH 
V-A Saturation in Velocity Feedback Path. 


The system chosen to analyze is given in Figure V-A. 





6. + aA Ge 
fe 
Y 
FIGURE Y-A 
SATURATION IN VELOCITY FEEDBACK .PATH 
V-A-1 Mitrovic's Method. , 
Assigning N as the instantaneous gain of the nonlinearity, the 
characteristic equation can be written as, 
ae “ ee + a + ( 8 + 120N)S + 30 = 0. (V-A-1.1) 
Assuming the last two coefficients variable, the characteristic 
equation is rewritten as, 
a + 78° + Wee + B,S + BY = 0 (V-A-1,2) 
where 
By = 8 + 120N (V-A-1.3) 
and 
BY = 30. (V-A-1.4) 
From Appendix A, Mitrovic's equations for BS and B, are, 
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By = -LI4 wn? g, (4) + 7wne Po (4) + Wy $(%) | (V-A-1.5) 
and 
B;= 14 Wr Do (¥) + 7 wn Pal) + Wn? De (¥) (V-A-1.6) 


Substituting the values of the @ functions from Appendix A for & = 0 


gives the parametric equations of the stability curve as 


Cae Wht wy 4 Slag) 
and 
iS 7 ee (V-A-1.8) 


From equations V-A-1.3 and V-A-1.4 the M-point locus is seen to be 

a straight line segment parallel to the By axis. At the intersections 
of the M-point locus with the stability curve, equation V-A-1.4 eguals 
V-A-1.7 and equation V-A-1.3 equals equation V-A-1.8. Thus at the inter- 
sections 

Wo, 14 4W, = BO (Vases 
and 

7W),7 = S+ I20N. (V-A-1.10) 


Solving equation V-A-1.9 for Wy gives 


by, = 2.6411 


ie 11.3589, 
from which 
On) = 1.625 rad/sec 
and 


Wn, = 3.370 rad/sec. 
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Since the system is fourth order, the stability curve is a para- 
bala with its axis parallel to the BY axis. Consequently, both values 
of Wry must be considered since the M-point locus, a straight line 


parallel to the BE, axis, can intersect the parabola in two places. 


1 


For Ly, = 1,625 rad/sec, the average gain of the nonlinearity 


) 
over a cycle is found from equation V-A-1.10 to be 
N = 0.0874, 
lave 
from which the approximate amplitude of signal X in Figure V-A is deter- 
mined to be 
X = 34,325 volts. 
lmax 
For Oo,= 3.370 rad/sec, the average gain of the nonlinearity over 
a cycle is determined to be 
No ave = 0.5958, 
. from which the average amplitude of signal X is determined to be 
X =O eV OltcE 
2max 
However, the signal X is proportional to the output velocity and 
not the output angle. To obtain an approximate amplitude of the limit 
cycle in terms of output angle, Syn » one must assume the signal X is 


relatively sinusoidal. With this assumption, the amplitude of the output 


angle can be obtained from 


Ae X max 
oS CRON (V-A-1.11) 


From equation V-A-1.11, the approximate amplitudes of the output 


angle for both limit cycles are 


Se: = 3,52 volts 


and 


Oo oma = 0249 volts. 
maxX 
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The graphical presentation of the solution by Mitrovic's method 


is shown in Figure V-A-1.1. 
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FIGURE Y-A-!.] 
GRAPHICAL SOLUTION FOR SATURATION IN 
VELOCITY FEEDBACK PATH 


Mitrovic's method predicts that the limit cycle at A,,= 1,625 rad/sec 
is an unstable one and that the limit cycle at Wny.= 3.37 rad/sec is a 
stable limit cycle by the following arguments. For the stable limit cycle 
at Wry,= 3.37 the same arguments used in Chapter I, Section I-A-1 can 
be applied to Figure V-A-1.1. For the unstable limit cycle assume Ox 
is such that the first peak of signal X causes the instantaneous M-point 
to be at M The system is then stable and the DI GEEN ESE decrease. AS 


the oscillations decrease N increases moving the M-point to the right. 
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Eventually, the “average'' M-point will lie at the intersection of the 
stable limit cycle at Je 3.37. Now assume that ©; is large enough 
to cause the first peak of signal X to be large enough to place the 
instantaneous M-point at Ma. The system is unstable and the oscillations 
increase. As the oscillations increase, N decreases causing the M-point 
to move even farther to the left into the unstable region. Consequently 


the oscillations increase indefinitely. 


151 





V-A-2 Root Locus Method, 
Substituting § = JW into equation V-A-1.1 and evaluating the 


higher powers of .) gives, 


ee) 70> = 140 +J(8+120N)W+30= 0, (V-A-2.1) 


By requiring that the real and imaginary parts of equation V-A-2.1 
go to zero independently, then 
a 14 Ww +t30=0 (V-A-2,2) 
and 
-)7w> +)(8+]2AON)w=O., (V-A-2. 3) 
Dividing equation V-A-2.2 by JW the above two equations are 
exactly the same as equations V-A-1.9 and V-A-1.10 developed in Mitrovic’s 
method. The solutions of these simultaneous equations are again given 
Q) , = 1.625 rad/sec 


() »= 3.370 rad/sec 


N = 0.0874 
lave 
Ny = 0.5958 
ave 
x = 34,325 volts 
lmax 
and 
x = 5), 035) volts. 
2max 
Also applying equation V-A-1.11 gives 
Onna = 3,52 volts 
and 
Op oma = 0: 249 volts. 


The graphical solution of the root locus solution is given in Figure 
V-A-2.1. In the figure, the "open loop poles" are the poles cf the tran- 
sfer function obtained when the position feedback loop is closed on the 
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KOOT LOCUS SOLUTION For 
SATURATION IN VELOCITY ‘ 
FEEDBACK PATH 


\ 


The root locus method predicts the limit cycle at UW), = 1.625 rad/sec 
will be an:unstable limit cycle and the limit cycle at &.,= 3.37 rad/sec 
will be a stable one. The arguments of Chapter I, Section I-A-2 can be 
applied to the W, = 3.37 rad/sec crossover to determine a stable limit . 
cycle at that point. At the QO), = 1.625 rad/sec crossover, assume the 
input signal, OQ, , is large enough to cause the first peak of signal 
X to be large enough to place the instantaneous root at th The system 
is stable and the oscillations decrease. As the oscillations decrease, 


N increases causing the instantaneous closed loop root to move up the 
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locus toward the linear closed loop root. 

Eventually, the average root location over a cycle will lie at the 
stable limit cycle crossover. If ©); is large enough to cause the first 
peak of signal X to be large enough to place the instantaneous root at 
ras then the system is unstable and the oscillations increase. As the 
oscillations increase, N decreases, driving the closed loop root farther 
into the right-half plane toward the open loop root. Consequently, the 


oscillations will increase indefinitely. 
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V-A-3 Describing Function Method. 

The describing foneeren for saturation is given in Equation T-A-3.1 
of Chapter I. 

To obtain an equivalent loop transfer function, the gain of the 
error channel was first made unity by moving the gain block to the input 


channel and to the unity feedback path as shown in Figure V-A-3.1. 


Reduction of the outer position feedback thus yields 





ae Ze 
G6) = 54, 7-24 )4s74 654 30 © (V-A-3.1) 
Or + + 7 _. 20 _- ’ Oo 
15 BC AC Hoe: st Seasa 





FIGURE Y-A-3.| 


A digital computer program was used to find the roots of the de- 


nominator of equation V-A-3.1, which are, 


a 
It 


0.2168 - 41.3854 


1 
ry = 0.2168 + j 1.3854 

res 3. ee 9 {12001 
r, => 3.7168 + j 1.2011. 


Thus the equivalent open loop transfer function becomes, 


Pa OUs 


_ -A-3.2 
G(s) = (5-0,2168 £3 1.3854 )(S + 3.1768 £31,201!) wu my 


nee 





With S = JW , the magnitude and angle of equation V~A-3.2 were 
calculated using a digital computer program with (OU varying between 0.1 
rad/sec and 100 rad/sec. The results of this program with angles in the 
neighborhood of -180° are given in Table V-A-3.1. These results are plot- 


ted on the gain-phase of Figure V-A-3.2. 


TABLE V-A-3.1 
ae _ [eons] se ZeTHW aay 
| ~182.57 
-181.51 
-180.50 
“179.54 
-178.64 


“L779 


-179.02 


=N7 9a 


-179.86 
St eS. 
-180.134 
aie, 27) 


-180.41 





Values of the 1/6, locus are given in the two areas of interest in 


Table V-A-3.2. These values are calculated from Equation I-A-3.1. 
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TABLE V-A-3.2 


0707655 


0.08397 
0.08651 
0.08905 
ORIEN To, 


0.09413 





By interpolation of Figure V-A-3.2 and the tables, two limit cycles 


are seen to exist at 


U), = 1.625 rad/sec 
R, = 0.0686 
Xi may 43.7 volts 


and 
eur = 3,37 rad/sec 
Ro = 0.488 


Xo may= 6, 15-VOLES wy 
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To determine the approximate amplitude of the output angle, sinusoidal 
variations of the signal X must be assumed. Thus equation V-A-1.11 can be 
used to give, for the limit cycle at W, = 1.625 rad/sec 

‘S) 


and for the limit cycle at UW), = 3.37 rad/sec 


Cape 4.48 volts 


Q5 4 ee 0.304 volts. 
The describing function method predicts a stable limit cycle at 
GJ9 = 3.37 rad/sec and an unstable limit cycle at WJ, = 1.625 rad/sec. 
At LW, = 3.37 the limit cycle is determined to be stable by applying the 
Same reasoning given in Chapter I, Section I-A-3 to Figure V-A-3.2 using 
the signal X. At WJ, = 1.625 the limit cycle is determined to be unstable 
by the following reasoning. If the input, O, , is large enough to 


cause the first peak of signal X to be large R will be small and the in- 


stantaneous critical point will lie at say R, in Figure V-A-3.2. The 


d 
system is unstable due to a negativephase margin and the oscillations in- 
crease. As the oscillations increase, R becomes even smaller driving the 
critical point farther up the -180° axis and increasing the negative phase 
margin. Thus the system is always unstable and the oscillations increase 
indefinitely. If Oa is small enough such that the first peak of signal 

X is small enough, the instantaneous critical point will lie at say Ro of 
Figure V-A-3.2,. The system is stable due to positive phase margin and the 
oscillations decrease. As the oscillations decrease R becomes smaller 
driving the critical point down the -180° axis and increasing the positive 
phase margin, Eventually the positive phase margin begins to decrease due to 
the shape of the eu, C) LJ) curve and the system reverts to the 


stable limit cycle at QJ, = 3.37 rad/sec. 


WS, 
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V-A-4 Analog Computor Simulation. 

The analog computor simulation of the given system is shown in 
Figure V-A-4.1. The coefficient pot settings and the associated resist- 
ances and capacitances for real time and magnitude scaling are given in 
Table V-A-4.1. 

TABLE V-A-4.1 


Pot Setting Associated Elements 


For input of Ei Volts and Ke = 1.5 





R, = 0.5 Meg, Rey = 1 Meg 


R, = 0.5 Meg, Rey = 1 Meg 


Meg, Reo = 1 Meg 


Meg, Reo = 1 Meg 


Re =- 1 Meg, Coy = Lue 


Meg, Coy = | uf 


R, = 0.5 Meg, Ceo == eis 


R. = 0-5 aler, Ce9 =) lame 


Ry = 0.1 Meg, Ce, =a 


R,. = 0.1 Meg, Coy = JL ue 


Ri = 1] Meg, Res = 1 Meg 


= 0.1 Meg, Rey, = ] Meg 


13 1 Meg, Res = 1 Meg 


Ri, = 1 Meg, C ey, =e tf 


For a 3 Volt Saturation 
Voltage 
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Initial Setting 





To accurately set the back-bias pots "a", a sine wave generator 


was connected to the input of the nonlinearity simulation and the output 


monitored. Pots "a" were then adjusted until the output was limited to 
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3 volts. The static characteristic curve is the same as that in Figure 
T-A.4.2, except the saturation voltage is 3 volts instead of 5 volts. 

A step input signal of 1 volt was first impressed on the system and 
the system was found to build up into the limit cycle shown in recorder 
trace V-A-4.1. A step input signal of 5.8 volts was then impressed and 
the system was found to go into the same limit cycle as shown in recorder 
trace V-A-4.2. As expanded portion of the stable limit cycle is shown in 
recorder trace V-A-4.3. 

A visual presentation of the unstable limit cycle was difficult to 
obtain. For a step input of 6.21 volts the system was only slightly 
damped before reverting to the stable limit cycle as shown in recorder 
trace V-A-4.4. For step input of 6.22 volts the system was slightly un- 
stable with the oscillations eventually increasing indefinitely as shown in 
recorder trace V-A-4.5. 

From the recorder traces, the stable limit cycle was determined to be, 

C(t), = 3.35 rad sec 


X eee 6.2 volts 


= 0.30 volts. 
Op lmax 


The unstable limit cycle was determined by averaging recorder traces 
V-A-4.4 and V-A-4.5 with measurements taken on the first few cycles only, 
UG), = 1.618 rad/sec 
Ximax= 42.3 volts 
Oo, max= 4.29 volts. 


The phase portrait of the system is shown in Figure V-A-4.2. 
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‘Recorder Trace V-A-4.1 
Saturation in Velocity Feedback Path, 1 Volt Step Input. 
e °e 
SCALES; £-0.1 Volts/Line E& - 0.1 Volts/Line X - 0.5 Volts/Line 


Y - 0.2 Volts/Line Time - 1 Second/Division 


163 





CIE aap ; ry E : iH me === ===) 5 h- 





PER EP CEE 
eS Sriet ates fetecsctatscse war 
SEER CCEE EE EEE EE EEEE EEE 
aeceeeceeace PC EECEELEEELL 
_- SS SS =S= Sear SEER EERE ELL 


Reece 


eee lee ELL 
S665 55 255. eee R eee SSeS SeSSeaee 
VL a 
Seeley ate sa ueaieaezcascs7eseces 
I S Eaue 2a\\ (mae BGS ee eee GeanEeaaur 
SSS6s===— sss eC RE EGES==SScca= ===Slec 
(ZEEEEEEEEELT CULE LETT LE LEEEE LL 


. wee 6 ame 


CLE Pledad tHE ELE 
=3 Spe ENE EPP aia eet 








NEAL ee SSnB2! 
: = 7aee SNe eee NEsSH==seeeaSscesese 
Zoe eae/cany. SAPP PEPE CERES 

EEL NAECMECEECT pee (LEEEEEEEE 
PEEEEEELECEEE wEEEESee= mH CEEEEEEELE 


















EEEEECPER TC Ce EEE EEE 
EEL EERE ETE a DU} PTET TTT aaa 
EPar eee arante ie mentees 
= = e/egie CEA I 

SEEISEIBE FHP ot i * { He Fon ca GVEA E REFS en OE 
Ree RUN St eee ate ee ate 
= SS aReesaa Ee eS SSeS. 


Recorder Trace V-A-4,2 
Saturation in Velocity Feedback Path, 5.8 Volt Step Input. 
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SCALES:€=0.5 Volts/Line €¢ -0.5 Volts/Line X - 5 Volts/Line 


Y - 0.2 Volts/Line Time - 1 Second/Division 
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Recorder Trace V-A-4, 3 


Saturation in Velocity Feedback Path, Limit Cycle. 
» 
SCALES: & -0.02 Volts/Line & - 0.1 Volts/Line °* X - 0.5 Volts/Line 


Y - 0.2 Volts/Line Time - 0,2 Seconds/Division 
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V-A-5 Comparison of the Four Methods in Predicting Limit Cycles of the 
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V-B Saturation in Acceleration Feedback Path. 


The system chosen to analyze is shown in Figure V-B. 





On 4 E Bo 
Figure V-B 

Saturation in Acceleration Feedback Path 

V-B-l. Mitrovic's Method. 

Assigning N as the gain of the non-linearity the characteristic 

equation is: ) 
4 3 ez 
S + 7s + (14 + 10N)s” + 88 + 20 =0, (V-B-1.1) 


; 1 2 : F oe 
Assuming the s” and s° coefficients to be variable the characteristic 


equation becomes 


. + 78° + Bys° + Bys + 20 = 0 (V-B-1.2) 
where \ : 
By = 8 (V-B-1.3) 

and 
By =14+10N , (V-B-1.4) 





From Appendix A, Mitrovic’s equations for 8 


1 and E. are 


ies &-[20 @.() = Twi, (= wn Bo (4) | 


(Veb-~ j 2 5) 


B, = 


n 


sh ~204,(4)+7 ww Gal$) +n $s) | 


(V-E-1.6) 


¢ 
Substitution of the @ fenstians from Appendix A for XZ = QO gives 
the parametric equaticns of the stability curve as 
R 2 


eee 5 (V-B-1.7) 
and 





The shape o 
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easily recogni 


Lou 


must equal equation V-B-l1 


determined tea be 


Also at the intersection equati 


v4 


t 
bs 
_— 
A 


</ 


qvation V=B-1.4 


- . ) 8 ee es 
or OU) yields the 





X = 10.82 volts 
max 


To obtain an approximate amplitude of the angular oscillations, one 
must assume X to be sinusoidal. Then the amplitude of the error signal 


may be determined by 


§ = X mar 2 GoTo) 
"EO Soe 


or 


s 18.94 volts. 


The graphical solution by Mitrovic's method is sketched in Figure 


max 
V-B-l.1. 


oF 
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Mitrovic's method predicts an unstable limit cycle by the following 
reasoning. If the input signal is large enough such that the first 
peak of signal X is large enough, the instantaneous nonlinear gain will be 
small enough to place the M-point at M, in Figure V-B-1.1 the system is 
unstable and the oscillations increase. As the oscillations increase N de- 
creases moving the M-point further into the unstable region. Thus the 
average M-point can never be in the stable region and the oscillations 
increase indefinitely. If the input signal is small enough such that the 
first peak of signal X is small enough, N will be relatively large placing 
the instantaneous M-point at M in figure V-B-1.1. The system is stable 
and the oscillations decrease. As the oscillations decrease N increases 
approaching one and moving the average M-point farther into the stable 
region. Thus the system is always stable and all oscillations eventually 


die out, 
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V-B-2. Root Locus Method. 
Substitution of S$ = JW into the characteristic equation and 


evaluation of the higher powers of J yields 


wt —\7w2— (14 +ION)Ww2 +) 8W4+20=0. (V-B-1.1) 


Requiring that both the real and imaginary parts of equation V-B-1.1 
go to zero independently 
wi (14410N)w* 420 =0 (W-B=1,2) 
and 
— \7 ws +8W=0., (V-B-1. 3) 
By dividing equation V-B-1.3 by JW , the frequency of the limit cycle 
is obtained directly as 
soe— oH = 1,069 rad/sec, 
Substituting (&) into equation V-B-1.2, the average gain of the non- 
linearity is obtained as 
N = 0.4614 
ave 
from which the approximate amplitude of signal X is obtained as 
X =40.82 volts. 
max 
To get the approximate amplitude of the error signal, X is assumed 
sinusoidal and equation V-B-1.9 is used, 
& = 18.94 volts. 
max 
The root locus solution is sketched in Figure V-B-2.1. Before 
sketching the figure, the unity feedback loop was first reduced and the 


poles of the reduced transfer function were used as the poles of the 


equivalent loop transfer function. 
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The root locus predicts an unstable limit cycle by applying the same . 


reasoning given in Section V-A-2 to Figure V-B-2.1. 


LFA, 





V-B-3. Describing Function Method. 

The describing function for saturation is given in equation [-A-3.1 
of Chapter I. 

To obtain an equivalent loop transfer function the unity feedback 


loop was first reduced which gives, 


20 
SS Sue Sh 
6) St 47554 1494827485 +20 : Soa 


A digital computor program was used to find the roots of the de- 


nominator of G‘(s). 


r= 0.0869 + j 1.201 
r, = 0.0869 - j 1.201 
r, = -3.5869 + j 0.9652 
r, = - 3.5689 - j 0.9652 


Thus the equivalent open loop transfer function becomes 


GH) = Le 
S- 0.0869 +)1,201)(S +3. 5869450,9652) 
(V-B-3.2) 
With s=Jw , the magnitude and angle of equation V-B-3.2 were 
calculated by a digital computor program with QQ) varying between 0.1 
rad/sec and 100 rad/sec. The results of this program wit! angles in the 
neighborhood of -180° are given in Table V-B-3.1. These results are also 


plotted on the gain-phase plane of Figure V-B-3.1 as the linear loop trans- 


fer function. 
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TABLE V-B-3.1 


/G' HOw) deo 





The values of 1/G) in the area of interest are given in Table V-B- 


TABLE V-B-3.2 


OS 7 62 
0.4364 
0.4601 
0.4625 
0.4648 


0.4719 


0.4954 





By interpolation of Figure V-B-3.1 and the tables, the limit cycle 
iS at 


68) 1.609 rad/sec 


and 
RiavO237 4 
from which 


Eo. S87-VOlese 


PS 
it 


max 
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To obtain an approximate amplitude of the error signal, X is 
assumed sinusoidal and equation V-B-1.9 used, 
ce 23.4 volts. 


The describing function method predicts an unstable limit cycle by 


applying the same reasoning given in Section V-A-3 to Figure V-B-3.1. 
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V-B-4. Analog Computer Simulation. 

A good simulation of the acceleration signal was difficult to obtain. 
The simulation that was finally decided upon was to use the velocity signal 
and approximately differentiate it by simulating s/(s + 100). By using 
this transfer function an extra gain of 100 had to be simulated in the 
feedback path in order to maintain the D.C. gain constant of the path un- 
altered. Thus the simulated gain in the feedback path is 100 K, or a0 
instead of 0.5. 

The analog computor simulation is shown in Figure V-B-4.1. The co- 
efficient pot settings and associated resistances and capacitances for 
real time and magnitude scaling are given in Table V-B-4.1. 

The diode back-biasing pots were set in the same way described in 
Chapter I. The static characteristic curve is the same as Figure I-A-4,.2, 
An input signal of 17 volts was first applied and the system was 

found to be stable but only slightly damped as shown in recorder trace 
V-B-4.1. An input signal of 23 volts was then applied and the system was 
found to be unstable as shown in recorder trace V-B-4.2. The unstable 
limit cycle was found to occur for an input step of 19.8 volts and is 
shown in recorder trace V-B-4.3. 

From recorder trace V-B-4.3, the limit cycle is determined to have an 
amplitude of ae = 10.6 volts and a frequency of WW = 1.03 rad/sec. 

The phase portrait of the system is shown in Figure V-B-4.2. 

The simulation of the acceleration signal is admittedly poor, especial- 
ly since the amplitude of X differs considerably from that predicted by 


the describing function. 
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Table V-B-4.1 


Pot Setting Assoctated Elements 


For a step input of Ei Volts 


1 Meg, Rey 1 Meg 


1 Meg, R., = 1 Meg 


fl 


1 Meg, R 1 Meg 


£2 


1 Meg, R 1 Meg 


PZ 


1 Meg, Coy Lut 


1 Meg, Cry 1 uf 


1.0 0.5 Meg, Ceo 1 uf 


1.0 0.5 Meg, Ceo 1 uf 


1.0 0.1 Meg, c 1 uf 


£6 


0.4 O.1 Meg, Cre late 


1.0 1 Meg, C., = 0.1 uf 


1.0 = Meg, Cry, = 0.1 uf 


P20 -l Meg, Re 1 Meg 


1.0 -oMege Rea 1 Meg 


as, Meg, Rey, 1 Meg 


ae) 1 Meg, Ree 1 Meg 


1.0 1 Meg, Ces i pe 


1.0 1 Meg, Ree 1 Meg 


0.05 Initial Setting saturation of 5 volts 





Although the nonlinearity input X is not strictly sinusoidal, it is 
not too far from a sinusoid. The large discrepancy can not be accounted 
for solely by wave shape. However, if a root locus solution similar to 
that of Section V-B-2 is done on the system assuming 50ea + 100) in 


the inner feedback path the following results are obtained: 
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When the results are compared with those of the root locus solution, 
Section V-B-2, it is seen that the extra pole in the feedback path causes 
a slight increase in frequency and a significant increase in the approxi- 
mate amplitude. Consequently, whatever it was in the actual analog simu- 
lation that caused the frequency to be less (1.03 rad/sec) than the pre- 


dicted value (1.069 rad/sec) also caused the amplitude to be considerably 


less than predicted by the describing function. 
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Saturation in Acceleration Feedback Path, 17 volt Step Input 
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Saturation in Acceleration Feedback Path, 23 volt Step Input 
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Limit Cycle, Saturation in Acceleration Feedback Path, 19.8 Volts Step Input 
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V-B-5. Comparison of the Four Methods in Predicting Limit Cycles of the 


System.of Figure V-B. 


Method Amplitude (Volts) 
| i = rad/sec X-Max | 


Mitrovic’s 
Root Locus 


Describing Function 


/Anglog Simulation 
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PART II 
A STUDY OF LIMIT CYCLES OF A SYSTEM WITH 
TWO GAIN VARIABLE NONLINEARITIES USING MITROVIC'S METHOD 
Introduction 

The author had no experience in the problem of a system with two 
non-linearities and consequently had no idea of what performance to expect. 
Thus a system was arbitrarily chosen and an analog computer study of the 
system was started. 

The system chosen is shown in Figure VI. thie particular configura- 
tion was chosen primarily because each of the nonlinear gains appears alone 
in two of the coefficients of the characteristic equation, which is, 


4 + ii + ary + (8 + KK, N.)8 + KK.N, = 0 (VI-a) 





Figure VI 


System with Two Saturation Nonlinearities 
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Assuming the last two coefficients variable, the characteristic equa- 


tion becomes 


: + 7 = + ee 4+ B,s + By = 0 (VI-b) 
where 
BS = KKLN, (VI-c) 
and 
By = 8 + KK LN, (VY-d) 


The equations of the stability curve are the same as those in Chapter 


¥, Section V-A-1 and are repeated below 


i 


2 4 
BL = 14 Wa - Wy (VI-e) 


and 


This particular arrangement easily lends itself to presentation on the 


BY versus BS plane in that 


1, The stability curve is independent of the gain K since the 


. coe and a terms. 


stability curve only depends on the coefficients of the s 
2, Omce K is chosen a second rectangular grid of constant K. and 
K, can be drawn. 
3. Once K, and K, are chosen then any instantaneous or average 
M-point is easily plotted using the K, - K. grid. For example, if K = 6 
17 0.5 and No = 
0.25, then this M-point is easily plotted at the intersection of the K, mi 3 


and K, = 4 and by some means an M-point is determined by N 


line and the K. = 1 line of the secondary grid. 
In the study, three terms are used to describe certain parts of the 
system shown in Figure VI and need to be defined. These are: 
Plant - that portion of the system represented by the linear trans- 


eect ion GS) = wenyspeyseay 
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Linear System - the equivalent linear closed loop transfer 
function with both feedback loops reduced and no nonlinearities present, 


or the linear system described by the characteristic equation VI-a with 


N, and N, equal to 1.0. 


Linear M-point - that point on the BY versus By plane which is 


determined by equations VI-c and VI-d with Ni and N, equal to 1.0, 
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Chapter VI 


ANALOG COMPUTER STUDY AND STATEMENT OF RESULTS 
Besides the poles of the plant, there are six parameters which directly 
or indirectly affect the two assumed variable coefficients of the character- 


istic equation. These are Q&, K., E As a starting point 


. 51° K, K,, and Es 


t a 
the plant gain, K, was made 20 and both saturation voltages were made equal 
to 3 volts. For these conditions the analog simulation is shown in Figure 
VI-l. The settings of the coefficient yoke and their associated resistances 
and capacitances are given in Table VI-l. 

The step input signal, ES or © , was obtained by a pot reduction of a 
100 volt d.c. source through an isolating amplifier. To change K. various 


combinations of Ais Ris Bos and R, were used, however, both channels had to 


Z 


be the same. To change K. various combinations of 2143 and Ri 3 were used. 

Because the diode resistance is not zero during conduction and pot 
loading effects of pots "a" and "a''', the back bias pots were readjusted 
slightly as the amplitude of the limit cycle changed. The readjustments 
were made by eye such that on the average the output signals of both satura- 
tion simulators were 3 volts. A typical output signal of a saturation simu- 
lator for a large input signal is sketched in Figure VIi-2. 


First the stability curve was drawn on the BY versus B, plane and then 


1 
a second grid of constant K, and K. lines was drawn as shown. in Figure VI-3. 
From this figure, it is apparent that one should get a good indication of the 
limit cycles possible by taking many linear starting points located by K. and 
K.. and observing the resulting limit cycles. 


The initial study took exactly this form, that is, a K, was chosen and 


several responses to a step input were observed for various values of K. on 
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Pot Setting 


a, =k, for K, Z 1.0 
Ke 
4, * Io for K.. PZ 0 


a, = K, for K, & 1p 


ay -= for K,, > 1.0 
a, = 1.0 

a, = 1.0 

ae = 1.0 

ap = 1,0 

a, = 1.0 

ag = 1.0 

ag = 1.0 

aig = 1.0 

417 0.4 

212 = 1.0 

ai3 7 K. for K. Z 1.0 
a.3 * +. for K, > §.0 
ai, = sO 

ais = 1.0 

46 = 1.0 


a = 0.03 Initial Setting 


a‘ = 0.03 Initial Setting 


TABLE VI-1 


Associated Elements 




























1 Meg, 


O.1 Meg, 


1 Meg, 


0.1 Meg, 


1 Meg, 
1 Meg, 
1 Meg, 
1 Meg, 
1 Meg, 
0.5 Meg, 
0.5 Meg, 
= 0.1 Meg, 
O.1 Meg, 


1 Meg, 


1 Meg 


O.1 Meg, 


Meg, 





uf 





Meg, 


Meg, Meg 





volt saturation level 


volt saturation level 
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a horizontal traverse across the plane. In the study, K was varied from K, 
= 0.5 to K, = 6.0. K. was varied from K. = 0.5 to K. = 10.0, About every 


other value of Ks the input step, E,, was varied over several values to see 


i 
if the resulting limit cycle was dependent on the magnitude of the step in- 
put signal. The only limitation placed on E was overloading of any one of 
the computer amplifiers; consequently, for certain linear starting points 
oF te was as large as 200 volts. The results of these computer runs are 
tabulated in Table I of Appendix B. 

From this first portion of the analog computer study several things 
were evident. These were; 

l. Wo limit cycles were obtained from a step input signal when 
the linear system M-point was inside the stable region on the BO versus By 
plane no matter what the input signal magnitude. 

2, Limit cycles were only obtained for step input signals when the 
linear system M-point was outside the stable region and these limit cycles 
were independent of input signal magnitude and only dependent on the linear 
system M-point position. 


3, The wnstable region of the B. versus B, plane can be divided in- 


1 
to three regions such that only one of the nonlinear elements saturate or both 
elements saturate for the linear system M-point lying in a given region. 
Several runs were made for various linear M-points inside the stable 
region and for various step input signal magnitudes, From the results of 
Section V-A of Chapter V, one might expect for linear M-points close to the 
left portion of the stability curve that the system could be driven unstable 
for large step input signals. However, all linear M-points chosen showed 
some damping even for the linear point extremely close to the stability 


curve. A typical response for a linear M-point very close to the left 
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portion of the stability curve is shown in Recorder Trace Vi-l] for a 

small step signal input and Recorder Trace Vi-2 for a large step signal 
input. Typical responses for @ linear M-point very close to the right 

side of the stability curve are shown in kecorder Trace VI-4 for a large 

step signal input. Typical responses for a linear M-point in the well 

damped portion of the stability region are shown in Recorder Trace VI-5 

for a small step signal input and Recorder Trace VI-6 for a large step signal 
input. 

Although it has been stated thac no limit cycles were obserwed for 4 
linear M-point in the stable region, the system could be driven inte a limit 
cycle for repeated disturbances in the input signal after the initial step 
input. This type of response was only noted for linear M-points with 4 very 
small damping ratio and close to the left portion near the maximum of the 
stability curve. A typical response of this type is shown in Recorder Trace 
VI-7 for extremely little damping. In this trace, a limit cycle was achie- 
ved with only one disturbance after the initial step. Recorder Trace Vi-8 
shows the same type of response but with more damping. In this trace three 
disturbances were required to achieve a limit cycle. 

To obtain these responses an initial step input signal was applied. 
Then with the system in 4 very oscillatory slightly damped condition the 
disturbance was applied. Te apply the disturbance, the input pot was given 
a small quick twist to approximate a step change in inpet signal. It was 
found that the disturbances had to be (1) generally small in amplitude such 
that the error channel did not appreciably saturate, (2) applied at a defin- 
ite time during a cycle and in a specific direction. It was found that the 
maximum effect was achieved when a positive disturbance was applied when the 
error was passing through zero from minus to plus. Alsc, the maximum effec 


i) 
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was observed for a negative disturbance as the error was passing through 
zero from plus to minus. This type of system performance was only found 
late in the experimental phase of the study and was not investigated 
thoroughly. 

Recorder Traces VI-9 and VI-10 show that the resulting limit cycle for 
a linear M-point in the unstable region is independent of the magnitude of 
the step input signal. These are typical responses, and at all points where 
the effect of the step input signal was investigated, the limit cycle was the 
same for different input magnitudes. Comparison of Recorder Trace VI-9 or 
VI-10 with Recorder Trace VI-11 illustrates the effect of the relative posi- 
tion of the linear M-point and the stability curve on the resulting limit 
cycle. Recorder Trace VI-1ll shows a limit cycle of higher frequency and much 
lower amplitudes than the one in Trace VI-10. 

The most striking result of the initial computer study was that the 


unstable region of the Bo versus B, plane can be divided into three sub- 


1 
regions. These subregions are shown on the coefficient plane of Figure 
VI-3. For the linear M-point in Region I, only the error channel saturates 
and the limit cycle can be predicted from the intersection of the K, co- 
ordinate of the linear M-point and the stability curve. For the linear 
M-point in Region II both channels saturate, but the resulting limit cycle 
frequency always indicated an intersection of an average M-point locus with 
the stability curve at a point to the left of the maximum point of the 
stability curve. For the linear M-point in Region ITI, only the tacho- 
meter channel saturates and the limit cycle can be predicted from the inter- 
section of the Ke coordinate of the linear M-point and stability curve. 

The dividing line which proceeds upward and to the left from the sta- 


bility curve will be referred to as dividing line No. 1. The dividing line 
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which proceeds almost horizontally to the right from the top of the 
stability curve will be referred to as dividing line No. 2. 

The approximate position of dividing line No. 1 was obtained from the 
results tabulated in Appendix B. To fix more accurately the position of 
the dividing line, more runs were taken on the analog computer simulation 
by choosing a K. and varying Ree A point on the dividing line was then 
determined when the amplitude of signal X5 in the resulting limit cycle was 
equal to 3 volts, the saturation voltage. The results are tabulated in 
Table 2 of Appendix B. Typical changes in the limit cycle as K. was in- 
creased across No. 1 dividing line are shown in Recorder Traces VI-12 and 
VI-13. The main feature to note from these two traces ig the slight in- 


crease in the amplitude of signal X, from just under 3 volts to just over 


Z 
3 volts. 

The approximate position of dividing line No. 2 was obtained from the 
results of Table I in Appendix B. To fix more accurately the position of 
the dividing line, more computer runs were made by fixing K. and varying 
Kae The results are tabulated in Table 3 of Appendix B. Typical changes 
in the response as Ke was increased above the dividing line are shown in 
Recorder Traces VI-14 and VI-15. Recorder Trace VI-14 is a couple of 
oscillations of the limit cycle resulting from a linear M-point at K = 8.0 
and K. = 2,38. However, sufficient time was allowed to ensure that the 
system did not exhibit the same response as shown in Recorder Trace VI-1L5. 
Recorder Trace VI-15 is a series of three sections taken out of the response 
for a linear M-point at KR. = 8.0 and Ke = 2,39 to show the salient features 


of the entire response. The time intervals between the sections are in- 


dicated on the trace. Section one shows the build up to an apparently 
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stable limit cycle with only the tach channel saturating. However, closer 
inspection reveals that the amplitudes of signals X, and Xn are increasing. 
Section 2 starts at approximately the same amplitudes as occur at the end 

of Section 1, but quickly builds up to the final stable limit cycle as shown 
in Section 3 at an expanded time scale. The salient features of these two 
recorder traces are as follows: 

1. The sudden build up from an apparent limit cycle to a final 
stable limit cycle for a linear M-point just above dividing line No. 2. 

2. The large difference in signal magnitudes between the limit 
cycles of the two traces and the difference in their frequencies. 

3. The change in frequency of the oscillations during the build 
up to the final limit cycle in Recorder Trace VI-15. This point is not 
very noticeable on Trace VI-15 but it is noticeable on traces with an ex- 
panded time scale using dividers to measure the period. 

4, The frequency and amplitude of Trace VI-14 is very nearly the 
frequency and amplitude of the apparent limit cycle of Trace VI-15. 

Since the foregoing results allow one to predict the limit cycles for 
a large number of linear system M-points for 3 volt saturations, the rext 
logical step was to see if the dividing lines are independent of saturation 
level. Thus both error channel and tach channel saturation voltages were 
adjusted to 10 volts. (pots "a'' and "a'" equal to 0.1). The same analog 
computer procedures used to determine the dividing lines for Eo = E50 = 3 
volts were again used. The resulting limit cycles for the linear starting 
points are tabulated in Tables 4 and 5 of Appendix B. A comparison of the 


dividing lines for the 3 volt and 10 volt saturation levels is given below 


in Tables VI-2 and VI-3. 
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TABLE VI-2 


1 Dividing Line 


10 volt Saturation Level 3 volt Saturation Level 





TABLE VI-3 


No. 2 Dividing Line 


Level 
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Tables VI-2 and VI-3 indicate that both dividing lines are essential- 
ly independent of saturation level, although dividing line No. 1 is more 
doubtful than dividing line No. 2. In order to definitely show this 
point, more saturation voltages will have to be used and more sophisticated 
measurement procedures should be used. 

Practically all of the data given in Appendix B is graphically dis- 
played in Figure VI-4 through VI-9. Figure VI-4 is a plot of the limit 


cycle amplitude of signal X, as K. is varied with K, fixed, Figure VI-5 


2 


is a plot of the limit cycle amplitude of signal X, for like K, variations. 


1 
Figure VI-6 is a plot of the limit cycle amplitude of signal X, as K, is 
varied with KR. fixed, Figure VI-7 is a plot of the amplitude of signal 

xX, for like variations of Kae Figure VI-8 is a plot of the limit cycle 
frequency as K. is varied with K, fixed, Figure VI-9 is a plot of the 

limit cycle frequency as K, is varied with K. fixed. 

These six figures essentially show the results as already stated. By 
correlation of these figures with the location of the linear M-point in 
Figure VI-3, one can see the regions in which only one element saturates. 

For example, the upper right portion of Figure VI-8 shows Region III of 
Figure VI-3 in which one would expect the frequency to remain constant as 

K. is increased. Also, the transition across dividing line No. 2 is very ap- 
parent in all the figures because of the large discontinuity in the curve. 
Some curves for relatively small K, and K. abruptly stop and then start up 
again in a different portion of the graph. This is due to the variable 

being inside the region of stable linear M-points in Figure VI-3, thus no 


limit cycle. An important result to note which has not yet been pointed out 


rs» that for K, greater than 3 and K. greater than 3, the frequency of the 
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resulting limit cycle is essentially independent of variations in either 
Ke or Kee 

Figure VI-10 shows various instantaneous and "'average'' quantities plotted 
against the number of elapsed oscillation peaks for a linear M-point just 
above No. 2 dividing line. A better abscissa name would be "number of 
elapsed half cycles" although both are related to the elapsed time by the 
characteristics of the system. The quantities plotted in the figure are 
defined below: 


B min = The minimum instantaneous value during a half cycle 


of the as coefficient calculated by equation VI-d 


using the minimum instantaneous gain N, ... N,_. 
2min 2min 


occurs when the signal X, is maximum during the half 


Z 


cycle, 


Bi ave = The "average" value of the si coefficient over a half 


cycle as calculated from equation VI-d using the "average" 


gain of the nonlinearity N . The describing function 


2ave 


was used to relate N to the maximum value of X 


2ave 2 


during the half cycle. 


Lae = The minimum instantaneous value of the s° coefficient 


calculated by equation VI-c using the minimum instan- 


N occurs when the signal X 


taneous gain Ny nin’ aS 


1 


is maximum during the half cycle. 


Bees = The "average" value of the s° coefficient calculated 


from equation VI-d using the "average" gain of the non- 


linearity N The describing function was used to 


lave 


relate N to the maximum value of X, during the half 
lave 1 


cycle, 
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W von The average frequency of the oscillation over a half 
cycle obtained by measuring the half-period of the 
oscillations. 

Justification for using the describing function in the manner indicated 
above will be indicated in Chapter IX. 

If the “average"' M-point is assumed to indicate the stability of the 
system over the half cycle, then it is interesting to note that the waver 


age" M-point, as defined by B and Bayes moves through a portion of 


lave 
the stable region on a constant B line. It is.also noted that under limit 
cycle conditions at the extreme right side of Figure VI-10, the values of 
B, and B calculated using W ave in Mitrovic's stability curve equations 


are approximately equal to the end values of B and B ee curves which 


lave 


were calculated using the describing function and the amplitude of the limit 


cycle. A comparison is made below. 


From Stability Curve From Figure VI-10 
= (.” B, = 30.9 Bi ave = 31.4 
B= 43,2. B= 42.8 
oO oave 
Towite =< —— 





AVC let a eee) 


FIGURE YW-2 
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Recorder Trace VI-2 


to Left Side of Stability Curve, K, = 2.4 | 


Linear M-Point Close 


K. = 1.7, 62.5 Yolts Input Step 
SCALES: 


x - 10 Volts/Line Y, - 0.2 Volts/Line X, - 0.5. Volts/Line 
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Time - 5 Seconds/Division 
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Recorder Trace VI-3 


Linear M-point Close to Right Side of Stability Curve, 
K, = MET el K. = 3.9, 2.6 Volts Input Step. 


SCALES: 

X, -0.2 Volts/Line Y, -0.2 Volts/Line Y.,-0.2 Volts/ 
-0.01 -0.01 Volt s/Line Line 

X,-0.2 Volts/Line 


Time - 5 Seconds/Division 
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Recorder Trace VI-4 





Linear M-point Close to Right Side of Stability Curve, K. = Ihe, 0) 
= 3.9, 95 Volts Input Step. 


SCALES 

Xy - 5Volts/Line Y1 - 0.2 Volts/Line ¥570.2 Volts/Line 
-0.05 Volts/Line -~0.05 Volts Line ’ 

Xy - 0.3 Volts/Line Time - 5 Seconds/Division 
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Recorder Trace VI-5 


Linear M-point in Well Damped Region, K, = Ome K. = 2,0,% 5.44) voles 
Input Step. 


SCALES: : 
Ky -0.1 Volts/Line ¥,70.2 Volts/Line X,-0.1 Volts/Line 
Y,-0.2 Volts/Line Time - 5 seconds/Division 
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Recorder Trace VI-6 


Linear M-point in Well Damped Region, K. = 0.5, Ko = 28 
60 Volt Input Step. 
SCALES : 


Xi- 1 Volt/Line Y,70.2 Volt/Line X,-0.2 Volt/Line 


Y,-0.2 Volt/Line Time - 5 Seconds/Division 
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Repeated Sisturbances, M-Point at Ke = 2.3, K ~=e122 
SCALES: 


xy -0.2 Volts/Line Y,70.2 Volts/Line X,-0.5 Volts/Line 
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Repeated Disturbances, Linear M-Point at K, = 2.3, K, = 2.0 
SCALES: 
X,-0.2 Volts/Line Y,70.2 Volts/Line X,-0.5 Volts/Line 
Y,70-2 Volts/Line Time - 5 Secnnds/Division 
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Recorder Trace VI-9 
Unstable Linear M-point at K ee K = 4, One Volt Step Input 


SCALES: 
K,-0.5 Volts/Line ¥,70.2 Volts/Line X,-1 Voit/Line 
Y,70.2 Volts/Line Time 5 Seconds/Division 


-0.2 Second/Division 
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Recorder Trace VI-10 --—- 









Unstable Linear M-Point at K, = 3, KR, = 4, 10 Volt Step Input. 


SCALES: 


X = 2Volts/Line ¥,-0.2 Volts/Line X71 Voit/Line 


-0.5 Volts/Line . 
Y70.2 Volts/Line Time - 5 Seconds/Division 
-0.2 Seconds/Division 
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Recorder Trace VI-ll 


Unstable Linear M-Point Close to Stability Curve at Ke ed ie ee 


t 
SCALES: 
xy - 1 Volt/Line Y.-0.2 Volts/Line X,-0.5 Volts/Line 
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Recorder Trace VI-12 
Limit Cycle for Unstable Linear M-Point to Left of No. 1 
Dividing Line, K, = 4, K. = 0.8 
SCALES: Xx, -0.5 Volts/Line ¥,70.2 Volts/Line X,-0.2 Volts/Line 
Y,-0.2 Voits/Hine Time -0.2 Seconds/Division 
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Recorder Trace VI-13 


Limit Cycle for Unstable Linear M-Point to the Right of Neo. 1 
Dividing Line, K. oe K. = 0.82, 
SCALES : 


xX, -0.5 Volts/Line Y,70.2 Volts/Line X,-0.2 Volts/Line 


Y,70.2 Volts/Line Time-0.2 Seconds/Division 
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Recorder Trace VI-14 


Limit Cycle for Unstable Linear M-point Below No. 2 

Dividing Line at K, ng 4 See K, eG 

SCALES: 

X)-0.1 Volts/Line Y, 70.2 Volts/Line X71 Volt/Line 
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CHAPTER VII 
EAFLANATION OF THE RESULTS AND SOME METHODS OF 
PREDICTING THE RESULTS 
an attempt fs meade in this chapter to explain the system responses 
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Insight imto t*< systér can be cbtained by first looking at the system as 
only having error channel saturation and then looking at the system as only 
Waving tach channel saturation, First, assuming only error channel setura- 
tion an? applying the arguments given in Chapter I, Section [-A-l to Fig- 


wre Vi-3, it is determined for any KS 4,5 the error saturation produces a 
steble limit cycle For K, >» 4.5 the linear system can never be made stable 


by warying K 5; consequently, there is no possibility of a limit cycle with 


(6 


error charnel saturation only and x. > 4.5. This fact is easily shown by 
te 


a2 quick root locus analysis. The root locus of the inner loop is sketched 


in Figure VEI-!. For a ae large enough to place the inner losp closed loop 


roots in the xigsht half plane, no amount of F can make the root locus swing 
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ints the left half plese. Starting with the closed loop roosts of the inner 
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loop a5 the open looxp roots of the su ster loop, the assymptotes remaiz 
And -180° altheugh their cantroid is meved to the right. A rough chetk of 


the angie of emergence from the wpper pole shows that this angle is a first 


Guadrant angle, The jw axis corgsover of the inner loop root losus i¢ 


if 


determined to be K 4,5, Thus the outer loop root locus can only procede 
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@s indicated in the figare and for & 7 4.5 the system is always vwastakile, 
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Figure VII-l 
Root Locus of the Tachometer Feedback Path 
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Assuming only tach channel saturation and applying the arguments 
of Chapter V, Section V-A-1 to Figure VI-3, it is determined that for a KG 
245: 

1. The side of the stability curve to the right of the maximum 
(Wy > 2.645) yields stable limit cycles. 
2, The side of the stability curve to the left of the maximum 
(Wr, < 2.645) yields unstable limit cycles, at least for 
0.75 < KX 2.45. 
3. Which limit cycle governs the response is entirely dependent 
upon the input signal magnitude. 
For K, > 2.45 Figure VI-3 indicates the linear system is always unstable 
no matter what the Ke» thus no limit cycle is possible. For K, < 0.75, 
the absolute minimum gain of the nonlinearity, i.e., N, = 0, the M-point 
is still within the stable region, and thus an unstable limit cycle is 
impossible. 

From the above discussion, one can safely say that with both satura- 
tions present a stable limit cycle will always result no matter where the 
linear M-point lies, even though individually each saturation has certain 
areas in which the single nonlinearity system is absolutely unstable. This 
statement is true because for large K, and K, the tach channel saturation 
moves the "average" M-point to the left such that the error channel satura- 
tion can move it downward to lie on the stability curve, or vice versa. 

A result of the Chapters of Part I will also be helpful in being able 
to predict the results of Chapter VI. The magnitude of the describing func- 
tion at the intersection at the intersection with the linear loop transfer 
function curve on the gain-phase plane is essentially equal to the magni- 


tude of the "average" gain of the nonlinearity as obtained by Mitrovic's 


Zao 





and root locus methods. Table VII-1 gives a comparison of Nee and Gy 


under the limit cycle condition for all the systems analyzed in Part I. 


TABLE VII-1 


Mitrovic's 


Describing 
Function (G_) 


Weriincarity No. of Root Locus 
Zeros N N 
ave 


Saturation 


Dead Zone 


Ideal Relay 


Relay with 
Dead Zone 


Tach FB 
Saturation 


Acceleration 


0.60 
0.8051 


0.5951 
0.4089 


T,.1053 
0.5573 


0.4568 
0.6208 


0.08737 


0.5958 


0.4614 


0.60 
0.8052 


0.590 
0.4072 


1.1032 
Oe >/ 3 


0.4543 
0.6208 


0.08739 


0.5959 


0.4614 


0.604 
Oe esas 


0.5989 
0.3987 


1.0976 
0.5629 


0.4225 
0.5757 


0.08802 
0.5957 


0.4648 





FB Saturation 


The number of significant figures listed for the describing function 
in Table VII-1 is misleading since Gh was calculated using at best three 
figure accuracy from a graphical plot of an intersection. Even with this 
uncertainty in Cys the comparisons only differ greatly in the last two or 
three significant figures. With more accurate calculation of the inter- 
sections of the three methods, the numbers should be exactly equal. A 
non-rigorous proof of this equality follows: 

The describing function method applies the Nyquist criterion and is 
basically a frequency response technique involving a variable system criti- 
cal point. For this proof assume that the symbol, Gy» and the words, des- 


cribing function, are undefined in terms of the amplitude of the input 


signal to nonlinearity. Then for gain variable nonlinearities, the Nyquist 
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criterion becomes 
G6, GWW) — (VII-1) 
This general equation can then be applied to any one of the systems anal- 
yzed in Part I. 
For convenience apply the root locus criterion to the system of Section 


I-A, Chapter I which yields 


s(s + 1)(s re ae. (VIT-2) 
But in equation VII-2 
10 = G(s) 
s(s + 1)(s + 2) 
which is the linear loop transfer function. Thus 
Noe - 1 (VII-3) 


where G(s) can now be any linear loop transfer function. 

In the root locus method of predicting limit cycles for gain variable 
nonlinearities, one is only interested in the limit of stability of the 
linear loop transfer function. Another way of saying the same thing is that 
one is interested in finding the point on the linear system root locus at 
which self sustained oscillations are possible. Of course, this occurs when 
the linear root locus crosses the JW axis or when s = UJ. Sub- 
stituting \W for s in equation VII-4, 

NGC Sw) =cl. (VII-4) 

It is a well known fact that both the frequency response technique 
and the root locus accurately predict the absolute stability of a system. 
Consequently, equations VII-1 and VI-4 must be equivalent and therefore N and 
Gy must be equal since the G(JW ) is the same in both equations for any given 
system. 
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The fact that Noe and Gy are equal under limit cycle conditions is 


really only useful when G_ as a number is related to the amplitude of the 


D 
input signal to the nonlinearity. Of course this relationship is strictly 
valid only when the input signal is sinusoidal. Thus when the use of the 
describing function can be justified in a given system, root locus and 
Mitrovic'’s methods will also accurately predict the limit cycle amplitude. 
The limit cycle can then accurately and quickly be determined by substitut- 
ing s = JU) in the characteristic equation, calculating the amplitude 
from the describing function relationship. 
PROBABLE QUALITATIVE EXPLANATIONS 

1. Stable responses to step inputs for which the linear system 
M-point lies within the stability region. 

a. For small signal inputs in which the error channel does 
not Saturate, Recorder Traces VI-1, VI-3, and VI-5 are representative. 
Two of these traces show the first peak of signal Xo is large enough to 
saturate the tach channel. However, after the first velocity peak, the X, 
and Xy signals are well below the saturation levels of 3 volts and the 
oscillations die out according to the position of the linear M-point. In- 
sight into an explanation is obtained if the instantaneous minimum M-points 


are calculated from X x of the first peak in Traces VI-3 and VI-5 and then 


2ma 
plotted on the By and Versus By plane. These instantaneous M-points are shown 
as M, and M, respectively in Figure VII-2. The starting point of the arrow 


is the linear M-point and the head of the arrow is the instantaneous mini- 
mum M-point. One trace which was not included in the results of Chapter VI 


showed that the first peak of X, saturated the tach channel sufficiently to 


2 


drive the instantaneous minimum M-point into the unstable region as shown 


Zoe 





by M, in Figure VII-2. However, the tachometer channel did come out of 


3 
saturation and the system was stable. Thus for linear M-points close to 
the left portion of the stability curve the tach channel may saturate but 
the magnitude of signal X, is not large enough to cause the instantaneous 
M-point to be in the unstable region for a sufficient portion of the first 
oscillation to drive the average M-point to the stability curve. This was 
thought to be due to the relatively small values of Ke As K. is increased 


in an attempt to raise the magnitude of the first peak of signal X,, two 


a3 
effects combine to keep the response stable. First as K increases, the 
linear M-point moves to the right and the linear system damping is increased, 
at least the damping is increased until the M-point moves approximately across 
the axis of the parabolic stability curve. The increased ane ine tends to 
nullify the expected increase in Xo max’ Secondly and most important, the 
relative distance the average M-point must move to reach the left side of 

the stability curve increases, The left portion of the stability curve is the 
governing stability limit for the linear M-point in the stable region because 
only by tach saturation can the system be driven unstable. Thus, the increase 
in Xo ax is not sufficient to offset the increase in this relative distance 
such that the instantaneous M-point can lie in the unstable region for a 
sufficient portion of the first oscillation to drive the average M-point to 
the stability curve. 

As the linear M-point moves to the right past the approximate axis of the 
stability curve the damping recat to decrease but the relative distance be- 
tween linear M-point and the left’ side of the stability curve is quite large 
and obviously overrides the increase in peak X5 and keeps the system stable. 


This relative distance is in some manner related to the degree of saturation 


necessary to cause limit cycle operation since this is the graphical distance 
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the average M-point must move to reach the unstable limit cycle governing 
the tach channel. 

b. For large input steps in which the error channel in- 
itially saturates, the same explanations given above also apply after the 
error channel comes out of saturation. Recorder Traces VI-2, VI-4, and 
VI-6 are applicable. It is constructive, however, to analyze the response 
prior to the error channel becoming linear. The effects to be discussed 
are most apparent in Recorder Trace VI-4. With the error channel initially 
saturated, the initial driving voltage of the plant is 3 volts. Since Xo 
does not change instantaneously (it is proportional to velocity) the drive 
ing voltage remains at 3 volts for a short instant of time. The trace in- 
dicates that Xo increases very rapidly and that Yo essentially jumps to the 


saturation value of 3 volts. Also, the output angle and X, have not changed 


1 
during this short time interval. Thus because of the signs of the feedback 


voltages, the driving signal of the plant goes to zero after a very short 


ry 
é < 
am * 


time, This leads to the initial plant driving signal shape shown in Figure 
VII-4a as a first approximation. Second approximations might be as shown 
in Figure VII-4b and c. The solution for the output angle when the first 
signal approximation is applied to the plant is of the forn, 

GO, =C, t+ oe +Ca a Cy ee es) toyz oa 


-~2(t-t, —4 t-t 
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Figure VII-3 


M-point Loci for Stable Linear Systems for Large Step Input 
Magnitude 








Differentiation gives the output velocity as 
-(t -t,) 


° -t b wea2c ; —4C i 
O5 = C,+Cz e + Ca e. + Ge = = Cet+G (as 


-2(t-¢, ~4(t- 
-t,e i) a e 4(e-t) 


@ 
Since Xy = Ke Go , the initial portion of signal X, is easily explained 


2 
by the sum of the expontential terms given above. Time ty is the time at 
which the tach channel saturates. At time to, indicated on Recorder Trace 
VI-4, the response of the plant to the initial input has ended and the 
system is coasting because the plant driving signal is zero. At time ts 
the velocity of the coasting system has decreased such that Xo is less than 

3 volts and the tach channel becomes linear. As Xo slowly decreases below 

3 volts, a slowly increasing positive voltage is applied to the plant be- 
cause the error channel is still transferring a positive 3 volt signal. 

This is indicated on the Recorder Trace by the fact that Xo levels out at 
about 2.8 volts instead of remaining in an expontential decay. The relative- 
ly sudden decrease in Xo after time t., and before Xo levels out can not be 
explained, except by noting that the saturation simulations do not strictly 
limit their outputs to three volts. This is due to a €inite diode resist- 
ance and pot loading effects when the diode conducts. Signal yy is initial- 
ly 3.6 volts and signal Y5 is 3.1 volts after saturation. Thus the plant of 
the analog cimulation is actually being driven by about 0.5 volts or less 
during the period in which the author has assumed zero driving voltage. This 
D after time Cy 
Eventually the output drives to within 3 volts of the input and the error 


fact may account for the relatively sudden decrease in X 


channel comes out of saturation at time t)- From this time on the system 


responds according to the linear M-point. 
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Possible Waveshapes of Initial Plant Driving Signals for Stable. 
Linear M-point with Large Input Steps Applied. ; 


The same response in Recorder Trace VI-4°can also be qualitatively ex- 
plained strictly by the system M- point: motion. With the error channel in- 
itially saturated the M-point intiediately jumps downward along the linear 
M-point K, coordinate. For Recorder Trace VI-4, this initial M-point jump 
moves the M-point to K, = 0.0316 or By = 0.632 by equation VI-c. This M- 
point motion increases the system damping a small amount over the damping 
of the linear system. Within a very short time interval, the tach channel 
saturates and rapidly moves. the M-point to the left, essentially along a 
constant Ke line since the error has not appreciably changed to cause an 
foward component in the M-point motion. The rapid movement of the M-point 
to the left greatly increases the system damping, possibly close to criti- 
cal damping. With the M-point moving into a highly damped region within a 
short time, very few oscillations would be expected and would be demoed out 
very quickly by say time t. in Recorder Trace VI-4. The peak velocity is 
more difficult to explain but is certainly dependent upon how soon the tach 
channel saturates since this is the controlling factor of the time required 


for M-point motion into the highly damped region. After time t it is 


2? 


difficult to explain the response strictly by M-point motion until time t, 
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is reached when the error channel comes out of saturation. At this 

time the system is completely linear and the oscillations occur according 

to the Y and Wn of the linear M-point. These oscillations will be small 
since the disturbance as the error channel comes out of saturation is 

small. The estimated M-point motion corresponding to Recorder Trace VI-4 

is sketched in Figure VII-3 along with those of Recorder Traces VI-2 and 
VI-6. The corresponding times from recorder trace VI-4 are indicated along 
the path. The lines of éonstant Y were not calculated but their shape was 
sketched from those calculated for a different fourth order system. For the 
given fourth order system, the constant ¥% lines will not differ greatly 
from those sketched. The objective of the constant 4 lines is to show 
that the instantaneous damping can increase or decrease either a small 
amount or a large amount for both horizontal and vertical motions of the 
system M-point. The entire figure is not to scale but the numbers listed 
are correct. 

c. The fact that the system could be driven into a limit 
cycle by repeated disturbances and linear M-points close to the left portion 
of the stability curve cannot be explained satisfactorily by M-point motion. 
However, from Recorder Trace VI-8, it was noted that between disturbances 


the signals X, and X, were 


and after the first saturated peak of X 1 m 


2? 
slightly damped and below 3 volts but at a generally higher amplitude than 
before the disturbance, With little linear system damping the magnitude of 
Xo oscillations decreases very little between disturbances. After several 
disturbances the X, oscillations are just below the saturation level. The 


2 


next disturbance causes the first peak of X, to saturate N, to the degree 


2 
required to move the average M-point from the linear point to the stability 


curve and the oscillations build up since tach saturation gives an unstable 
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limit cycle. The final limit cycle is stable because the error channel 
eventually saturates as the oscillations build up. It was stated in 
Chapter VI that the magnitude of the disturbance had to be small such that 
the error channel does not saturate appreciably, otherwise the magnitude of 
Xo oscillations after the disturbance was less than the magnitude of the 
oscillations before the disturbance. This result is easily explained by 
the M-point motion. When the error channel saturates the M-point moves 
rapidly downward, greatly increasing the damping. The increased damping 
damps out the oscillations before the output has time to come into corres- 
pondence with the new input demand. Consequently, the disturbance causing 
the oscillations after the error channel comes out of saturation is much 
less than the disturbance which caused the error channel to saturate. The 
response is very similar to the response to a large step input applied at 
time zero as described previously. 

2. Limit cycles were obtained only for linear M-points in the 
unstable region. 

a. The result that the limit cycles were independent of the 
magnitude of the input signal was suspected because of the arguments given 
in the introduction to Part [. For saturation in the error channel limit 
cycles are only possible for an unstable linear system. Thus the existance 
of the limit cycle is entirely dependent on the linear system and not the 
step input signal because the linear poles and zeros determine the possibik 
ity of instability and the linear gain determines whether or not the linear 
system is stable or unstable. For saturation in the tach channel, an un- 
stable limit cycle is possible with a stable linear system. However, the 
existance of this limit cycle is only dependent on the linear system. The 
magnitude of the step input signal determines only whether or not the non- 
linear system response is unstable or stable. Thus if the existance of a 
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limit cycle for each of the nonlinearities individually is indepdent of 
the magnitude of the input step function, then the existance of a limit 
cycle with both nonlinearities in the same system is independent of the step 
input. | 

b. The result that the limit cycle amplitudes of Xy and Xo 
are dependent on the relative distance between the linear M-point and the 
stability curve is explained by degrees of saturation. Degree of satura- 
tion means the relative magnitudes of NF a and eve which are required 
to place the average M-point on the stability curve. For systems with 
linear M-points close to the stability curve, the average M-point has a 
small distance to travel either horizontally and/or vertically to reach the 
stability curve. For sinusoidal input waveshapes of signals Xy and Xo» Lt 


has been shown that Nae. can be related to the amplitude of the input signal 


through the describing function. For small average M-point motions N 


lave 

and No ave will be just less than 1.0 from equations VI-c and VI-d because 

B and B will be very close to the linear M-point values. By the 
oave lave 


describing function relation for values close to one, the input amplitude is 
just slightly greater than the saturation voltage. Thus the amplitudes of 


Mm and X 


1 2 will be slightly greater than 3 volts. For linear M-points re- 


latively distant from the stability curve, the values of B and B 
oave lave 


will be greatly different from the linear M-point coordinates. Thus Niave 


and N are small fractions and through the describing function relation- 


2ave 


ship the amplitudes of X, and Xo are correspondingly large. 


i 
c. The result that for K. y 3 and K, », 3 the frequency of 

the limit cycle is essentially independent of KR. and K. variations can be 

explained only in terms of other results which are apparent in Figures VI-4 


through VI-7. From Figure VI-5 it is noted that for constant K, above the 
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stability curve and for K, > 2, the amplitude of is independent of K, 
variations and only dependent on Ke variations. Thus as K increases above 
K. = 2, the vertical distance the average M-point must move to reach the 
stability curve is relatively constant for K, constant. This fact can be ex- 
plained as follows. No matter what the K, or Kee the initial driving 
voltage of the plant remains at 3 volts, thus the initial changes in @o 
and Sa are unaffected by K, and Kee Also ©o is subject to an extra 
time lag which Bo is not. Thus for K > K, > 3 , it can 

safely be said that the tach channel saturates before the error channel and 
the resultant average M-point motion is to the left instead of downward. 
Consequently, it is difficult to concieve of a limit cycle occurring on the 
right side of the stability curve. This fact is shown in experimental re- 
sults since all the limit cycles obtained from linear M-points in Region IT 
exhibited frequencies on the léft side of the stability curve. Thus the 
tach channel first saturates to move the average M-point to the left past 
the maximum of the stability curve along a constant K, line. Once the 
average M-point is in this position the vertical distance it must move to 
reach the stability curve by error channel saturation is relatively con- 
stant as long as the Kecoordinate of the linear M-point is kept constant. 
As K, increases, Figure VI-7 shows that for K, > 3, the limit cycle ampli- 
tude increases linearly with K, and that K, variations have little effect 


when K > 2. The linear variations of X i with K is reasonable since 


1lma 
it has been shown that the relative vertical distance the average M-point 
must move is independent of K, for a given Ko Thus the nae coordinate of 
the average M-point of the limit cycle is relatively independent of the 


variations in the linear M-point for K and K. greater than 3. Since the 


Bye coordinate is always the same, it can intersect the left portion of 
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the stability curve at only one point and therefore the frequency of any 
limit cycle resulting from a linear M-point located by K, 7 3 and K. ae 


is unaffected by variations of K, and Ke These same conclusions could be 


obtained by using X and Figures VI-4 and VI-6. However, it must be 


2max 
pointed out that the limit cycle frequency and the signal amplitudes 
are not necessarily independent in Region II of Figure VI-3, thus the above 
explanation is somewhat trivial. In Regions I and III, the frequency is 
independent of amplitudes and the above result indicates that it is also 
essentially true in Region II for K, 7 3 and K, me Le 
3. Division of the Unstable Region into Subregions. 
The signals Xy and X, are not independent of one another. 


That is, except for constants Xy is the time derivative of Xy- Consequent- 


ly it is entirely reasonable for the system in limit cycle operation that 
the limit cycle amplitude of one of the signals is not large enough to 
allow the limit cycle amplitude of the other signal to be great enough to 
saturate the other channel. For example, a limit cycle with small velocity 
excursions will not drive the output angle far enough before the velocity 


changes sign to allow signal X, to be greater than the saturation level of 


sl 


the error channel. This is especially apparent in Recorder Trace VI-14 
where the limit cycle amplitude of Xo is 14,3 volts whereas the amplitude of 


Xy is only 1.61 volts. 


a. The explanation of dividing line No. 1 is this. Record- 
er Traces VI-12 and VI-13 are applicable. Recorder Trace VI-12 shows the 
limit cycle resulting from a linear M-point just to the left of dividing 


line No. 1 at K. = 4 and K. = 0.8. The amplitude of X, is 8.1 volts and 


the amplitude of Xy is 2.95 volts. Recorder Trace VI-13 shows the limit 


cycle resulting from a linear M-point just to the right of dividing line 
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No. 1 at K. = 4 and K. = 0.82. The amplitude of X, is unchanged at 8.1 


1 


volts within the readable accuracy of the trace whereas the amplitude of 


X, has increased to 3.02 volts. In this case the amplitude of the velocity 


and output angle probably did not change enough to be readable had these 

° a 
signals been recorded because — = 2,95/0.8 = 3.69 volts and 2 = Sys) 
0.82 = 3.69 volts. However, the increase in K. was enough to cause Xo max 
to be greater than 3 volts. 

b. Dividing line No. 2 is due to an entirely different pheno- 
menon since from Figure VI-7 for K, Z 2.4 the amplitude of Xy is well below 
saturation at about 1.5 or 1.6 volts. Also, since the maximum of the stabil- 
ity curve is at K, = 2.45 and the amplitude of Xy is about half the satura- 
tion voltage at K, = 2.38, one would expect something quite drastic to occur 
as K is increased above 2.45 because the error channel must saturate for 
larger values of Kos The phenomenon which occurs can be explained as fol- 
lows. For a linear M-point just below No. 2 dividing line, R2corder Trace 
VI-14, only the tach channel saturates and the limit cycle average M-point 
lies on the right side of the stability curve. A quick calculation using 


the amplitude X, and equation VI-d shows the minimum instantaneous M-point 


2 


of the limit cycle is at B = 41.5 and K, = 2.38 where as an unstable 


lmin 


tach channel limit cycle exists at B, = 41.0 for K. = 2,38. Thus the mini- 


1 
mum instantaneous M-point never crosses into the left unstable region and a 
dynamic equilibrium condition is established about an average M-point on 
the right side of the stability curve at the K, = 2.38 line. As Ko is in- 
creased to place the linear M-point above No. 2 dividing line, two effects 
add together to cause the response and limit cycle of Recorder Trace VI-15. 
These two effects are sketched in Figure VII-5 along with B amin of Recorder 
Trace VI-14,. The first effect is that the increase in K, increases the re- 


lative distance the average M-point must travel to reach the stability curve 
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from Aato fa‘. As indicated previously, this causes the limit cycle 
amplitude of X5 to increase slightly. Secondly, the increase in K, causes 
the graphical distance between a stable tach channel limit cycle and an 
unstable tach channel limit cycle to decrease from MAb to S50 

This change in Qb is easily calculated from Equation VI-e. For K = 
2.4, Equation VI-e gives 


yn = Ae ~48=6 


from which = ip 

Wry = 

——— 
on = 8 
a 
The smaller of the values is the desired one. Substituting for UW in 
Equation VI-f gives 
By = 42, 

fonk = 2.38, Equation VI-e becomes 


LW, Se ipers ute © 
or Don = 5, Sis 

n= 3,188, 
Substituting the smaller value in Equation VI-f gives 

By = 40.705 
Since a parabola is symetric about its axis the change in A b is then 
Change in Ab = 2(42 - 40.705) = 2.59 

It can be seen from the figure that the instantaneous minimum M-point for 
K. = 2.4 would lie in the unstable region even if the amplitude of X, did 
not change by the first effect. Consequently, the small increase in K, 
causes the instantaneous minimum M-point to move past the left portion of 
the stability curve into an unstable region. This in turn causes the dynamic 
equilibrium that would be attained ©2::1¢ at point 1 in Figure VII-5 to be 
upset and the signals increase in amplitude. As the signal X, increases 
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the average M-point moves through the stable region along a constant K, 


line. Eventually the error channel saturates allowing the system to 
reach a new limit cycle described by an average M-point on the stability 


curve somewhere to the left and below point 2 in Figure VII-5. 
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: ; ( Ke=2.38 
O Pn Bama 
A 
v0) | 4b AQ 
) L-PORTIONS OF STABILITY 2 
an | CURVE PARABOLA = 
1 3 
a gs Figure VII-5 B, ov Ke 


Sketch of Phenomenon Occurring at No. 2 Dividing Line 


| Another indication that the above explanation correctly describes the 


phenomenon is taken from Figure VI-4. It was noted that the constant K 
lines are linear with varying K. for M-points below No. 2 dividing line. 


The line for K, = 2.4 was extrapolated to K. = 8 which should indicate 


Xo max if the limit cycle occurs on the right side of the stability curve. 
If this Xo nax is used to obtain a Bi min from equation VI-d, one gets 
Bees en ene e 
Imin X 
~ 2max 
If this Bimin is now plotted on the coefficient plane, it falls just to 
the left of the left side of the stability curve which has a B, coordin- 


1 
ate of 42 for kK. =z 2.4, Although this Bimin is further to the left than 


‘would be expected it could be considerably in error because the extrapolated 
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point is far from the last experimental point and the lower portion of 
the K, = 2.4 line in Figure VI-4 is based on only four experimental points. 
The above explanation either greatly conflicts with certain concepts 
already used or the author has failed to see the points that would resolve 
the conflicts. For instance, the concept that the average M-point over a 
cycle or half cycle indicates the stability of the system over the period 
or half period which the average M-point is valid. The average M-point has 
to move through the stable region to reach the left side of the stability 
curve since the linear value of K, is the upper vertical limit any M-point 
can attain. Average M-point motion through the stable region indicates 
that the signals should decrease slightly in amplitude before increasing to 
their final limit cycle values. But all of the Recorder Traces exhibiting 
the same response as Trace VI-15 showed increasing amplitudes until the 
final limit cycle amplitude was reached. None of the traces showed even a 
slight decrease from one peak to the next. One explanation that might 
resolve the conflict is that the fictitous average M-point jumps from the 
right side to the left side of the stability curve. This is not very reason- 
able either since the average M-point must in Some way be related to the in- 
stantaneous amplitude of signal X, and this does not change appreciably from 
peak to peak. Also the Biave curve of Figure VI-10 indicates a smooth trans- 
ition across the stable region, although the use of describing function when 
the system is not in a limit cycle is questionable. 
c. Dividing lines independent of the magnitude of the satura- 
tion voltages. 
In Region I just to the left of dividing line No. 1, 
only the error channel saturates. For linear M-points in this region the 
frequency of the limit cycle is completely determined by the stability curve 


which is independent of the saturation voltage. For linear M-points in 
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Region I, the M-point locus is always the vertical line which is the K. 
coordinate of the linear M-point which again is independent of the satura- 
tion voltages. Thus the limit cycle is completely determined by the inter- 
section of a curve which is dependent on W,yand a straight line which is 
only dependent on the linear M-point. From the intersection, the fre- 


quency of the limit cycle and a quantity called Nees are determined. Only 


when N is related to X does the saturation voltage play a part. 
lave lmax 


Thus the relative distance the average M-point must move to reach the 


stability curve is unaltered by E. as long as the linear M-point remains 


1 


unchanged and in Region I. As shown in the statement of results and in 


the explanation of No. 1 dividing line, the dividing line occurs when MD ae 


equals E52 and K is increased from low values. Whatever the exact re- 


lmax ae coma 


unaltered as long as Eel = E 52° 


lationship between X ” in the limit cycle is, it should remain 
Consequently, No. 1 dividing line should 
be independent of the saturation voltages as long as they remain equal. 

The same arguments can be applied in explaining why dividing line No. 
2 is apparently independent of the saturation voltages. In explaining the 
phenomenon that occurs at dividing line No. 2 it was stated that in the pro 
cess of building up toward a stable tach channel limit cycle on the right 


side of the stability curve, the amplitude of X, became large enough with 


2 
increasing K to drive the instantaneous minimum M-point past the left por- 


tion of the stability curve. The pertinent relative distances in this ex- 


that is distances Aa + &b and Aa 


planation are independent of E52. 
in Figure VII-5. No ave as a number is completely determined by the dis- 

A ; ; ‘ ae : 
tance a. No. is defined as E o9/Xomax and a specific No nin Ls defined 


by an infinite number of combinations of E52 and X it is notedethat 


2max’ 


N is defined exactly the same as the quantity R in the describing func- 


2min 


tion of Equation [-A-3.1l in Chapter I. Thus once No ave is determined for 
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a limit cycle the N or R can be uniquely determined by the describing 


2min 


function. The dividing line is then determined when Nomin as found in the 


above manner is equal to the N required to place the instantaneous 


2min 
minimum M-point on the left side of the stability curve. Since both Nonin © 


are independent of E. and only on the characteristics of the linear sys- 


2 
tem, No. 2 dividing line must be independent of E 52° 
PROPOSED METHODS OF PREDICTING THE RESULTS OF CHAPTER VI 
1. Linear M-point Inside Stable Region. 
a. Small step inputs which do not initially saturate the error 
channel. 

In order to predict the response of the nonlinear system 
under these conditions, the initial peak of Op for the linear system 
must be found. This can always be found by factoring the total characteris- 
tic equation with N,and N, equal to 1.0 and then solving the linear differ- 


1 2 


ential equation for a step of Esl or less. Certain linear system predic- 


tion aids may be available such that the actual solution is unnecessary. 


Then calculate X from . 
2max -K 6 
7max ae: Omax 
If X is greater than E then calculate B, , from N, .. If B,., 
2max s2 lmin 2min lmin 


plots to the left of the left side of the stability or on it, a limit cycle 


could result and the linear system should be redesigned. If plots to 


Bac 
the right of the left side of the stability curve then the response of the 
nonlinear system is essentially that of the linear system. 
b. A large step input such that the error channel initially satur- 
ates. 
The first velocity peak again must be found. This could easily 


be done by using one of the shapes of Figure VII-4 as applied to the plant 
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if the time t) could be determined. However, more study of the problem is 
required to determine reasonable values of tie Since time ty is indefinite 
an alternative is to solve the differential equation by a two part piece- 
wise linear solution. In this method assume a step input of magnitude 

Ea is applied to the plant with the linear tach channel feedback included. 
The differential equation obtained applies up to time t, Or when Xo = E 52° 
This differential equation must be solved for a and its derivatives 
up to the point of the previous criteria. The values of Oo and its 
derivatives at this time are then used as the initial conditions in the 
second differential equation which is obtained from the plant alone. The 
plant driving voltage is zero when the second differential equation applies. 


With the first velocity peak known B can be calculated and plotted on 


lmin 


the coefficient plane. If Ba min plots to the left of the left side of the 


Stability curve a limit cycle may be possible. If B ymin plots in the 
stable region the response will be similar to those of Recorder Traces VI-2, 
VI-4, and VI-6. 

A second approximate method to find the first velocity peak is now 
proposed. The author has not had time to think through the method but it 
does show possibilities. The method basically defines a new equivalent 
instantaneous transfer function which should be applicable, approximately 
at least, through the first velocity peak. First the lines of constant ¥ 
with Wr as a parameter must be calculated from Mitrovic's equations of 
Appendix A and then plotted on the coefficient plane. Knowing the step in- 


put magnitude, the initial degree of error channel saturation can be cal- 


culated from 








Knowing N,, the initial M-point which governs the system can be calculated 


li 
by equation VI-c, the other coordinate is the K. of the linear M-point 
since Ny = 1.0 initially. Plotting this M-point, the g and Wy which 
initially govern the response can be found from the constant g lines, 


This 4 and WJ, are for the dominant closed loop roots only. To find 


the remaining closed loop roots the instantaneous characteristic equation, 


4 3 2 
s + 7s + 14s” + (8 + KK, )s + KKN = 0) 


and the quadratic factor, 


2 2. 


s +2 4 Wy s+ Wy 
must be formed and the first divided by the second. Thus all of the roots 
of an equivalent instantaneous closed loop transfer function are known. 
The gain of the equivalent closed loop transfer function is just KK, Nj. 
This equivalent transfer function is shown in Figure VII-6. Two assumptions 
must now be made in order to use the initial equivalent transfer function. 
These are: 


1. The nonlinear error channel gain N,, cannot change apprecia- 


li 
bly until the first velocity peak is reached. 

2. The nonlinear tach channel gain No cannot change appreciably. 
The first assumption is shown to be approximately true in the recorder 
traces. The second assumption is not true but must be made in order to 
proceed, The primary question now becomes what value of E. to use with 


the equivalent transfer function of Figure VII-5. The author has no answer. 


A guess might be 
E,’ =E,—KE 
i i e sl 
since the initial error channel saturation has already been taken care of by 


formulating the transfer function. 
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Figure VII-6 


Initial Instantaneous Equivalent Closed Loop Transfer Function 
for Stable Linear M-point with Large Input Steps Applied. 


The reason for including the above method was to illustrate an idea 
which may or may not be applicable to the two nonlinearity problem. onl 
‘further study and development will show if the idea has merit. Certainly 
the instantaneous equivalent closed loop transfer function is easily ob- 
tained from Mitrovic's method. 

c. The author was aRanie to predict when the system oem 
be driven into a limit'cycle by repeated disturbances. However, it was 
noted that for any reasonable damping the system should remain stable. One 
trace was taken for K. = 2 and K. = 1.5, which is still quite close to the 
left side of the seuetitey curve. The “4 at this point was roughly 
Y = 0.1. The disturbances applied were commidered quite drastic. ‘They 
consisted of 6 to 10 disturbances alternating in Rien on succeeding half 
cycles, The tach channel saturated to quite a high degree but the system 
did remain stable, | 
| 2. Prediction of Dividing Lines. | - 

a. No. 1 Dividing Line. 

The criterion governing No. 1 dividing line stems from the 


\ 


relationship that must be true between signals X, and X, under steady state 


limit cycle operation. This relationship is 


Kt 
X(t) = xX 


d 
e dt NSD, 
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The only quantity available from Mitrovic's or the root locus methods 
which is related to the input signal to the nonlinearity is Nes Without 
knowing the waveshapes of x, and Xo» 


less for prediction. If sinusoidal waveshapes can be assumed, the above 


the above equation is practically use- 


equation can be stated in a useful form. That is, on an amplitude basis 


Kt 
so ee OSes | 


For sinusoids, Xa and Xo occur a quarter cycle apart but this is im- 


material in the dividing line prediction because of the manner in which 
the dividing lines are defined. From the amplitude relationship, the 


criterion for No. 1 dividing line is easily formed as 


Kt 
Ke Ww Aim a ED 


Consequently, it is apparent that an accurate relationship between Xan and 


N is needed. If the relationship used in Part I, N 


lave / 


Xan? is 


used, the amplitudes obtained differ greatly from the observed amplitudes. 


lave ~ aol 


However, using this relationship in the criterion, the dividing line is 
predicted to be a vertical line starting from 4) = 2.0 on the stability 
curve, although the use of UW) in the criterion then becomes questionable. 
In a search for a more accurate relationship, the equivalency of No and 


ve 
Gp in Part I was noted. Also the use of the describing function is indicat- 
ed by the fact that sinusoidal signals had to be assumed to put the criter- 


ion into a useable form. Thus the criterion for No. 1 dividing line can be 


written as 
Kt -l 
Ke W) Sp, OS ae ae, 
| 
where C1 1D = N and the symbol Coy (xX, ) is the inverse function of 


lave 
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: ; ; A : 
of Gy aa? or just X 1 m (Fp To find the point on the stability curve 
where the dividing line starts the same criterion is used but written in a 
different way, that is 


Kt QO Esl -—£ 


K s2 
e 


Eel replaces Xy 


close to the stability curve, the degree of saturation required to move 


= by the following reasoning. For linear M-points very 


the average M-point to the stability curve is very small. Thus Renae is 


essentially equal to E. It is then just a matter of trial and error to 


1° 
find the point at which Ko Ke and W satisfy the criterion. It should 
be noted that for this particular point on the dividing line waveshapes are 
unimportant if Ut) is considered the average limit cycle for frequency 
obtained from the average period. 

In applying the criterion to calculate the dividing line, the follow- 
ing procedure was used. A linear M-point was chosen and error channel satura- 
tion only was assumed. Fromthe intersection of the KR. coordinate of the 


linear M-point and the stability curve the (4) and N of the resulting 


lave 


limit cycle were read. A table of Xa on and the resulting Gy was calculated 


by a digital computer program. Entering the table with Nave? the limit 


cycle amplitude X, was obtained. Using ® and X 


along with the linear 
Im m 


1 


M-point coordinates in the criterion, a value of X, was obtained. If this 


2m 


K. was increased and a new X calculated 


value of X was less than E Om 


2m 2 


using the same procedure. This procedure was continued until E52 was 


bracketed by X The procedure rapidly converges to give the K for the 


2m 


chosen K. which results in Xo = E 5° Table VI-1 gives a comparison of the 
calculated and experimental linear M-points which define dividing line No. 1. 


The results of Table VI-1 were considered quite good since the analog computer 
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TABLE VI-1 


K, a ee 
t t 

Zeek 0.92 

253 0.89 

Zs5 0.85 

a0 0.82 

4.0 Oro 

5.0 0.80 

6.0 0.785 





results are subject to 3-5% error. 

The criterion can be put into a form in which once choosing a linear 
M-point Ro» the linear M-point K of the dividing line can be found direct- 
ly but the algebra is quite involved due to getting Xs Pp? from equation tom 
I-A-3.1 of Part I. The inverse sine and square root terms of Equation I-A- 
3.1 can be expanded in an infinite series form. If the two infinite series 
are added, reversion of the series sum can be accomplished to give Xan as 
a function of Gps At the intersection of the linear M-point Re coordinate 
and the stability curve, N, assumed 1.0, Equation VI-c equals Equation VI-e 
and Equation VIi-d equals VI-f. From these equations expressions for 


and Nove? which equals G_, can be obtained in terms of K, and Ree Substi- 


D? 
tution for the various quantities results in a statement of the criterion 
in terms of K, and K. of the linear M-point only. The expression is very 
complicated and includes a power series whose argument contains both K, and 
K.. 


t 
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b. Prediction of No. 2 Dividing Line. 
From Figure VI-5 and the explanation of the phenomenon 
occurring at the dividing line, the criterion governing the dividing line 
can be stated as follows: 


"For the linear M-point on the dividing line, B , as 
calculated for a stable tach channel limit 4 pa the 
right side of the stability curve must equal the B, co- 
ordinate of the stability curve on the left side where 
the linear M-point K, coordinate intersects." 


Calculation of the dividing line is again a trial and error procedure. One 


such procedure using By coordinates is as follows. Choose a K. and a Kos 


Reading the B, coordinate of the K line intersections from the graph was 


st 


found to be too inaccurate to give adequate results. Thus the By co- 


ordinates of the intersections of the K line with the stability curve had 


to be calculated from Equations VI-e and VI-f. From Equation VI-e, Q)y, 


is found by solving, 


2 
14 OD», - wo, 4 =s5 = KK, 


O 


Substitution of the two values of Boye into Equation VI-f gives the de- 


sired By coordinates. The following symbols will be used: 


Bie = The coordinate of the right intersection at which a 
stable tach channel limit cycle exists, 
Bat = The coordinate of the left intersection at which an 


unstable tach channel limit cycle exists. 


Using B N 


3 can be calculated from Equation VI-d. From N 
lr’ 2ave a 


ve 


= Gy and using the pre-calculated table, the value of N can be deter- 


2mi 


mined, It is noted that the definition of N / ke 


2min’ Noma a Bg: nae 


exactly the definition of the quantity R used in the describing function in 


Part I, Equation I[-A-3.1. Having No min? B amin can be calculated using 
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Equation VI-d. A comparison of B and B,. can be made, If 


Imin LL Le 


ee 


greater than B the linear M-point was below the dividing line and K, 


LL* 


must be increased. If is less than B the linear M-point was 


Li 
above the dividing line and K, must be decreased, Table VII-2 gives a 


Fine 


comparison of the experimental dividing line and the calculated dividing 


line using the procedure outlined above. 


TABLE VII-2 


Calculated K, Experimental K, 





The last significant figure shown in the calculated K, was interpolated 


between the differences existing in the comparison of B with B 


Imin 1L 

for the calculations at K, = 2,42 and Ka = 2.43. The results are surpris- 
ingly good considering the calculated K is subject to slide rule accuracy 
and considering the analog results are inherently subject to the 3-5% ac- 


curacy. The worst percentage difference between calculated K. and experi- 


mental K is less than 47%. 
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CHAPTER VIII 
SOME PREDICTIONS OF RESULTS FOR DIFFERENT SYSTEMS 
WITH TWO GAIN VARIABLE NONLINEARITIES 

An attempt is made in this chapter to extend certain of the ideas 
developed in study of the system of Figure VI with equal saturation voltages. 
However, it must be pointed out that the extension of the ideas is mostly 
speculation with no specific experimental results to back up the ideas. 

An attempt will first be made to indicate the expected results when 
different parameters of the system of Figure VI are changed. In Chapter 
VII it has already been shown that both dividing lines are independent of 
the saturation voltages as long as they remain equal. It is then reasonable 
to expect the same types of responses for linear M-points in the stable re- 
gion as the equal saturation voltages are changed. The argument applied is 
that the relative distances between a given linear M-point and the stability 
curve remain unaffected by the saturation voltages, thus the degree of sSatura- 
tion required is unchanged. Only the amplitudes of the signals are changed 
to cause the required degrees of saturation. 

For the same linear system with E52 > Esl certain results seem evi- 
dent. When the linear M-point lies in the stable region and for small in- 
put steps in which the error channel does not initially saturate, it should 
become increasingly difficult to find a linear M-point at which the tach 
channel saturates as the difference between E52 and E. is increased. A 


1 


point should be reached at which E52 exceeds E. by such a large amount 


1 
that the response remains linear throughout. By the same reasoning that 


the tach channel saturates less and less as the difference between E62 and 


Esl is increased, it should be more difficult to drive the system into a 
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limit cycle by repeated disturbances, For large signal inputs such that 
the error channel initially saturates, it should be found that the entire 
response approaches error channel saturation only as the difference between 


E and E. is increased. 


s2 l 


For E52 > Ee. the dividing lines in the unstable regions should be 
predictable by the criterion as stated in Chapter VII. For dividing line 


No. 1, the criterion can be rewritten as 





Kt Ud = Es2 
K X 
6 lm 


Assuming only error channel saturation, X,_ is independent of Ee since 


lm Z 


the vertical distance from linear M-point to the stability curve is indepen- 
dent of E 2° Thus for a given Ko» one would expect the dividing line to 


shift to the right where KR. and (Q) are larger as E62 is increased and 


Eo is kept constant. A minimum of calculations show that for Eel = 3 and 


Eo = 5, the dividing line starts at (j),= 2.42 on the stability curve by 


applying the criteria in the form 


NE 


K 
e 


Another point for the same saturation voltages was calculated for K, = 5 
and was found to be about K. = 1.4, The criterion in the form useful at the 


stability curve, indicates that the difference of E . and E. is unimportant 


2 1 


to the placement of the dividing line. It does, however, indicate that the 
ratio of Ee? to Ee is the quantity which determines the placement of divid- 
line No. 1. It has already been shown in Chapter VII that for equal satura- 


tion voltages, a ratio of 1.0, the dividing line is independent of their 


magnitude. 
This same reasoning also applies to the case where Eo < Esl only 
the dividing line would shift to the left of the position for Fol = E 52° 
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Thus for any given linear plant and any saturation voltages, dividing 

line No. 1 should be predictable. The only restriction is that the plant 
and the feedback paths have the same configuration as Figure VI. Postulat- 
ing even farther, if the linear plant can be normalized, which has already 
been done for 4 third order Pilsner then a master set of curves should be 


possible defining dividing line No. 1 in terms of the ratio of E 52 to E. 


w 


1 
For the system of Figure VI, described by the coefficient plane of Figure 


” to Eel becomes large 


enough such that the starting point of dividing line No. 1 moves to the right 


VI-3, it is unclear what happens when the ratio of E. 


of the stability curve maximum. At the maximum for E = 3, the criteria 


1 


gives E 2 = 6.64 or a ratio of 2.213. For this ratio and K = 5, another 
point on the dividing line is at K. = 1.72. The criterion applied at the 
stability curve keeps giving reasonable values of E52 for points to the 
right of the maximum. For example, for Esl = 3 and Wpy= 3, the required 


E is ll volts. 


s2 


It was shown in Chapter VII that dividing line No. 2 was independent 
of the saturation voltages. Although it was not stated explicitly the 


criterion for the dividing line is only dependent on N and N, .. as de- 
2ave 2min 


grees of saturation. Since the dividing line is independent of the error 


channel and it has already been shown to be independent of E . in Chapter 


Z 
VII, it is therefore the same for any combination of saturation voltages. 


Since No. 2 dividing line depends only on the degree of saturation in 


becomes greater than 


the tach channel, the error channel may saturate as E52 


lvhaler, G. J. and Brown, R. G., Analysis and Design of Feedback Control 
Systems, McGraw-Hill, 2nd Ed., 1960: 380-383. 
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E It was implied for equal saturation that No. 2 dividing line 


sl° 
separated linear M-points with tach channel saturation from linear M-points 
with both channels saturating. A better description of the dividing line 
No. 2 is that it separates linear M-points for which limit cycles occur on 
the right side of the stability curve with the tach saturation producing 
the stable limit cycle from linear M-points for which limit cycles occur 

on the left side of the stability curve with the error channel producing 
the stable limit cycle. Thus for linear M-points below No. 2 dividing line 
it is entirely possible to get both channels saturating for E> Es 
although the error channel saturation is incidental to the stable limit 
cycle. Consequently, a third dividing line must be postulated, The des- 
cription of the third dividing line would be that it separates linear M- 
points for which only the tach channel saturates from linear M-points for 
which both channels saturate. For a linear M-point in the region between 
No. 2 dividing line and No. 3 dividing line the limit cycle will occur on 
the right side of the stability curve with the stable limit cycle produced 
solely by the tach channel. The criterion for No. 3 dividing line has to 
be the same as that for No. 1 dividing line only applied in reverse. In 
general terms the limit cycle relationship between xX, and X, is 


2 
He Koi seme) 


Assuming sinusoidal waveshapes of KX, (t) and Xo(t), then 


K 
Ud) Ke Xom = Aim 





Then forming the dividing line criterion yields, 


Ke 
W Ke 


At the stability curve the criterion would simply be 


Ke f- . 
LD K+ Fs2 = Es 





Xom= Es, 
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Since dividing line No. 1 and dividing line No. 3 are obtained by differ- 
ent ways of applying the same criterion, it is reasonable to say they are 


the same dividing line at least for an E52/E. ratio large enough for the 


1 
right hand portion to be below No. 2 dividing line. With this assumed true, 
for a few points on the right side stability curve the E52 voltage required 
to just make X = E ,. = 3 were calculated by the criterion. Then a few 
lmax sl 

points on the dividing line at three of these voltages were calculated at 

K. = 8, A very rough sketch of what the complete set of curves might look 
like is given in Figure VIII-l. The points actually calculated are circled. 
The dividing lines are labeled in terms of the ratio E oo/E.y> but were cal- 


culated assuming E = 3 and E. necessary to satisfy the criterion at 


1 2 


the stability curve. The dividing lines with no calculated points other 
than at the stability curve are estimated from the shape of the ones with 
calculated points. The shape of the dividing lines from the left side of 


the stability curve were estimated by the shape for E .o/E. = 1.0 although 


1 
the value of Eo /Esy was calculated by the criterion at the stability curve. 


If the above postulation of dividing lines proves correct, then for a divid- 


= 8.77 the only meaningful portion of dividing line 


ing line such as E oo/ Eo 


No. 2 is the portion to the right of the 8.77 dividing line. The shape of 
the dividing lines close to stability curve has not been determined. Fig- 
ure VI{I-2 shows three possibilities which may occur. If these postulated 
dividing lines can be confirmed by analog computer results, then it should 
be interesting to see if the dividing lines for E oo /E.y Li 75233 Stop at 
the stability curve. For oe Ee ? 2.33 the dividing lines at worst have 
a discontinuity at the stability curve. Thus it may be reasonable that the 


other dividing lines may enter the stable region. If this is true then they 
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Figure VIIT-1 
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Figure VIII-2 


Possible Shapes of the Dividing Lines at the Stability 


Curve for Large E ,/E , Ratios. 
s2° sl 
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should have more than a casual association with possible limit cycle opera- 
tion of systems whose M-point lies in the stable region. In the small 
amount of investigation done in driving a stable linear system into a limit 
cycle with repeated disturbances, all the points investigated were close 

to the upper part of the stability curve and not much lower than K, =e2 20 

or Ly = 2.0. This may indicate that the dividing line for E/E 6} = Laie 
may circle through the stable region and leave at dividing line No. 2. The 
author just did not have the time to see if the dividing line criterion 
could give answers for M-points in the stable region. Certainly from the 
concepts of degrees of saturation and relative distances it should be easier 
to drive the instantaneous M-point past the unstable tach channel limit cycle 


for Eoo/E, < 1.0 on the first velocity peak. 


1 
When the plant gain K is changed, not very much should change as far as 
Figure VI-3 is concerned. The stability curve does not change since it is 
independent of K. The scales on the K7R, grid will vary inversely with K. 
That is, if K is doubled, the scales will be halved. The all important re- 
lative distances between linear M-point and the stability curve will not 
change as K increases, they will only have different K, and K. coordinates. 
Since the relative distances determine Nave? then No. 1 and No. 2 dividing 


lines are unchanged as far as the BoB coordinate system is concerned. If 


1 
the above extension of the dividing lines proves to be valid, then they will 
be unchanged. 

The author can concieve of no reason why most of the concepts discussed 
in connection with the plant of Figure VI can not be applied with different 
plants in the block except for second order systems. Mitrovic's method 


does not lend itself easily to the analysis of second order systems. Fora 


third order plant the stability curve is a straight line and thus dividing 
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line No. 2 could not exist. Also, with a third order system, the filter- 
ing may be insufficient such that the signal X, is not sinusoidal in which 
case the criterion for dividing line No. 1 would give poor predictions since 
it makes use of the describing function. 

With different types of gain variable nonlinearities in the blocks 
the concept of dividing lines and a knowledge of the type of limit cycles 
produced by each nonlinearity individually should give a good indication of 
the actual responses. For example, assume two elements with dead zone in 
the nonlinear blocks of Figure VI. In this case the entire stability curve 
represents unstable limit cycles for the error channel dead zone. For the 
dead zone in the tach channel, the right side of the stability curve represents 
unstable limit cycles and the left side represents stable limit cycles. It 
cannot be definitely concluded that a system with a linear M-point in the 
unstable region will always go into a stable limit cycle. For instance, 
consider a linear M-point well above the stability curve. If the average 


over a half cycle is initially large, then N can 


value of signal X lave 


1 
never be smail enough to move the average M-point downward below the maxi- 
mum of the stability curve where the tach channel Mave can act to move the 
average M-point to the left side of the stability curve. At this same linear 
M-point, it should be possible to obtain a stable response for small input 


Signals, that is, N and N never get large enough to move the in- 


lmax 2max 
stantaneous M-point outside the stable region. Thus it can be seen that 
the analysis is much more complicated than for two saturations primarily 
because the initial M-point starts at the point (0,0) of the K-K, grid. 
However, with more study, one should be able to relate certain M-point 
motions to the position of the linear M-point and thereby be able to 


describe certain dividing lines. If dividing lines can be found, then 
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their position should be dependent only on a ratio of the two amounts of 
dead zone because Na and Nw would only be dependent on the appropriate 
relative distances between linear M-point and stability curve. 

In all of the previous discussion, a reliance on the describing func- 
tion is strongly implied. It should be point out at this point that in 
general the describing function and the requirement of sinusoidal wave- 
shapes are not necessary. The only necessary quantities are the relative 
distances involved in the coefficient plane since these distances can be 
related to various phenomena occurring in the experimental results. The 
relative distances can be described and measured entirely by the somewhat 
mebulous quantities Nes and Nain or Nnax’ The author has only been able 
to give ads and N in meaningful interpretations for prediction purposes 
through the use of the describing function. It should be stressed again 
that SO ye is just a symbol which is used for lack of a better symbol. Even 
for a sinusoidal signal, a straight integral time average of the instantan- 
eous gain, N , over a cycle does not corrolate with the observed maximum 
Signals. However, for sinusoidal signals the describing function does corro- 
late re with maximum observed signal. Consequently, the relation between 
the amplitude of the input signal and ae is much more complicated than a 
time average of the instantaneous gain and for non sinusoidal waveshapes 
this relation must certainly be dependent on the waveshape. If a relation- 
ship between N ag and the input signal can be found which is independent 


of waveshape, then all of the prediction methods should be valid. 
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CHAPTER IX 
MISCELLANEOUS IDEAS AND CONCEPTS 

In studying the two nonlinearity problems, the author had a few 
ideas which were not particularly concerned with the problem or which did 
not give pertinent results. These ideas are considered to have possibil- 
ities but were not developed. 

The first idea to be discussed concerns a general method for cal- 
culating the limit cycles in any system with any number of nonlinearities 
if any limit cycles exist in the system. If no limit cycle exists, then 
the method should give no non-trival solutions. The only restrictions on 
the method are these: 

1. Each nonlinearity must be such that it can be described by 
a describing function. 

2, The linear system must provide enough filtering between the 
nonlinearities such that the describing function description can be justi- 
fied. 

The basis of the method is application of the stability criterion to 
the total nonlinear system. In Chapter VII it was indicated that the crit- 
erion for the root locus method degenerates into the Niquist criterion at 
the limit of stability or when s = JW) . In applying the root locus 
criterion to the total nonlinear system, each nonlinearity must somehow be 
included as a magnitude and angle. At the stability limit or under limit 
cycle conditions, the describing function provides the description of the 
nonlinearity as a magnitude and angle. In other words, the describing func- 
tion can be thought of as linearizing the nonlinearity on the average over 


a cycle. Thus the root locus criterion can be formed for the total system 
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including the nonlinear effects by the magnitude and angle of their 
describing function. Since the root locus criterion is just a method 
of rearranging the total characteristic equation, one may as well start 
with the characteristic equation. It then becomes convenient to express 
the describing functions as a real part and an imaginary part. 

For illustration purposes consider the system configuration of Figure 
VI. Also let the plant be fourth order, although any order plant can be 
used, and let the two nonlinear elements be any general nonlinear elements. 
That is, the nonlinear elements can be gain variable only, hysteretic, or 
frequency sensitive. The describing function of the nonlinearities can be 
written as 

Wy 2 YO = NO 2D) + Ny sg OR WD 
and 
Ny pm? WO) = Noy Ko W) + ND» W) 

The total characteristic equation is then 


4 3 2 
48 + a,8 + as + (a, + KK_No)s + KKN, 0 


substitution of the above describing functions for Ny and No 
Og 3, 14 a2 5? +Qg Sagi Gy + K Ke] N2- (Xam w) +J Naz (Xan ,w)_] € 5 
+ KKe | N,-(X1 mW) + JN, (Xim ; ) J = O 


At the stability limit for the total system, s = JW , thus 
Agu) SO.4u)®— Oy We + J WIA) + KKe[ Nar Nong w) +N 2k Kayla) Te 
+ KKe [Nie (Xm, 20) 4 Nic (Kim, &) | =9Ou 


Thus two independent equations can be obtained from the above equation 


each containing the three variables, (), Xi and Xo 
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Clearly a third equation must be obtained and this equation must come 
from the physical configuration of the system and must relate Xan to Xo 
For other than gain variable nonlinearities, the phase difference must be 
maintained between Xi and Xo Thus for the configuration of Figure VI, 


the relationship is 


: AMG IE a x, (s) 
e 


or for limit cycle operation, 


JW a Manne eon 


Substitution of this relationship into the characteristic equation will 
give one equation in two unknowns which can be broken up into two indepen- 
dent equations in the same two unknowns. Of course the nature of the des- 
cribing functions will make the algebra very messy especially manipulating 
the terms containing J into real and imaginary terms. Conceivably this 
can be done and then the single equation can be broken up into the two indep- 
endent equations, the solutions of which will not be simple. 

Thus, in theory anyway, the method is not restricted to the number of 
nonlinearities as long as a sufficient number of relationships between the 
individual nonlinearity input signals can be determined from the physical 
configuration of the system. The only limitation is the justification of 
the use of the describing function for each nonlinearity, which in itself 
will eliminate a vast number of multiple nonlinearity system configurations. 

A second idea came about in trying to predict the limit cycle for the 
configuration of Figure VI for a linear M-point in Region II where both 
channels saturate. The limit cycle must surely be described by the inter- 
section of an average M-point locus and the stability curve. The problem 


becomes one of predicting the average M-point locus from the linear M-point 
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to the stability curve. It has been demonstrated that in limit cycle condi- 
tions the amplitude of the nonlinearity input signal can adequately be relat- 
ed to ie through the describing function. Consequently, it was reasoned 
that in a highly oscillatory condition other than the limit cycle condi- 
tion, the system can be thought of as being in a quasi-limit cycle at 

least over a half cycle of oscillation. Thus, it appeared logical that 

the describing function could be used to relate the average M-point to the 
maximum amplitude over the half cycle in question. The envelope of the 
oscillations then becomes important since this determines amplitude of the 
oscillations. It was then theorized that the total envelope could be ap- 
proximated by a series of exponential segments with each segment dependent 
only on the real part of the average root over the half cycle as determined 
by the average M-point of the half cycle. 

For linear M-points in the unstable region, the highly oscillatory 
condition is assured since for the system studied a stable limit cycle al- 
ways results. Thus the following method was evolved. For a chosen linear 
M-point, the resultant closed loop transfer function can be formed and its 
denominator factored. The linear differential equation is easily solved 
by operational methods to give the initial linear response. An impulse 
forcing function of such a magnitude was assumed such that the linear enve- 


lope of signals X, and Xo essentially retained symmetry about zero. The 


1 


amplitude of successive peaks of both signals X, and X, were calculated 


1 
until a peak of one of the signals was greater than the saturation voltage. 
When a signal peak was found that exceeded the saturation voltage, the 

average M-point was calculated for that half cycle by using the peak magni- 


tude in the describing function to determine the ae for the saturated 


channel, This new average M-point was plotted on the coefficient plane. 
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In order to obtain the information necessary to calculate the next exponent - 
ial segment of the envelope, a curvilinear grid of constant % and UJ, lines 
had to be superimposed on the coefficient plane. This curvilinear grid was 
calculated from Mitrovic's mapping equations of Appendix A for negative 
values of 4, From the position of the average M-point with respect the 

n and Wy, grid, a new average 4 and Wy was read which were assumed to 
govern the response over the next half cycle. Thus the envelope equa- 


tions for signals X, and X, for the next half cycle can be written as 


1 2 
sie 8 spi 
and 
-Y4 Wnt 
X, = Xana 
The quantities X,. and X,. are the values of the envelope calculated from 


li 2i 
the linear solution,or the end points of each preceding exponential seg- 
ment. The average frequency over a half cycle changes as the governing 
% and W, of the half cycle change as 
W= WrVi- $? 

From this relation a new period, T ve? can be obtained. Since the next 
peak of the signal under consideration occurs a half cycle later, the 
magnitude of the new envelope segment at the end of the next half cycle 
was calculated by setting t = Teac in the envelope equations. Thus new 
values of Nina se and Xo max were determined which were used in the describ- 
ing function to determine a new M-point to be used over the next half cycie. 
Thus by successive calculations the average M-point locus was constructed. 

Several average M-point loci were calculated in this manner but none 
of them ever crossed the stability curve. A typical locus calculated in 


this manner is shown in Figure IX-1l. As can be seen the locus approaches 


but never reaches the stability curve. It appears to continue indefinitely 
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past the experimental limit cycle while never quite reaching the stability 
curve. The behavior of the locus calculated in the manner described is 
reasonable since the average M-point can approach the stability cutve as 
close as one pleases but the Z to be used in the next half cycle is 
always negative causing the next peak to be slightly greater than the 
previous peak. The method as outlined does not take into account the 
phase difference between the signals. One improvement would be to take 
the phase difference into account possibly by using quarter cycle incre- 
ments instead of half cycle increments. 

An observation made in the application of the method for the several 
M-points tried is that the manner of making the transition from the linear 
M-point to the first average M-point does not affect the locus appreciably 
for succeeding average M-points. That is, it apparently made no difference 
whether the exact linear M-point equation was solved completely to calcu- 
late the first peak greater than the saturation voltage or that the first 


-¥uw,Tt 


peak greater than the saturation was calculated by X =e where 
the linear system “4 and Wr are used. Thus an initial amplitude of 
the dominant channel signal greater than saturation could be assumed as 
long as this amplitude does not cause saturation in the other channel 
through the relationship 


Re 4. X, Ct) = Xa) 


Another observation to be noted was that after both channels saturate, 

the average locus departs from a horizontal or vertical line as the case 

may be. If a straight line is drawn through the next couple of calculated 
points, this straight line intersects the stability curve quite close to 

the experimental limit cycle. This observation is indicated in Figure IX-1l. 
This observation is probably coincidental, but some interpretation may be 


made in the future. 
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APPENDIX A 


General Equations of Mitrovic's methods 
If the overall system characteristic equation is 
n n= n-2 
COnS tay..> + Oqso 2 eee poe Cene me 


then Mitrovic's general mapping equations required in this thesis are: 


1. For coefficients a_ and ay assumed variable, 


Bo= — lhe wn, () +s wn? bol4) +~- -- 7 
+n. Wy 'Dy a (8) +n ey Py) | 


and 
Bi = Ar wn, (4) +43 wn7G,(4) + -- - - 
Sea ne Dna (4) +A, wy", (4) 
2. For coefficients ay and a, assumed variable, 
alne = Gy Wy “$e ~< cm Wr eer 
A-l 
and 


B2 = 2 | Ao as AW? Pp (4)+G 4 Wr Ba (4) + Ses 


+O y-) Wy 7 ,,-2 (4) tay Wn’ Py-1 (4) Z, 


The @ functions in the above equations are various functions of % 


only which occur in the derivation.of Mitrovic's mapping equations 


The @ 
functions are as follows: 
b, (4) =0 
OUES = 
f, (4) = 28 : 
p, (4) = |- 4% 
3 
dy (4) = -44484% 
$5($) = -1+/247-16% 
| 
| 
; | 
Thaler, G. J. and Ohta, T., Mitrovic's Method - Some Fundamental 
Techniques, United States Naval Postgraduate School, Research Paper, No. 
9°, Jan 1964. 
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From the above listing, a general formula evolves from which next ¢@ 
function can be formed by knowing the two preceding @ functions, This 


formula is: 


bn (¥) = —~| 24Qn-,(4) oF Dine Sra! i Cae 


For convenience, the values of the @ functions for Ub = 0 are given 
Been 
p, (0) = -1 
Qr(0) = 0 


(0) = ual 
@, (0) = 0 


P;-(0) = ao 
D,(0) = 
$, (0) =+1 


} 
>) 


| 
| 
| 


2 
Thaler, G. J. and Brown, R. G., Analysis and Design of Feedback 
Control Systems, McGraw-Hill, 1960, Ch 10. 
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APPENDIX B 


Tabulated Results of Analog Computer Study with Ey = E 52 = 3 Volts. 


Table {I gives the results of the first set of computer runs travers-~ 


ing the Bo versus B, plane with a linear system M-point. 


1 


TABLE I 


inear System Resulting Limit Cycle 
M-Point 
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Stable| for E. up |to 75 volt 
Stable|for E, up|to 75 volts 
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Oo Oo OO Ww Oo CO UT) tne Go =) 
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Oo nwnuu & we HS Ceo eS 
sg =f o ce © ¢- 2 o © © 
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One Gy2 0.96 

0.5 6.25 1.98 

O.7 teeo 2.6 

0.2 4.85 2.88 

0.9 4.6 Sale 

1.0 44 3.4 

1.5 B96 4.85 

Z20 S702 6.45 

Cae 3276 8.0 

3E0 Sa 2 a 

332 Oe et 

4.0 Se 

a20 306 0.1887 
6.0 Sie 0.1555 
G0) 3.68 0.119 
0.0 3.68 





| 


0.0937 


—— 
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Table I (Cont'd 









Linear System | Resulting Limit Cycle 
M-Point 
i Nain 

















0.2 Lia 1.0 
Dees 7.65 1.0 
0.7 6.4 1.0 
0.8 6.3 0.974 
1.0 5.7 0.815 
1.5 5.5 0.54 
3.0 2.0 5.45 0.413 
ous tc25 0.333 
a5 Se 0.278 
bea, 5.25 One 
4.0 555 0.211 
5.0 B35 0.164 
6.0 533 0.139 
8.0 5.25 0.1095 
10.) S22 0.0847 
3.3 0.216 
0.1445 
0.1065 
0.836 
0.444 
4.0 0,216 
0.1055 
4.5 


OO. Oo Ore 
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Linear System Resulting Limit Cycle 
M-Point 
c [| Sone] Minin | Yemen 
O5e27 
0.42 
6.0 0,211 


0.1405 
ele 
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Table YI gives the resulting limit cycle for linear M-point close 


to dividing line No. 1. 


TABLE ITI 


Linear System Resulting Limit Cycle 
M- point 
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fable Lil gives the resulting limit cycles for linear M-points close 
to dividing line No. 2. 


TABLE ITl1 


| Resulting Limit Cycle 
X X 







Linear System 
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